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Abstract 


The dual issues of modal decomposition for tonal sound fields and the 
temporal coherence of the modal amplitudes are investigated for the case of the 
central Arctic sound channel at very low frequencies (15-80 Hz). A detailed study of 
the Arctic modal structure for these frequencies reveals the central role played by 
the strong Arctic surface duct. The performance of each of four different modal 
beamforming algorithms when applied to the vertical array deploved during the 
FRAM IV Arctic Acoustic Experiment is analyzed. A multiple beam (or decoupled 
beam) least squares processor produces the most acceptable results for Arctic 
conditions. The modal decomposition is sensitive to vertical array tilt caused by. 
hydrodynamic drag; a technique for its estimation from the acoustic data 18 


developed. 


Tonal data taken from both the horizontal and vertical arrays deploved 
during FRAM IV is analyzed. Horizontal array results confirm the modal 
amplitudes generated from vertical array data. The rough surface scattering from 
the 1ce canopy places an upper limit of 40 Hz on efficient surface duct propagation. 
Attenuation measurements for the first mode show excellent agreement with 
predictions made for ice scattering using the method of small perturbations aud 
experimental ice statistics. The high levels of coherence observed (0.95 to 0.99) 
show that tonal signal propagation in the Arctic channel is essentially deterministic 
for time periods well in excess of one hour. The various modes may then be 


considered to maintain a constant phase relationship over time. 
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Chapter 1 


Introduction 


It is the spatial structure in a wave field that carries the great majority of 
what might be called the immediate information about the surrounding 
environment. Exploitation of this structure can provide answers to questions of 
object existence, location, identification, and arrangement. Making instrumentation 
sensitive to the spatial structure of a wave field 1s accomplished by providing the 
detection system with a directional response. While it is possible to build detectors 
having inherent directionality, it is often easier to synthesize the desired directional 
response through the employment of arrays of simpler detectors. The outputs of 
these detectors are then combined, usually in а linear fashion, through a process 


known as beamforming. 


The art and science of beamforming has a rich history of application. The 
human body employs the concept by incorporating pairs of both eves and ears. Tle 
fundamental physics of the technique is similar to that of such diverse areas of 
science as diffraction grating theory, lens optics, x-ray crystallography, and radio 
antennae design. Тһе underlying principles of beamforming are applicable 
whenever one is dealing with either the directional transmission of energy or the 


directional reception of energy propagating ın a Wave field. 


This thesis deals with a fairly new and unique method of applying the 
principles of beamforming in the complex vertical structure exhibited by thre low 
frequency acoustic field found in the world’s oceans. Since the electromagnetic 
spectrum encounters unacceptable levels of attenuation In seawater, the acoustic 


spectrum Is the wave field of clroice for all oceanic sciences. ln particular, the very 
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lowest frequency acoustic waves (those below a few hundred Hertz) are capable of 
propagating efficiently over hundreds of kilometers in water, making them a prime 
candidate for use in long range ocean surveillance and communication applications. 
An understanding of the beamforming techniques appropriate for use with this wave 
field is a central engineering issue in all but the simplest efforts to use it both 


scientifically and practically. 


1.1 Motivation 


Traditional acoustic beamforming theories typically consider the detection (or 
generation) of plane waves in an unbounded fluid medium. The primary properties 
of the medium that impact sound transmission (its sound speed and density) are 
typically assumed to be constant throughout the medium. The constant medium 
assumption is made to avoid an overly complicated propagation problem, while the 
plane wave assumption is attractive for two reasons. First, it is physically realistic 
in the case of a true unbounded uniform medium, since the spherical wave 
generated by a point source may be considered to be locally planar at long ranges 
from the source. Second, it produces mathematically tractable results, since, once 
the assumption is made, the beamforming problem can be interpreted in terms of 
spatial Fourier transforms of the observed sound field. A large body of theory and 
experience involving Fourier techniques can then be borrowed from other 
disciplines. 

The underwater acoustic research of the last two decades has increasingly 
pointed towards the conclusion that an unbounded uniform medium is a poor choice 
of model for the world’s oceans. This is particularly true at the low frequencies of 
interest here, where even the deepest ocean depths may correspond to only a few 


hundred acoustic wavelengths. This conclusion is well understood within the sound 
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propagation community, where more appropriate layered waveguide models have 
gained wide acceptance. However, its impact on beamforming theory has been 
considerably more limited. This is primarily because the incorporation of a non- 
uniform medium in beamforming theory makes the plane wave assumption 
physically unrealistic; it is the propagating normal modes that are the fundamental 
spatial sound structures in a waveguide. Giving up the plane wave assumption in 
turn requires that the Fourier transform interpretation and its attendant 
mathematical elegance be abandoned. Instead, the traditional techniques have 
generally been extended through the assumption of a medium that is locally uniform 
in the vicinity of the array. The plane wave assumption may then be maintained: 
the resultant beamformer output is interpreted in terms of the spatial Fourter 


transform of the sound field present at the array. 


This approach has three advantages. First, the procedure is reasonably 
robust. Second, it is well understood, from both the theoretical and application 
viewpoints. Third, it still works well for horizontal arrays in oceanic waveguides, 
since, for this particular geometry, the vertical structure of the normal modes 15 
effectively hidden; for a horizontal array, each mode appears to be no more than a 


plane wave with a particular grazing angle. 


The technique also has a number of serious drawbacks. First, it Is not very 
insightful, since, as mentioned earlier, the modal decomposition is much more 
physically relevant in a. waveguide than the Fourier decomposition. Second, It Is 
not a particularly efficient estimation. scheme, for it requires a large number of 
parameters (the amplitudes and phases of the incoming plane waves from all 
possible directions) to be estimated їп order to characterize the total sound field at 
the array. Representations that minimize the number of parameters needed to 


completely describe the field generally make better estimation techniques because 
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each parameter's estimate includes at least a small amount of error. Finally. and 
most importantly, the approach obscures the impact that the waveguide model can 


have on the issue of target parameter estimation. 


Consider the information about a point source that an array can extract from 
the sound field in a uniform unbounded medium. The direction to the source (its 
three-dimensional bearing) is easily found by estimating the direction of the signal's 
wavenumber vector, but an estimate of target range can only be made through the 
measurement of the curvature of the spherical wave. At anything other than short 
ranges from the source, this implies an impractically long array, so that range 
estimation is generally not considered feasible. This limitation on target range 
estimation is closely linked to the unbounded medium assumption; if one adopts a 
lavered waveguide model, though, direct estimation of the source range (and depth) 
from the observed sound field is at least theoretically possible for an array of finite 
aperture deployed vertically across the waveguide. Source range information can be 
obtained most directly from the relative phases of the various propagating modes, 
while the source depth can be extracted from the relative modal amplitudes. Other, 
less direct (and perhaps more robust), target range and depth estimation techniques 


are also possible. 


The great potential value of source range and depth estimation serves as the 
motivation for studying the nature of the modal amplitudes and phases. One of the 
fundamental scientific issues that must be addressed in assessing the practical utility 
of range estimation in a Waveguide is whether or not the relative modal amplitudes 
and phases are temporally stable. If both types of parameter can be considered to 
be constant over reasonable lengths of time, then there ts some chance that range 
estimation techinques might be feasible; conversely, if one or both show significant 


random behavior, then the chance for praetieal success ts small Note that it is the 
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stability of relative modal amplitudes and phases that is of interest here; it is 
entirely possible for these parameters to be essentially constant even if tlre 
equivalent absolute measurements display a great deal of random behavior. All that 


is required is that the absolute measurements not be independently random. 


The ultimate purpose of this thesis is to study the temporal stability of 
relative modal amplitudes and phases tn one specific instance; that of long range 
sound propagation tn the central Arctic Ocean. There are a number of reasons for 
choosing to study the Arctic channel. First, a reasonably extensive data set, 
including vertical array data, is available from the FRAM IV Arctic Acoustic 
Experiment. Second, the unique nature of the Arctic simplifies investigation In 
some important respects, although it complicates the effort in others. Finally, the 
Arctic channel possesses characteristics that make it a prime candidate, for 


successful application of the range and depth estimation techniques of interest. 


Intimately related to any such attempt is, of course, the ability to estimate the 
modal amplitudes and phases directly, which is equivalent to the development. of 
beamformers more appropriate to the waveguide nature of the low frequency ocean 
acoustic channel. In addition to supporting the ultimate objective of the thesis, 
these are of obvious interest in their own right. An ancillary purpose of this thesis 
is to characterize some of the different ways in which the modal decomposition 
night be implemented and to study the performance of these methods in the Arctic 


sound channel. 
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Figure 1-1: Signal Replication for Hydrophones in a Horizontal Array 
| (0 to 80 Hz Band) 


1.2 Preliminary Concepts 


Consider the set of time series displayed in Figure 1-1. These traces are the 
signals received on various hydrophones of a horizontal array from an explosive 
source several hundred kilometers away. The important observation to make ts 
that the signal received on any sensor can be considered to be Just a time shifted 
replica of the signal received on any other sensor. Most linear beamforming 
techniques, at some point in their derivation, assume that the effect demonstrated 


in Figure 1-1 is true, Whether or not it actually holds in practice. 


Because traditional beamforming techniques are so dependent on the shifted 


replica assumption, the methods experience difficulty when employed in situations 
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where the the assumption is not valid, as is the case with the set of time series 
shown in Figure 1-2. "These traces are the signals received from the same explosive 
source as in the previous figure, but now taken from various elements in a vertical 
array, so that the elements are distributed above and below each other in the water 
column rather than to either side, as was the case in Figure 1-1. Vertical arravs are 
of natural interest in a waveguide, since they sample the modal structure of the 
sound propagation much more fully than do horizontal arrays. Figure 1-2. for 
example, clearly highlights the modal structure of the channel much better than 
Figure 1-1. In particular, the long coda for the signal is the arrival of the first 
mode, which travels with a slower group velocity and with more dispersion than the 
other modes. The gradual shortening of the coda with depth demonstrates quite 
dramatically the variation of the shape of the first mode with frequency. Since the 
propagation path for the signals of Figures 1-1 and 1-2 is the same central Arctic 
channel examined throughout this thesis, the reader will find the actual shape of the 


first mode for various frequencies in Figure 4-4. 


If it is the modal propagation structure that tis to be studied, as it is here, then 
the shifted replica assumption is clearly inappropriate. A review of the channel 
model assumed by traditional beamforming techniques provides a better 
understanding of the alternatives. This model may be partially described in the 
vernacular of signal processing as the known signal in noise model. This ineans 
that the sound energy of interest, called the signal, is considered to be a known, 
deterministic waveform, with some unknown but non-random parameters that do 
not vary over time. The sound еу not of interest is called the no/se, and is 
assumed to be an additive random process. A tonal signal of known operating 
frequency typically has two unknown parameters, its amplitude and its absolute 


phase: these can be incorporated into a single unknown complex amplitude if the 
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Figure 1-2: The Effeet of Modal Structure on the Signals Received by 
Different Hydrophones in a Vertical Array (0 to NO Hz Band) 
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complex exponential notation Is used. If one first quadrature demodulates the 
received signal at the nominal tonal frequency and then performs a low pass 
filtering operation, then the resultant output for any single hydrophone may be 
written as 


p(t) = A + n(t), (1.1) 


where p(t) is the complex demodulated form of the received signal, A is the 
unknown complex amplitude, and n(t) is the complex noise process that results from 
processing the noise process received at the hvdrophone through the demodulator 


and filter. 


This deseription tells only part of the full story, however, since it includes no 
information about the spatial nature of the field; that is, it does not indicate how 
the signal at one hydrophone is related to the signal at some other hydrophone 
located nearby. The shifted replica assumption is used to define this relationship. 
The model of equation (1.1) can then to be extended to cover all № sensors in the 
array. Grouping all the different demodulated hydrophone outputs into one V X 1 


complex vector, and their noise processes into another, one may write 


p(!) = ЕА + п(/), (152) 


where the .V X I complex steering vector E indicates the phase shift that each 


replica undergoes in propagating to its particular hydrophone 


Е = С И) 


When the sensors ın the array are distributed at varvınz depths in a 





S 
waveguide, as they are in Figure 1-2, the modal nature of the propagation must be 
acknowledged, making the simple shifted replica assumption an inadequate 
description. Instead, each normal mode of the total field can be considered to be an 
independent waveform, so that the the shifted replica assumption must be applied 
on a mode by mode basis rather than to the full signal. A vector of unknown 
complex amplitudes now exist, one for each of the M modes included in the model. 
Note that the selection of the number of modes M to consider is essentially a 


modeling decision that must be made by the user. The mathematical signal 


description E now becomes an N X Af complex steering matrix of the form 


- zs Jð 1 = 18 12 о • o | 2 : ie ) A 
ale : jb. Pal =) у ju. 9 • e 21/21) x "i i 
= ee е” ee 
Е = ' | | (1.4) 


where the AP column of the matrix represents the generalized steering vector for 
the 4% mode, which now includes amplitude factors 2(2;) that reflect the size of the 
mode shape at the various sensor depths in addition to the modal phase shifts 2. 


. 
B 
t 


The resulting model is then 


plt) =E A -« n(i), ШЕ 


where A is the M X l vector of unknown complex modal amplitudes. In Chapter 


2, this model 15 derived rigorously from both range-independent and range- 


dependent normal mode descriptions. 


'The model mmplied by equation (1.5) forms the basis for all the different modal 
beamforming algorithms developed in Chapter 5. In its temporal aspects, this 
model is similar to the one used ш the plane Wave beamforming development. The 


two differ only in their spatial aspects, and then only in the number, and not the 





type, of the unknown signal parameters. 


Even though they have been modeled here as non-random constants, it is 
reasonable to assume that the complex modal amplitudes will actually vary 
randomly over time in any real ocean. The size and structure of these fluctuations, 
particularly the size and structure of the phase fluctuations for the different modes, 
is of fundamental importance to the proper characterization of the acoustic 
propagation from a signal processing viewpoint. Indeed, the measurement and 
analysis of these fluctuations is the ultimate objective of this thesis. The critical 
issue here is Whether or not the phase variations found for any one mode are 
independent of the phase variations found for the other modes. If the fluctuations 
occur independently of one another, then they destroy the phase relationship 
between the various modes; in this case the modes are said to be ¿ncoherent with 
respect to each other. On the other hand, if any fluctuations that occur do so 
simultaneously across all the modes, then the modes remain phase locked with 


respect to each other, and may be considered to be coherent. 


Examination of the coherence of the different modes is important for a 
number of reasons, three of which are mentioned here. First, the issue affects 
propagation modeling. If the modes are incoherent with respect to each other, then 
the individual modal phases can be ignored, and the total field energy should be 
computed by summing individual modal energies. Conversely, if the modes are 
coherent, then both their amplitudes and their relative phases have to be predicted 
accurately, since they are needed to properly compute the total field. Second, the 
issue Is central to the question of the feasibility of direct target range estimation. 
which, as discussed earlier, requires a strongly coherent mode field in order to 
produce accurate results. Finally, from a signal processing point of view, a measure 


of the mode coherence Is a critical element in any complete description of the multi- 
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path environment of the channel. An incoherent mode field indicates that the 
various transmission paths are independent and phase random; the often used 
WSSUS assumption (wide sense stationary, uncorrelated scatterers) implicitly 
assumes this characterization. On the other hand, a consideration of the linear 
relationship between rays and modes leads to the conclusion that the multiple paths 


cannot be considered independent if the mode field remains coherent. 


A good measure of mode coherence can be made from the M Xx M modal 
cross-coherence matrix, which exists in two versions. The unnormalized form of the 


matrix can be defined as 


= _ z | 

5, =Е [АА | (1.6) 
where the + symbol stands for the conjugate transpose operation. The notation 
E[-] nominally represents the average over an ensemble of different trials; 
however, time averages rather than ensemble averages are almost always used when 


implementing the process. Each element of the matrix can then be written as 


S ik =. (А; A, = | A; АН £ = «a (1.7) 


where 

a le: (1.8) 
is the magnitude-phase representation of the complex modal amplitude. The 
elements of the normalized modal cross-coherence matrix are scaled versions of the 


unnormalized matrix elements, namely, 


> Е[А, А] 


Т = ===. (1.9) 
—Mik : Te > 
| ВИА а Еј 





The normalized matrix has two useful properties: its diagonal elements are all unity: 


and the magnitude of any non-diagonal element is not larger than 1.0. 
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To understand how the cross-coherence matrix reflects the phase randomness 

of the various modes, consider a situation where the magnitudes of the modal 
amplitudes are constant and only the modal phases are allowed to be random. 
Then the elements of the normalized modal cross-coherence matrix assume the form 
we az), (1.10) 

If the two random phases are now completely independent of each other, then their 
difference 1s also completely random, so that the expectation takes on a value of 
zero. On the other hand, if the two modes are phase locked, then the phase 
difference takes on a constant value, even though the individual phases themselves 
may be totally random; the resulting magnitude of the matrix term Is one. The 
magnitude of any term of the normalized modal cross-coherence matrix, therefore, 


provides a quantitative measure of the coherence-of the two modes comprising the 


term. 


Two issues remain; that of estimating the modal amplitudes and their cross- 
coherence based on the model of equation (1.5), which is exactly the modal 
beamforming problem; and that of the applying a model that assumes a non-random 
signal to mode coherence estimation, which is a measure of the signal's randomness. 


These subjects are taken up in Chapter 5. 


1.3 Objectives 


The ultimate purpose of this thesis is to measure the coherence between 
different normal modes emanating from a single tonal source and then propagating 
through the central Arctic acoustic sound channel to ranges of several hundred 
kilometers. Such an investigation implies a great deal more than that which is 


explicitly stated. To begm with. if one wishes to measure the properties of the 





normal modes, then one must first understand the modal structure of the sound 
channel of interest; this, in turn, requires thorough investigation of some of tlie 
oceanographic properties of the local environment, including items such as tlie 
channel sound speed profile and surface and bottom descriptions. Obviously, there 
is a need to understand in detail the nature of the available data and the equipment 
and techniques used to obtain it. The received signal energy must be separated into 
the component modes while at the same time rejecting as much background noise as 
possible. This effort turns out to be non-trivial, and much of this thesis is spent 
understanding the methods available, their performance in various situations, and 


their sensitivity to the realities of field research. 


If one has managed to accomplish all of the above, then there is a reasonable 
chance of making some valid measurements of the coherence of the various modes. 
At the same time, greater insight is (hopefully) gained into the nature of sound 
transmission in the channel under consideration, and thus it is valuable to consider 


the acoustic propagation implications of the results. 


This thesis thus has five general objectives: 


1. To describe and understand the modal structure of the central Arctic 
sound channel encountered during the FRAM IV Experiment, including 


the effect of different environmental aspects on this structure and its 
implications for the source and receiver geometries involved: 


to 


. To develop the beamforming methods needed to make direct modal 
amplitude estimates for narrowband tonal signals, particularly when 
vertical arrays are employed, and to answer some of the performance 
questions regarding these techniques; 


3. To assess the operational utility of vertical arrays, and to understand the 
impact of some of the practical realities encountered in their use, 
particularly with respect to the modal decomposition process: 


t. To auswer the question of whether or not the amplitudes of the various 
modes generated from from a single harmonic point source remain 
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coherent after propagation to long ranges in the Arctic sound channel: 
and 


5. To assess the implications of the modal amplitude and coherence 
estimates made in light of the current understanding of central Arctic 
sound channel propagation. 


1.4 Contributions 


This thesis provides a number of significant contributions to the fields of 
underwater acoustics and signal processing. These contributions mav be classified 


into three areas. 


First, there are several contributions made to the discipline of signal 
processing and, in particular, to beamforming theory. The modal beamforming 
algorithms developed in Chapter 5 are not really-new. Hinich [42], [43], Clay [17]. 
and Bucker [9], among others, have all covered similar ground. In any case, the 
results are straightforward analogs of earlier plane wave processing results examined 
by such a large number of other investigators, such as Schweppe [7-4], Capon [10]. 
and Baggeroer [2]. Rather, it is the performance evaluation of these estimates 
which is unique in this work, especially because it is made for a real array and a set 
of modes developed from a real channel. The assessment of mode resolution 
provides some valuable insight into which modes can be distinguished and which 
cannot. The study of the numerical stability limits to the number of modes that 
can be simultaneously included in a multiple beam beamformer contributes in a 
similar way. The performance relationships of multiple beam and single Белш 
algorithms, especially for the MILM approach, is a subject not adequately addressed 
previously. Finally, the examination of the difficulties of the MLM algorithm in the 
face of coherent modal signals provides a different perspective on another important 


signal processing Issue that is all too often incorrectly ignored. 





The second area of contribution deals with the practical issues involved in the 
use of vertical arrays. Most important here ıs the finding that modal decomposition 
is extremely sensitive to array tilt. This result has tremendous operational 
significance, because it indicates that the shapes of vertical arrays must be known 
far better than they are presently, if the devices are to be exploited to anvthing 


near their theoretical potential. 


The area in which this thesis makes the largest number of significant 
contributions, however, is to the discipline of underwater acoustics, particularly to 
the understanding of how low frequency sound propagates to long ranges in the 
central Arctic Ocean. This thesis represents the first time that modal amplitude 
and phase measurements have been made for low frequency tonal signals using 
actual field data from a vertical array and direct modal decomposition techniques. 
Previous efforts involving modal amplitude estimation for tonal sources, such as the 
shallow water experiment conducted by Ferris [33] and the laboratory investigation 
of Hobaek, Tindle, and Muir [41], have all allowed the channel itself to accomplish 
the modal separation by utilizing pulsed sinusoids which then separate temporally 
into the various modal arrivals because of different modal group velocities. The 


present effort obviously has much wider scientific and practical application. 


Beyond the intrinsic value of modal beamforming as a tool for further 
scientific investigation, additional contributions are made in understanding the 
process of low frequency transmission in the central Arctic. In particular, both the 
modal amplitude and coherence estimates made here are the first of their kind to be 
attempted. and provide some very new and very different msights into Arctic Ocean 
propagation. The picture of an essentially determmistic sound. field. that. emerges 
from the coherence measurements and the stability observations not only verifies 


Miklialevsky s earlier findings [57]. but also expands on them considerably, since the 
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various propagation paths have now bcen at least partially separated. The 
classification of two significantly different propagation regimes, one below 40 Hz 
and the other above, has practical implications beyond its usefulness in cataloging 
Arctic propagation; these are discussed in the conclusions of Chapter 8. The ice 
scattering effects displayed in the results provide new insight into a very important 
problem that has not been satisfactorily solved. Of the most interest here are the 
indications that sound energy scattered by the ice remains coherent with the 
specular field, and may play a significant role in the overall sound transmission 
picture. Finally, the in depth study of the modal structure of the central Arctic 
conducted in Chapter 4 provides new perspective on the role that the strong Arctic 


surface duct plays in channel propagation. 


1.5 Thesis Organization 


The remainder of this thesis consists of seven chapters and two appendices. In 


general, each chapter covers one major aspect of the overall analysis. 


Chapter 2 lays the theoretical groundwork for modal beamforming techniques 
by providing a mathematical description of the sound field generated at a distant 
array by a point source that is located in a waveguide. This description. which is 
based on acoustic normal mode propagation theory, is then couched in a forin to 


which optimal estimation theory can easily be applied. 


Chapter 3 provides a full description of the FRAM IV Arctic Acoustic 
Experiment from which the data set of interest is drawn. Most of the important 
technical details of the experiment (navigation data, hardware descriptions, etc.) are 
also provided. A preliminary analysts of some of the practical issues involved in 
dealing with suspended array systems is included. This is done to highhieht the 


issues” relative importance to the modal beamforming problem when using both 
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horizontal and vertical arrays. The preprocessing scheme used to compress the data 
set to a usable size is discussed. Finally, some of the resultant time series are 
analyzed to provide preliminary insight into the modal energy distributions and 


coherence values that might be expected. 


Chapter 4 presents a detailed analysis of the significant environmental aspects 
of the central Arctic sound channel. An in depth study of the modal structure 
encountered during the FRAM IV Experiment is made for the source-receiver 
geometries of interest. This study is a necessary preliminary to any modal 
beamforming effort, since it defines the specific spatial structures of interest 


throughout the rest of the thesis. 


In Chapter 5, several different modal beamforming algorithms are developed 
theoretically from optimal estimation theory. Their theoretical performance is then 
analyzed, leading to the selection of a multiple beam least squares technique as tlie 


most suitable method for processing the data from the vertical array. 


Chapter 6 studies the most important practical issue that arises when 
attempting modal beamforming with vertical arrays; that of array tilt. The modal 
decomposition process is shown to be extremely sensitive to the effective tilt angle. 
Since no direct. tilt ТУОС were made during the FRAM IV experiment, a 


relatively simple method for its estimation from the acoustic data is developed. 


Chapter 7 presents the outputs of the modal beamforming processor when 
applied to the actual data set. Horizontal and vertical array measurements are 
contrasted, and both are compared against. theoretical predictions. Some 
corroborating shot data is also presented. Fhe modal amplitude and mode 
coherence results are analyzed to understand ther implications for central Aretic 


sound propagation. 
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Chapter 8 provides a summary of the results and a general discussion of the 


final conclusions of the thesis. 


Each of the two appendices deals with a particular mathematical aspect of 


modal beamforming theory needed to make the analysis rigorous. 
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Chapter 2 


A Modal Source Description 


The purpose of this chapter is to rigorously develop, in a form that is 
tractable within the context of signal processing theory, a mathematical model to 
describe the sound pressure field that a harmonic point source imbedded in a 
waveguide generates at a distant receiving array. The final result of this 
development has already been presented (in a somewhat simplified form) in 
equations (1.1) and (1.5). The present chapter is designed to provide a more 
detailed presentation of the implied assumptions and the line of reasoning that leads 
io this result. The development consists of two parts. First, normal mode 
propagation theory is reviewed for both the range independent and range dependent 
channels. Normal mode theory is the natural starting point for this thesis, since it is 
modal characteristics that are of interest here. Besides providing a solid link from 
the various forms of normal mode theory to the present effort, the review also 
provides an opportunity to introduce the nomenclature and notation that is used 


throughout the rest of the work. 


The remainder of this chapter is concerned with the conversion of the various 
modal propagation theories to forms that describe the sound pressure field in the 
immediate vicinity of a distant receiving array. This conversion is important, as it 
represents the boundary between theories in the underwater acoustics domain and 
those of the signal processmg domain. The results then become the basis for all 
further work in this thesis. The assumptions and approximations needed to 
complete the conversion are carefully presented. and the change in notation to the 


matrix forms typically used in signal processing is described. 
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In its simplest form, that of a hard or non-propagating bottom, range 
independent normal mode theory is a simple application of tle standard separation 
of variables approach to tne solution of the wave equation. Pekeris, in his classic 
paper [64], first showed how the concept could be extended to the case of a more 
realistic ocean bottom, that of an infinite fluid half space with a given sound speed 
and density. Today, of course, any number of standard references include 
discussions on the subject [79], [18]. Range dependent normal mode theory is 
somewhat more recent and not as well documented. Originally proposed by Pierce 
[65] and Milder [58], it is an area of active interest that has been used to investigate 
the effect of ne variations in. the surface, in the water column. and in the 
bottom. Examples of recent work include Rutherford [72], Dozier and Tappert [27]. 


[28], and Beilis and Таррег [5]. 


2.1 Normal Mode Propagation in a Range Independent Channel 


Consider the channel illustrated in Figure 2-1. Ап idealized oceanic 
waveguide is assumed to be horizontally stratified and symmetric with respect to 
the angle dimension of a cylindrical coordinate system, making only the range and 
depth dimensions of interest. The waveguide is of depth HM with boundaries 
ОПЕ те Of a pressure Telease surface at z — O and a bottom at = H. The 
bottom may be either hard (non-propagating) or soft (propagating), but the former 
Is assumed for ease of development. A discussion of the propagating bottom 
development Is presented in the sequel. An arbitrary sound speed variation that is a 
function of depth, but not of time or range, is assumed. Por simplicity, a density 
Which is constant over depth is assumed; the approach can be easily extended to 


melude densities that vary with depth. 


Let a harmonie point source of frequency ip ОТИ УТЕ Ор Ihe 
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Figure 2-1: An Idealized Range-Independent Oceanic Waveguide 


wave equation for the sound pressure may be written as 


9 [д ar eC T 
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where 
р == p(r,z,t) (БАУ 


and where S is the sound pressure tevel in «Pa at a distance of 1 meter from the 
source in an infinite uniform medium (the units of S are ¿Pa m). The factor of 2 


arises from the different normalizations of the impulse function im the evlindrical 
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and spherical coordinate systems. The appropriate boundary conditions are 


p(r,0,t) = 0, (2.3) 

O, 

E ок (2.1) 
0 


and a radiation condition at infinite range. 


Taking the temporal Fourier transform of equation (2.1) and evaluating it at 


the frequency of interest yields the Helmholtz equation 
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where 


Range independent normal mode theory for the hard bottom case then arises as a 
natural separation of variables solution to equation (2.5). The depth functions о42) 
which provide the various mode shapes are the eigenfunctions of the one- 


dimensional Helmholtz equation 
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with appropriate boundary conditions derived from equations (2.3) and (2.1). The 
horizontal wavenumbers &, of the various modes are related to the associated 
eigenvalues. Because of the nature of the boundary conditions, it Is easy to cast this 
problem into a Sturm-Liouville form [11], so that the mode shapes form a complete 
orthonormal (CON) set. This CON set has three important properties. To begin 


with, the functions within the set are orthogonal, so that 


—— (4: = 0 [for 7 x k. кх) 
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E 
The density gy is included here in analogy to the results obtained when the density 
is allowed to vary with depth; for that condition, the reciprocal of the density 
function acts as the weighting function in the orthogonality relation for the mode 
shapes of pressure. Note that this differs from the results obtained for the mode 
shapes of acoustic velocity potential, where the density term appears in the 


numerator of the orthogonality relation rather than the denominator. 


In addition to being orthogonal, the functions are also normalized, so that 
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Po 





Finally, and most importantly, the completeness or closure property guarantees 
that any arbitrary function of x can be completely represented by an appropriately 


weighted sum (perhaps infinite) of the various ЖЕР 


Decomposing equation (2.5) on the mode set and using the delta function 


expansion 


3 2210) 
=] 0 
leads to the range equation for each mode 
1d (4h) > 2 dr) 
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This is a Bessel equation of order zero, having the two types of Hankel functions as 
its solution. Such a result is expected, since the range dependence represents the 
cylindrical spreading of each mode. The radiation condition requires that the 
Hankel function of the second kind be selected. for strictly outward propagation. 


Note that if the temporal exponential in equation (2.1) were chosen to have the 





ES: 
opposite sign, the proper choice to meet the radiation condition would be the 
Hankel function of the first kind. At distances of more than a few wavelengths 


from the source, the Hankel function may be replaced by its asymptotic equivalent 


(eae) 





vielding 
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where 
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Each term of the result has been conveniently arranged in three parts. br) 
plays the role of a range dependent modal amplitude, while ó(z) provides the 
normalized depth dependence that represents the shape of the mode and tlie 
complex exponential provides the traveling wave phase advance. While the purpose 
of this development has traditionally been to compute a sound pressure level as a 
function of depth, the point of interest in this thesis is rather the internal form of 
each term and, particularly, the result for the modal amplitude given in equation 


(2.14). 


The resulting solution thus provides the following description of the 
propagation (shown in Figure 2-2): 


1. Тре source excites each mode to a level proportional to the size of the 
mode shape at the source depth; and 


2. Each mode then propagates Independently outward from the source as a 
non-homogeneous cylindrical wave (since the wave amplitude varies with 
depth) having its own unique wavenumber and depth dependence. 
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Figure 2-2: Range-Independent Modal Propagation 


Although the expression in equation (2.13) nominally extends over an infinite 
mode set, only a finite number of modes actually propagate. Beyond a certain 
mode number, the associated modal horizontal wavenumbers turn imaginary, 
vielding a decaying exponential in range rather than a propagating wave solution, 
Such modes are known as being in cutoff. At long ranges from the source, the 
contributions to the sound pressure field from these modes in entoff become 
negligible, effectively limiting the sum to a finite (although possibly large) number 
of modes. The number of modes that need to be included in the sum may be 


further reduced by other effects on an individual case basis. Two examples are 





-38- 
selective modai excitation due to the deployment depth of the source and physical 


attenuation processes which favor some modes over others. 


Having completed the development for the condition of a hard bottom, the 
more physically realistic case of a propagating bottom must now be addressed. A 
full development for a propagating bottom along the lines of Pekeris [64] has been 
avoided, as it requires much more mathematical complexity yet generates no 
significant differences in the final result. Instead, the propagating bottom results 
are simply contrasted with those presented above to provide some insight into the 


relationship between the two problems. 


The primary difference between the hard bottom condition and the 
propagating bottom condition occurs in the definition of the CON mode set. The 
lower boundary condition associated with equation (2.7) in the latter case no longer 
allows it to be classified as a Sturm-Liouville problem, so that a CON set of mode 
functions cannot be guaranteed. However, it is still possible to define a set 
consisting of a finite number of trapped mode shapes. These mode shapes can be 
made orthonormal through proper scaling, so that equations (2.3) aud (2.9) are still 
applicable if the range of integration is extended to infinity. This is reasonable, 
since, for an infinite. half space bottom, the resulting sound pressure field also 


extends infinitely in depth. 


The set of trapped modes is not complete, however, so that although the 
results presented in equations (2.13) and (2.14) are «applicable to the portion of 
energy carried in the trapped modes, it is possible for energy not resident in the 
trapped modes to also contribute to the sound pressure field. This additional 
contribution is carried by the continuum of modes that propagate into the bottom, 
aud a further tern: involving an integral over the infinitely dense continuum. mode 


set must be included in equation (2.13) to account for it. The contribution can be 
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expected to attenuate with range, however, as more and more of the continuum 
energy is lost due to propagation into the bottom. At long ranges, the continuum 
may be neglected, so that equations (2.13) and (2.14) still pertain; but with the sum 
limited to the trapped modes. Thus, the net effect of the propagating bottom is 
twofold: to change the range of the summation from the modes not In cutoff to the 


trapped modes; and to increase the minimum range at which the result is valid. 
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Figure 2-3: An Idealized Range-Dependent Oceanic Waveguide 
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2.2 Normal Mode Propagation in a Range Varying Channel 


The most restrictive assumption made in the previous section is the 
characterization of the channel properties, including boundary conditions, as being 
independent of range. To relax this assumption, one can assume a sound speed that 
is now a function of both depth and range (see Figure 2-3). This change makes the 
resulting partial differential equation inseparable, so that the approach of the 
previous section must be abandoned. Instead, a partial separation of variables is 
invoked by allowing the various mode shapes to vary in range as well. In other 
words, as the modes propagate outward from the source, their shapes are now 
allowed to change as they encounter channel range variations. To find the range 
dependent mode shapes, one must solve the equivalent of the one-dimensional 


Helmholtz equation (2.7), now parameterized in range 
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Suitable boundary conditions must also be included. The approach requires that 
the associated eigenvalues, and thus the modal horizontal wavenumbers, also be 
functions of range. If the discussion is again restricted to a hard bottom boundary 
condition, a CON set of mode shapes can be found at any particular range. 
Because the separation of variables is only partial, the resulting range equations 


remain coupled; following Rutherford [72], one obtains 
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The assumptions implicit in the derivation of these coupling parameters are a 
constant density and range invariant boundary conditions. Allowing the boundary 
conditions to vary with range affects the form of the coupling coefficients [72]. 
However, most problems mvolving range varying boundary conditions can be recast 
to have range independent boundaries by the proper redefinition of sound speed and 
density. The comments made previously concerning the inclusion of soft or 


propagating bottoms are, mn general, also applicable here. 


Because the range equations presented in (2.17) are coupled. they are not 
typically amenable to analytical solution, and further development beyond this 
point involves erther numerical evaluation or the application of approximations. As 


might be expected, the resulting general solution Is not very revealing: 





Be!) lo EOS (2.20) 


One particular approximation generates an important subset of coupled mode 
theory that is more analytically tractable. Adiabatic mode theory arises by ignoring 
the effects of the coupling coefficients in (2.17). This is equivalent to assuming only 
very gradual channel variations in range, which is often the case in practice, 
Dropping the coupling terms once again leads to independent range equations. ln 
addition. the assumption of gradual range variations requires that, in the large 


range limit, the solutions of these equations must be asymptotic to 
fg ke cd 
~ o es ә ә 
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The full adiabatic solution is then 
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Unlike the range independent case, it is not possible to find a general closed form 
solution for 6{r). The nature of the coefficient preceding the range delta function 
on the right hand side of equation (2.17) does guarantee that b(r) are proportional 


to both the source level S and the size of the mode shape at r= 0. 


Adiabatic mode theory thus generalizes the physical picture of mode 
propagation as follows (as illustrated in Figure 2-4): 
l. The source excites each mode to a level proportional to the size of the 
mode shape at the source depth, where the mode set 1s evaluated for the 
source sound speed profile: and 


2. Each mode then propagates independently outward from the source as a 
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Figure 2-4: Range-Dependent Adiabatic Modal Propagation 


non-homogeneous cylindrical wave (since the wave amplitude varies with 
depth) having its own unique wavenumber and depth dependence. Both 
of these are now allowed to vary with range. It is important to realize, 
however, that in this approximation the modes still propagate 
independently. 


Full coupled mode theory further complicates the picture by including mode 


conversion effects. In this case, the various modes are allowed to exchange energy 


in addition to varying their shapes and wavenumbers as they propagate down the 


channel, so that they no longer propagate outward In a completely independent 


fashion. 
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Figure 2-5: Relationship of Source and Receiver Coordinate Systems 
(top view) 


2.3 Modal Source Description at a Far Field Receiving Array 


While the descriptions developed above are quite uscful when dealing with 
underwater acoustic propagation theory, a further step Is necessary to convert them 
into a form useful in the context of signal processing. Array theory deals with the 
field at the receiver rather than at the source; therefore. it is appropriate to shift 


coordinates to a system with its origin at some arbitrarily defined. reference. point 





12 
near the receiving array, and then to expand the field representation in the vicinity 
of this new origin. Several approximations and assumptions can then be used to 


simplify the resulting expressions. 


Consider the diagram of the pair of two-dimensional horizontal coordinate 
systems shown in Figure 2-5. The (r,6) coordinate system at the left represents the 
polar coordinates in which the propagation problem has been solved. The (x.y) 
coordinate system on the right represents the new receiver coordinates, with an 
origin at the arbitrarily defined receiver reference point. This point is assumed to 
be located at range ry and angle 9, in source coordinates. "The angle 7 is the 
horizontal direction of propagation at the receiver reference point with respect to 
the receiver coordinate svstem. 

From the law of cosines, the source range г for any point (x,y) in receiver 


coordinates may be found to be 





1 + 2(2/r,) cos э + Ay/ro) sin 9 + (2/1, + (y/n. (2.23) 


In the immediate vicinity of the receiver reference point, the receiver coordinate 
offsets r and y are small when compared to rj For this condition, the terms 
containing these offsets are of first order, while the terms involving their squares are 
of second order. Expanding the square root as a binomial series and dropping all 


terms beyond first order leads to 


reasr, + 2 cos p + ysm om (9 2s) 


Neglecting the higher order terms is equivalent to approximating the cylindrically 

о о о “ ` 
propagating modes with ones that are planar in the immediate vicinity of the array. 
In the sequel, the approximation in. equation. (2.21) is applied to all the sensors in 


the receiving array. To guarantee that the approximation is valid. two restrictions 





Aen 
are needed. First, the receiver reference point must be chosen to be close to the 
аггау; second, the array aperture must be small compared to the source-receiver 


range. Both of these conditions are easily met in practice. 


Conversion of the range independent modal solution to the new coordinate 
system simply requires the replacement of r in equations (2.13) and (2.14) with the 
expression given in (2.24); retaining first order accuracy in the phase but limiting 


amplitude terms to order zero only, one obtains 
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is of magnitude kK. and points in the direction of propagation. The modal amplitude 


is now a function of only rp 
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For the range dependent cases, an additional assumption of the horizontal 
homogeneity of the channel over the aperture of the array is needed for tractability. 
Given this assumption, the adiabatie theory solution (equation (2.22)) can then be 


converted to 
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Here the magnitude of the wavenumber vector is the value of the horizontal 
wavenumber evaluated at the receiver reference point k (rg). Conversion of the full 
coupled mode solution (equation (2.20)) is more complicated, but results in the same 
form as for the adiabatic solution with a different definition of the modal amplitude 
a; Because the channel is assumed to have no range variations over the aperture of 
the array, the two sets of coupling coefficients A., and В. vanish in this region. А 
far field representation of the solution of equation (2.17) that is valid over the 
aperture of the array can be written as 

АУТО ЫЯ (2.31) 

1 1\ 0 , > 

where any energy scattered back to the array from channel range variations that 
are beyond it has been ignored. The resulting solution is then has the same form as 
equation (2.29), with the modal amplitude now defined as 


m lo). (2.32) 


Again, the magnitude of the wavenumber vector is that of У 


It can be seen that all three approaches lead to the same representation in the 
region of the receiver. In vector notation this may be written as 


Дт) (Е) Ае (2.33) 


И.М modes are included in the model, A becomes an. M X 1 column vector, each 


element of which is one of the complex modal amplitudes a. 9(r.z) i» also an 


Mx 1 columna vector with elements 
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„ (7,2) = 6 (2) e , (2.34) 


where all mode shapes and wavenumbers are those for the vertical channel structure 
found at the receiver. While the functional form of the modal amplitudes d; varies 
depending on the particular theory chosen, in all cases they may be considered to be 


constant across the aperture of the array. 


The bridge from equations (2.33) and (2.34) to equations (1.4) and (1.5) is now 
straightforward. The total signal received on any hydrophone is assumed to consist 
of the signal model just developed plus additive noise. If equation (2.33) is then 
evaluated at each of the N hydrophone locations in the receiving array, the results 
collected into an N X 1 complex vector, and a quadrature demodulator used to 
remove the harmonic time dependence, equation (1.4) 1s recovered exactly. Each 
row of the N X M steering matrix is comprised of the v* vector evaluated at a 


different sensor location 
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Comparison with equation (1.5) allows 8., to be identified as 
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2.4 Summary 


A mathematical model to describe the sound pressure field that a harmonic 
point source imbedded in a waveguide generates at a distant receiving array has 
been developed from first principles. This model has been presented in a form to 
which optimal estimation theory can be apphed. It has been shown that range 
independent mode theory, adiabatic mode theory, and coupled mode theory all lead 
to the same model; only the nature of the modal amplitudes is different (equations 
(2.28), (2.30), and (2.32), respectively). To achieve this result, both of the range 
dependent theories require an additional assumption of horizontal homogeneity of 
the channel characteristics over the aperture of the array. Subject to this 
assumption, the modal amplitudes can be considered constant across the array. For 
the range independent and adiabatic results, the modal amplitudes have been shown 
to be directly proportional to both the source level and the size of the mode shape 


at the source depth. 


The long range or far field assumption has been used to justify a number of 
approximations. These include the neglect of either the modes in cutoff or the 
continuum modes, as applicable; the use of asymptoticly equivalent forms for the 
Various range functions; and the acceptance of a local plane wave approximation in 


place of the actual cylindrically spreading mode. 


As seen by the array, each mode appears to be a non-homogeneous plane wave 
propagating horizontally away from the source with unique horizontal wavenumber 
and depth dependence. The depth dependence of each mode ts defined by its mode 
shape function; the various mode shapes have been shown to be orthonormal at airy 


particular range. 





Chapter 3 


Ihe FRAM IV Data Set 


The purpose of this chapter is to describe the details of the experiment during 
which the data set to be analyzed was taken. The nature and general layout of the 
experiment are first discussed. Following that are descriptions of the acoustic 
source and the two receiving arrays available, one a horizontally deployed two- 
dimensional array and the other a vertical line array. Soine important conclusions 
about practical differences between horizontal and vertical arrays are included. The 
details of the data acquisition and storage system and a description of the actual 
signals analyzed in this thesis are then presented. A discussion of the preprocessing 
employed prior to beamforming and an initial analysis of some of the preprocessed 


data complete the chapter. 


3.1 The FRAM IV Experiment 


ЕКАМ IV, conducted in the spring of 1982, was one of a continuing series of 
multi-institutional Arctic Ocean research projects sponsored by the Office of Naval 
Research. The expertment was designed to study the low frequency acoustical and 
geophysical properties of the central Arctic environment. Additionally., several 
physical oceanography and bottom geology experiments were conducted: besides 
characterizing the central Arctic region in their own right, these secondary 
investigations were designed to provide concurrent measurement of the 
environmental parameters necessary to characterize the acoustic channel from а 


theoretical standpoint. 


The experiment was conducted from two ice camps located on the central 
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Figure 3-1: General Locations of the FRAM and TRISTEN Ice Camps 
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Arctic permanent pack ice. The larger of these camps was named FRAM., and was 
established in the Barents Abyssal Plain, as shown in Figure 3-1. The hydrophone 
receiver arrays and the data acquisition systems were located here. A second canıp. 
known as TRISTEN, was established about 300 km to the west of FRAM near the 
Mid-Arctic Ridge, primarily for the purpose of deploying a high power. low 
frequency acoustic source used to transmit tones and other, more complex 
waveforms. Explosive charges were also employed to generate impulsive source 
data for both acoustic and geophysical analyses. These charges were set off from 


Various sites. 


The FRAM IV data set has formed the basis for a number of recent papers: 
several cover topics that are germane to this thesis. The report by Tiemann, Ardal, 
Allen, and Manley [78] provides a full analysis of the navigation data for the 
experiment. Dyer [32] and Makris and Dyer [50] summarize what is currently 
understood about the nature and causes of Arctic ambient noise. Duckworth 
[29] and Duckworth and Baggeroer [30] provide detailed investigations of the 
bottom characteristics in the region near the FRAM camp. КОК: and DiNapoli 
[55] and DiNapoli and Mellen [25] characterize the propagation loss for the central 
Arctic and attempt to quantify the effect of the Arctic ice canopy on this 
propagation loss. Mikhalevsky [56], [57] investigates the temporal stability of tonal 
signals propagated in the Arctic channel in some detail; although the data exaniined 
in these studies were taken from the earlier FRAM II experiment, tlie results are 
equally applicable — here. Polcari [67] characterizes the structure of the 
TRISTEN/FRAM acoustic transmission path. Finally, Yang and Giellis [St] have 
recently attempted modal decomposition from the FRAN IN vertical array using 


spatial Fourier techniques, with mixed success. 
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Figure 3-2: Drift of the FRAM and TRISTEN Ice Camps 


3.2 FRAM IV Navigation Data 


The FRAM receiving camp was established at 84.1 N, 23.1 E on 23 March, 
1932 (see Figure 3-2). It drifted approximately 260 km to the southwest over the 
course of the experiment, and was abandoned at 82.6 N, 6.9 Eon 11 May. The 
URTEN source camp was-established at 84.29"N, 4.6 W on 26 March; it was 
manned until 26 April, during which time it drifted 90 km south to 823.1. N, 5.8 М. 
The acoustic data of interest was transmitted from TRISTEN to PRAM during two 
separate periods of time, 6 to & April and 16 to 19 Apr. Figure 3-2 indicates 
representative positions for each of these periods, while Tables 3-1 and 3-H tabulate 
the full navigation data for the same positions. All the data presented here has 


been derived from reference [75]. 
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Date: 6 Арг 82 
Time: 1800 
ЕКАМ Рози: 83.807 М 16.82 ° Е 
TRISTEN Posit: 83.84°N 5.55 °\ 
Array Rotation: 15.6° east of true North 
TRISTEN Azimuth Angle: 86.5 ° 
TRISTEN/FRAM Range: 265 km 


Date: 7 Apr 82 
Time: 2100 
FRAM Posit: 83.76°N 16.08°E 
ЛЕТ ВОСЕМ 5.88" \W 
Array Rotation: 15.9° east of true North 
TRISTEN Azimuth Angle: 86.1” 
TRISTEN/FRAM Range: 262 km 


Date: & Apr &2 
Time: 1530 
FRAM Posit: 83.72°N 15.06 Е 
ЛЕ о Е мым 
Array Rotation: 16.2° east of true North 
TRISTEN Azimuth Angle: 86.17 
TRISTEN/FRAM Range: 262 km 


Table 3-I: TRISTEN/FRAM Navigation Data for the Horizontal Array Data Set 


3.3 Horizontal Azimuth Angle Conventions 


Before continuing, a short comment on the measurement conventions used in 
this work for horizontal direction angles is in order. Whenever Incoming signal 
horizontal azimuth angle values are quoted, the convention implied throughout the 
thesis is that the azimuth angle has been measured from the reference direction to 
the direction in which the signal is propagating ina clockwise direction. This 
makes the azimuth angle exactly the reciprocal angle of the more traditional 
bearing angle. As a simple example, a signal propagating from a source to the 


northwest arrives on a bearing of 315°; since the signal itself is traveling towards 
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Паје: 16 Арг #2 
Типе: 1430 
ЕКАМ Рози: 83.555" М 15.04°E 
TRISTEN Posit: 33.66 570 лу 
Array Rotation: 17.6” east of true North 
TRISTEN Azimuth Angle: 85.6 * 
TRISTEN/FRAM Range: 255 km 
Date: 17 Apr 82 
Time: 1930 
FRAM Posit: 83.49°N 1486 Е 
EDIEN Poit: 5960 N 5.69 W 
Array Rotation: 17.4° east of true North 
TRISTEN Azimuth Angle: 86.8 ° 
TRISTEN/FRAM Range: 254 km 
Date: 19 Apr &2 
Time: 0300 
ЕКАМ Posit: 83.48°N 14.77 ° Е 
ПРВИ А Роза 8361 N 5.59 W 
Array Rotation: 17.1” east of true North 
TRISTEN Azimuth Angle:. 87.2 ° 


TRISTEN/FRAM Range: 253 km 


Table 3-11: TRISTEN/FRAM Navigation Data for the Vertical Array Data Set 


о 


the southeast, however, its azimuth angle is 135 


Azimuth angles are reported in a range from 0” to 360 , and are referenced 


to the north leg of the horizontal array rather than true North. The reference ıs 
chosen so as to be fixed with respect to the array. Because the ice from which the 
array is suspended is free to rotate relative to the global coordinate system. the 
relationship between the array and true North is not static. In actuality, the net 
array rotation turns out to be small. The direction of the array’s north leg varies in 
a band of only 2.5 around Its median value of 16.4 east of true North during the 
period of the main acoustic experimentation: somewhat larger shifts are observed 
None the less, even these small variations 


towards the end of the experiment. are 





26 
significant, especially when high resolution beamforming techniques are employed. 


Tables 3-I and 3-JI both include appropriate array rotation data. 
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Figure 3-3: Source Levels for the TRISTEN Source 


3.4 The TRISTEN Low Frequency Source 


A modified TILF low frequency source was deploved to a depth of 91 m at the 
TRISTEN camp. Approximately 150 hours of signals were broadcast in a frequency 
range from 5 to 110 Hz. Several different waveform types are included in the signal 
set. Of interest here are a series of tonal signals m the range from 5 to 71 Hz. The 
typical signal format consists of one tone per hour. with the source on. continuously 
for the first 55 minutes and then off for the last 5 minutes of the hour... Scheduling 


constraints often. forced more abbreviated formats: most of the horizontal array 





N: 
data used here was broadcast at half-hour intervals rather in the full hour format. 
for example. The source has a rated frequency stability of better than one part in 
10°: it can therefore be considered to be temporally coherent over the full duration 
of any signal investigated in this thesis. Figure 3-3 provides a plot of source levels 
versus frequency in the range from 5 to 80 Hz for the TRISTEN source. This data 


was extracted from reference [26]. 
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Figure 3-4: Arrangement of the FRAM Receiving Arravs 
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3.5 The FRAM IV Receiving Arrays 


Two receiving arrays were available at the FRAM camp. A two-dimensional 
horizontal array consisting of 26 elements was deployed through the tce and 
suspended on cables to a depth of 91 m. The net gain to the recording system was 
-159 dB ге І v/,Pa, with a flat frequency response to below 2 Hz. As shown in 
Figure 3-4, the array was deployed in an cross-shaped configuration. Each leg of 
the array had an aperture of about | km, with individual elements spaced 


approximately logarithmically from the center. 


The second array was made up of 28 elements and extended vertically into the 
water column to a depth of approximately 1000 m. Since the bottom depth at the 
receiver camp averaged 3800 m, this г. represents the top 25° of the 
channel. The elements in the array were spaced approximately linearly.” The array 
was located about 30 m to the southwest of the center of the horizontal array, as is 
also indicated on Figure 3-4. The net gain from the water to the recording system 
was -121 dB re 1 v/„Pa, with the exception of a few particular hydrophones where 
the gain was reduced to -127 dB in order to reduce amplifier saturation. The 
frequency response was flat above 6 Hz. A 500 lb weight was hung from the bottom 


of the array to stabilize it agamst current action. 


Both arrays suffered to varving degrees from two related problems endemic to 
suspended sensor systems, sensor displacement and strum. Sensor displacement is 
the offset of a particular livdrophone from its nominal position caused by current 
drag on the hydrophone and its support cable. Since only the position at which the 
sensor is deployed though the ice is known, such displacements translate into phase 
errors In the beamforming process. Simple delay-and-sum) beamforming theory 


mdicates that sensor displacements become significant when they begin to reach 
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magnitudes of approximately one quarter wavelength; the problem obviously gets 
more severe with increasing frequency or more sophisticated beamforming 
algorithms. Modal decompositions may be even more affected than other techniques 
when, as for the vertical array, the errors occur across the array in a consistent 
fashion, causing the array to become tilted. This sensitivity is investigated in detail 
in Chapter 6. Time varying sensor displacements are possible due both to 
variations in the environmental conditions and possible pendulum oscillations of the 
sensors. <As discussed below and in Chapter 7, however, the time variations of 
sensor displacements observed in the data appear to be relatively small when 


compared to the static offsets. 


Strum is the vibration of the suspension cable of a sensor due to vortices shed 
bv the surrounding current flow. This vibration causes slight oscillations in. the 
sensor depth; since the ambient pressure is a function of depth. these depth 
excursions cause pressure field variations to which the sensor then responds. As 
opposed to sensor displacement, which occurs on time scales so large that its 
primary effeet is a phase modulation of the received. signal, strum energy is 
concentrated at the lowest frequencies of the acoustic domain, and forms an 
additive noise. source. from. which. the. true. signal. must. then. be separated. 
Strumming action generally appears to consist of a fundamental frequency and a 
number of strong harmonics. These harmonies typically occupy the frequency range 
below 10 Hz for horizontal array sensors: significant harmonics can be found at 
frequencies as high as 30 Hz during severe episodes, however. Because different 
horizontal hydrophones exhibit varying strum patterns, spatial processing 
discriminates well against strum energy. The greater effective diameter of its 
support cable causes strum effects for the vertical array to oceips a lower frequency 


range: strmum energy above 10 Hz is rarely observed tn vertical array. data. making it 
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much less of a problem than for the horizontal array. These observations agree 
with strum frequency predictions based on the Strouhal number, which vield 
expected strum frequencies of 28.5 Hz for the horizontal array and 5.1 Hz for the 
vertical array for an upper limit current of 0.8 kt [61]. Cable diameters of 0.125 in 
for the horizontal array elements and 0.7 in for the vertical array are assumed in 


these calculations. 


3.6 The Effects of Hydrophone Sensitivity Mismatches 


An important difference between the horizontal and vertical arravs involves 
the effect of variations in hydrophone sensitivity and amplification gain across the 
arrays. For a horizontal array of reasonable aperture, the signal received at each 
sensor is a time shifted replica of that found at other hydrophones. as is 
demonstrated in Figure 1-1. This assumption is often used in horizontal array 
beamforming to eliminate variations in the receiver hardware by artificial 
normalızation of each time series before beamforming. This can be done because a 


horizontal array uses the signal phase rather than its amplitude to develop spatial 


information. 


Since a vertical array samples the modal structure of the channel. however, 
signal amplitude variations between different hvdrophones are both expected and 
desired; a normalization scheme as might be used on a horizontal array is no longer 
applicable. Thus, the nature of the data and the processing make hardware 
uniformity between channels a much more important issue for Vertical arrays than 
for horizontal arrays. No indication of this type of problem is seen anywhere in the 
FRAM IV data set: none the less, the issue must be carefully considered whenever 


one Is dealing with data taken from vertical arrays. 





Ale 
3.7 Sensor Displacement Predictions for the Horizontal and Vertical 


Arrays 


In general, the effect of sensor displacement on the horizontal array is 
significantly less severe than its effect on the vertical array. This ts to be expected, 
because the vertical array extends approximately ten times deeper into the water 
column than the horizontal array sensors; the longer scope then makes it more 
sensitive to current effects. The purpose of this section ts to investigate the 
significance of the sensor displacements likely to be encountered in practice for both 
horizontal and vertical array sensors. As the purpose is not to make precise 
quantitative predictions, but rather to provide some feel for when sensor 
displacements must be considered and when they might be safelv ignored, some very 


:simple hvdrodynamie approximations suffice. 


The actual sensor displacement is the primary parameter of interest for a 
horizontal array sensor. For the vertical array, however, things are somewhat more 
complicated, as many sensors are suspended from the same cable. Obviously. the 
maximum sensor displacement achieved is important: this displacement will 
generally be found at the deepest sensor in the array, although more complicated 
geometries are possible Additionally, however, the distribution of displacement 
Done the vertical array is significant. One simple way of parameterizing из 
distribution is by some type of tilt angle, representing the slope of some notional 
line with which one has replaced the actual array shape. Several different tilt angle 
definitions are possible. The one used here is simply the angle made with the 
vertical by a line drawn through the array anchor point at the surface and. the 
bottom sensor. A more rigorous definition is introduced tm Chapter 6, where the 
justification for such an approximation is discussed. A similar parameter of interest 


for both horizontal and vertical arrays is the slope of the suspension cable at the 
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surface, since it provides a check against field experience. 


The sensor displacement problem may be modeled in a very simple fashion by 
assuming the sensor suspension cable to be a uniform line of length H with one end 
fixed to the ice at the origin and the other end attached to a concentrated mass of 
weight WW that represents the weight normally deployed at the end of the cable for 
current stabilization. This model is applicable for both the horizontal and vertical 
arrays, with one difference: for the vertical array, the mass of the line must be 
considered, whereas for the horizontal array it is negligible. Let w be this weight 
per unit length of the cable. To keep the hydrodynamics simple, the line is assumed 
to have a uniform diameter d, and the current V is assumed to be constant in depth. 
The first assumption is not strictly correct for the vertical array, since the actual 
arrangement of the support cable and the electrical connections from the 
hydrophones to the surface is rather complicated; this is compensated for by using 
some type of effective diameter over the length of the cable. Similarly, basic 
oceanography indicates that the second assumption is not very realistic either. 
However, both of these assumptions will more often cause the displacement 
estimates at the deepest sensors to be too large rather than too small. A secondary 
issue Involves the three-dimensional nature of the array shape, which 15 not reflected 
by the constant current assumption. This problem is not serious, however. because 
in the sequel it is not the actual array shape that ts important, but only its 


projection in the vertical plane of propagation. 


Assuming a steady state situation, the primary effect of the current Is to 
provide a constant horizontal drag force per umt length f on the line and a total 
drag force F on the concentrated mass suspended from it. The drag force per unit 


length of line is estimated from two-dimensional Reynolds number theory [61] to be 
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= Cr рута 


9 ' (3.1) 


where p is the density of water and Cp is a two-dimensional drag coefficient that is 


a function of the Reynolds number 


== =>. 32) 


Here. v is the kinematic viscosity of water. A similar but somewhat less reliable 


estimate may be made for F. 


Given estimates for the two drag forces, one can then solve the resulting 
statics problem for the shape of the line as a function of depth in a straightforward 


manner. The solution 15 


r(z) = alz — b In (1 — z/c)}, (3.3) 
where 
й == Í A (3.1) 
10 
F y - 
ђ == F — ПР И | ) )) 
and 
coe i (3.6) 
w 


For the vertical array, /7 is assumed to be 960 in, and w and Y may be taken to be 
17.8 N/m (4 ут) апа 2224 N (500 Ib), respectively. Using an effective cable 
diameter of 9.7 In and a current of 0.8 kt (a reasonable upper lint for the speed of 
ice movement observed in the Arctic), the drag forces f and F are estimated to be 
Fa N/mm and 10.5 XN. respectively, "These numbers give values of 0.087 Tor a, -LIX m 


lor bo. and 1085 m for c. 
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Figure 3-5: Two Models for the FRAM Vertical Array Shape 


(second model taken from (85]) 


[Note different depth and offset scales.) 
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Equation (3.3) is plotted as a function of depth in Figure 3-5 for the values 
just discussed. The results of a more sophisticated model developed at the Naval 
Research Laboratory (NRL) [85] are also shown for comparison. The two models 
agree quite well in the first 600 m, and then diverge slowly after that. A maximum 
difference of about 10 m between the two predictions occurs at the bottom of the 
array. Several different factors probably contribute to the difference between the 
two models, but the most important is the assumption of a current that is constant 


with depth made in the simpler model. 


As shown in Figure 3-5, the sensor displacements indicated for the lowest 
ћудгорћопез are in the 50 to 60 m range, making them a significant fraction of a 
wavelength at even the lowest frequencies of interest. The array tilt as measured 
from the anchor point at the surface to the bottom sensor is 3.6 '; the slope of the 
array at the surface is 5.0^, which agrees well with field experience, in which slopes 
of no more than about 5 or 6° have been observed. The equivalent values for the 
NRL model are 3.0° and 3.8 ^, respectively. These results are in general agreement 
with other field tests carried out by NRL in the Arctic for the same array [69]. It is 
important to note that the magnitude of the tilt angle is fundamentally proportional 


to the square of the current Y. 


The same approach outlined above can be used for the horizontal array 
sensors. In this case the weight of the cable wis negligible, and the resultant shape 


of the line takes on a somewhat simpler parabolic form 


A) = GIE + 70и - 2р8), an 


Per a horizontal array hydrophone, #7 is 91 m and W max be taken to be 1323.1 N 
fom) tora cable dnimeter of 0.125 tn and a current of 0.8 kt, the drag forces f 


and E are computed to be 0.277 N/m and 1.688 N, respectively. These numbers 
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vield an offset of about 10 m, which corresponds to a quarter wavelength at roughly 
35 Hz, and a cable slope at the surface of 11.4°. These results agree with with 
operational experience, which indicates that beamforming degradations are first 


noticable in the 30 to 40 Hz frequency band. 


Several conclusions can be drawn from this section. The most important is 
that sensor displacements for the horizontal array can be expected to be negligible 
below about 30 Hz, while those of the vertical array can never be totally ignored. 
Additionally, typical offsets for the vertical array can be expected to be about five 
times those of the horizontal array, primarily due to the much longer array scope. 
Finally. the angle that the vertical array makes with the true vertical can be 
expected to be limited to a range of about 5° based on the range of currents likely 
to be encountered in the central Arctic. All of these observations are of importance 


in Chapter 6, where the problem of vertical array tilt is examined in detail. 


3.8 The FRAM IV Data Acquisition System 


The recording system is illustrated in Figure 3-6. It consisted of a 
minicomputer-based digital data acquisition and storage system. Twenty-four 
channels of information, each sampled at 250 Hz, were digitized simultaneously; the 
resulting digital data were then stored on specially formatted 1600 BPI] magnetic 
tapes. An analog filter with a roll off of 48 dB/octave above the corner frequency 
of 80 Ilz prevented aliasing, and a floating point conversion scheme preserved the 
required dynamic range. Phe well matched phase responses of the various channels 
preserved the signal phase synchronization required for beamforming. The input to 
each of the 24 channels was selectable via a patch panel. Several different sensors 
sets were recorded during the course of the expenment. Further technical details of 


the recording system can be found in Prada. et al. [71]. 
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Figure 3-6: The FRAM IV Digital Recording System 


The FRAM IV acoustic data is classified according to the array configuration 
(Le., the sensor set) being recorded when the data was acquired. Two major types 
of array configurations exist. The data taken from 6 to & April were recorded from 
only horizontal array hydrophones. Over this period minor variations were made in 
the exact hydrophone set recorded, resulting їп slightly different array 
configurations for different signals. These variations are not significant when 
considered in. terms of theoretical array. performance, however. The data taken 
from 16 to 19 April was recorded from a mixture of horizontal and. vertical array 


lıydrophones. The details of the sensor set chosen for this configuration are 





Br 

| provided in Table 3-IH. Of the most interest here are the last eighteen channels, 
which were used to record selected hydrophones of the vertical array; the depths of 
the sensors recorded are included in Table 3-III. In processing the data taken from 
this vertical array configuration, only these channels have been used, so that an 
eighteen element fully vertical array has been effectively synthesized. This has been 


done to allow assessment of the operational utility of a true vertical array. 


Channel Description 
1 Horizontal array center phone 
2 Horizontał array northernmost phone 
3 Horizontał array southernmost phone 


4 Horizontal array easternmost phone 
5 Horizontal array westernmost phone 
6 Сеорћопе 

i Vertical array phone at 960 m 
8 Vertical array phone at 860 m 
9 Vertical array phone at 782 m 
10 Vertical array phone at 690 m 
11 Vertical array phone at 630 m 
12 Vertical array phone at 570 m 
13 Vertical array phone at 510 m 
14 Vertical array phone at 450 m 
15 Vertical arrav phone at 390 m 
16 Vertical array phone at 350 m 
17 Vertical array phone at 330 m 
18 N ertieal arrax рһопе аё 270 т 
19° Vertical array phone at 210 m 
20 Vertical array phone at 180 m 
21 Wertieal anıav phone ät 1-10 m 
2р Vertical array phone at 90m 
25 Vertical array phone at 60 m 
24 Vertical array phone at 30 m 


ж . . 
These hydrophones had net gain to recording system of 
-]27 415 ratlier than the nominal -121 dB (re I „Pa/v). 


Table 3-IH: Тһе PRAM IN Vertical Array Recording Configuration 
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3.9 Details of the Vertical Array Data Set 


The objectives of this thesis make the data recorded from the vertical array of 
primary interest. The available data set consists of tonals at fourteen different 
frequencies between 5 and 71 Hz recorded from 16 to 19 April. These are listed in 
Table 3-IV. Many of the frequencies were broadcast several times over the course 
of the four days; only the longest representative data segment at each frequeney has 
been chosen for analysis here. The three lowest frequencies were never detected in 
either the vertical array data set or the more extensive horizontal array data set. 
This is most likely due to insufficient signal to noise ratios at these frequencies. Not 
only does the source exhibit significantly lower signal levels below 15 Hz (as shown 
in Figure 3-3), but, in addition, Arctic ambient noise tends to increase with 


decreasing frequency in this band [32]. 


Frequency Time Date Duration 

3.00 Hz Not Detected 

10.00 Hz Not Detected 

ПЕЛ УН? Not Detected 

15.00 Hz 0310 19 Арг 82 10 min 
1 1430 16 Apr &2 55 min 
20.00 Hz 0225 19 Арг 82 35 mm 
23.50 Hz 0330 18 Apr 82 23 MIN 
30.00 Hz 0320 19 Apr 82 10 min 
© лш; 0130 18 Apr 8? 52 min 
35.50 Hz 1030 17 Арг 82 55 min 
47.00 117 1830 17 Арг 82 Пол 
93.29 Hz 0930 17 Арг 82 55 min 
55.00 Hz 1130 17 Арг 82 52 min 
71.00 [iz 1530 17 Apr 8? 55 min 


Table 3-IV: The FRAM IV Vertical Array Data Set 
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3.10 Details cf the Horizontal Array Data Set 


To confirm the results generated from the vertical array data set, a smaller 
number of tonal signals selected from the more extensive horizontal array data set 
have also been analyzed. The particular segments chosen are listed in Table 3-V. 
In general, the frequencies taken from the horizontal array have been chosen to 
closely match those available in the vertical array data set. Several of the segments 
are rather short, making them inappropriate for use in investigating long term 


signal stability issues; none the less, they are still valuable for comparing received 


signal levels with those of the vertical array. 


Frequency Time Date Duration 
15.00 Hz 2330 7 Apr &2 30 min 

7.79 Hz 1800 7 Apr &2 30 min 
23.50 Hz 1830 7 Арг 82 30 шт 
27.00 Hz 1530 8 Apr 82 30 min 
35.25 Hz 1930 7 Apr &2 30 шт 
47.00 Hz 2030 7 Apr 82 13 min 
вы? 2000 T Apr 32 s mm 
70.00 Hz 1800 6 Apr 82 30 min 


Table 3-V: The FRAM IV Horizontal Array Data Set 


3.11 Preprocessing of the Data 


The extensive nature of the data set (the signal set exceeds 1500 Mbytes of 
data as recorded in the field) makes analysis of the data in a raw form impractical. 
The preprocessing scheme shown in Figure 3-7 has therefore been developed to filter 
and compress the data set to a tractable size. The field data for each sensor 
(sampled MA Hz) IS ШО demodulated and = low pass filtered to produce 4 


complex time series of the signal amplitude and phase. For reasons of efficiency and 


phase stability, the demodulation is implemented mn the frequency domain. through 
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Figure 3-7: The Preprocessing Scheme 


the use of bin shifts in the FFT of the data. This approach has the disadvantage 
that nominal signal frequeney can not be achieved exactly: the 8192 point length of 
the FFT used provides resolution only to within +15 mHz of the destred frequency. 
The need to generate quadrature sinusoids with adequate phase stability over the 
signal durations involved is avoided, however; a typical signal of one hour in length 
spans about 10° cycles. With the frequency domain implementation. this. phase 
stability is controlled only by the frequency stability of the highly accurate sample 
clock. The scheme also increases computational efficiency by allowing. the 
simultaneous implementation of a. Parks-MeClellan. FIR. low pass filter [63] in the 


frequency domain. The amplitude and phase responses of this 3909 point. Filter are 
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Figure 3-8: Parks-McClellan FIR Filter Frequency Response 


presented in Figure 3-8, and its characteristics are provided in. Table 3-VI. The 
filter is designed to provide as narrow a passband as possible while retaining 
adequate cutoff levels in the rejection band. The filter length is limited by the 
particular Parks-McClellan algorıthm available, which сап handle no larger than 
400 point filters. After complex demodulation and low pass filtering, the data is 
then decimated by a factor of 40 (to a 6.25 Ilz sample rate); the FIR filter 


characteristics guarantee that the Nyquist criterion is met. 


Following this first stage of preprocessing, an identical second stage of 


complex demodulation, low pass filtering, and decimation is accomplished. The 





Im: 
second stage processing provides further refinement of the demodulation frequency 
and a tighter low pass filter, allowing reductions in the number of data points to 
convenient levels. Because of the reduced sampling period at the input to tlie 
second stage, a demodulation frequency resolution of better than 1 mHz can now ђе 
achieved. The FIR filter of Figure 3-8 is again applied, this time resulting in a 
passband of 12.5 milz (compared to the 500 mHz bandwidth available in the first 


stage). Finally, the data is once more decimated, this time bv a factor of 10. 


First Stage Second Stage 
Data Sample Rate: 250 Hz 6.25 Hz 
Passband: 0-250 шН? 0 - 6.25 mHz 
Passband Bandwidth: 500 mfz 12.5 mHz 
Nominal Passband Gain: 1.0 1.0 
Passband Ripple: 0.35 dB 0.35 dB 
Transition Band: 250 - 1500 mHz 6.25 - 37.5 mHz 
Pass plus Transition Bandwidth: 3.0 Hz 75 mHz 
Nominal Stopband Gain: 0.0 0.0 
Stopband Ripple: -47.75 dB -47.75 dB 


Table 3-VI: FIR Filter Characteristics 


Both stages of the preprocessing described above have been carried out on all 
24 channels of data simultaneously, resulting in preprocessed complex time series 
that are demodulated to within 1 mHz of the nominal signal frequency and then 
very closely filtered around that frequency. The final sample frequency 15 0.625 Hz 
(equivalent to a sample period of 1.6 sec), providing approximately 2000 data points 
per channel for a signal of 55 minute duration. As the correlation length of the 
second stage low pass filter is about 20 post-decimated points, tliese 2000 pommts 
represent roughly 100 independent degrees of freedom. Further details of the digttal 


signal processing involved can be found in reference [68]. 
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Figure 3-9: Preprocessing Output Time Series - Vertical Arras 
Sensor Magnitudes at 17.00 Hz (part L of 2) 
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Figure 3-9: Preprocessing Output Time Series- Vertical Array 
sensor Magnitudes at 47.00 Hz (part 2 of 2) 
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3.12 An Analysis of the Preprocessed Data 


A great deal of insight can be gained just by studying the preprocessed 
complex time series of the various vertical array channels. Figure 3-9 displays the 
log magnitude of the these time series for the 47 Hz tone from the the vertical array 
data set, while Figure 3-10 shows the corresponding phase traces for the same 
signal. Data from all eighteen vertical array sensors are included; the depths of the 
various sensors are indicated on the left. A 40 dB per inch scale is used on the log 
magnitude plots; the phase plots cover the range from -7 to z and display only the 
principal value of the phase. The time axis spans a range of 65 minutes (3900 sec). 
The data displaved here ın Figure 3-9 and subsequently ш 3-11 are both 
uncorrected for the none sensitivity, making the amplitude levels only 
relative. In Chapter 7, where the data set 15 studied in much more detail, all 
amplitude levels presented are absolute and referenced to 1 „Pa. Note that for the 
present displays, the decreased sensitivities of the vertical array hydrophioies 
located at 180 and 210 m cause their levels to appear 6 dB lower than they actually 


are. 


The signal begins to be seen about 100 seconds into the trace and has the 
expected duration of 3300 seconds, ending just before the end of of the trace. 1t Is 
easily distinguished on most of the channels; signal to noise ratios of about 30 dB (їп 
the 12.5 mHz bandwidth of the final preprocessing filter) are typically observed im 
the log magnitude traces. The general slope of the phase traces Is due to residual 
mismatch between the received frequency and the demodulator frequency. The 
negative value of the slope indicates that the demodulation frequency is greater 
than tlie received frequency, while the magnitude of the slope is about 1 cycle per 
hour, corresponding to a mismatch of 0.3 milz. Of tliis error, 0.1 mHz is due to the 


difference between the actual demodulator frequency and the nominal source 
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Figure 3-10: Preprocessing Output Time Series - Vertical Array 
Sensor Phases at 47.00 Iz (part 1 0f 2) 
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Figure 3-10: Preprocessing Output Time Series- Vertical Array 
Sensor P’-hiases at 47.00 Tlz (part 201 2) 
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frequency. The remaining 0.2 mHz is attributed to a small amount of Doppler shift 
caused by source-receiver motion, and corresponds to a relative opening range rate 
of about .01 kt. Тһе size of this result agrees well with the ice drift speeds 
experienced at the FRAM camp [78] (no drift speeds were available for the 
TRISTEN camp). The shght change tn phase slope towards the end of the signal 


indicates that the drift rates were not constant over the full duration of the signal. 


In general, the temporal stability of both the log magnitude traces and the 
phase traces is exceptional. What variations do exist may be broken into short term 
fluctuations that have time scales of minutes, and longer term fluctuations on the 
order of tens of minutes. That no significant fluctuations exist with time scales of 
seconds or smaller is evidenced by the demonstrated ability to coherently process 
the data in a 12.5 mHz bandwidth without loss of signal power or phase consistency. 
The short term fluctuations are almost completely negligible on traces of adequate 
signal to noise ratio, and generally become more pronounced on both tvpes of traces 
as the SNR = decreases. This supports the supposition that the short term 
fluctuations are probably due to background noise, since the other two likely causes. 
temporal channel variations and array movement, would not be expected to slow as 
much dependence on SNR. The cause of the longer term variations. particularly 
those in the log magnitude traces, is more problematic. One example is the deep 
fade exhibited on the 90 m hydrophone over the last third of the signal: other 
examples of smaller fades can be seen throughout the data. The cause of these Is 
not well understood. They are probably most attributable to channel variations 
over time, which can be generated in two ways: by actual temporal variations in the 
physical structure of the channel, and by vartations In souree-receiver range due to 
relative ice drift. It is possible that at least some of these variations are related to 


longer term array movement caused by variations m the effective current merdent 
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on the array. The change in the signal Doppler noted above is one indication that 
this current might not be constant over the full duration of the signal. Navigation 
data also indicates some changes in drift velocity over the length of the 
transmission. Theoretically, however, it is difficult to justify the signal level 


variations observed from just small changes in sensor offset. 


While the general quality of both sets of traces is excellent, some noticable 
exceptions occur. One example is the pair of hydrophones at 330 and 350 m. As 
above, the poor quality of these traces appears to be caused by the abnormally low 
signal to noise ratios exhibited on these hydrophones. After a study of all the data. 
it is difficult to support the hypothesis that these low received signal levels are, 
caused by instrumentation problems. First, Just the evidence of coherent (if noisy) 
phase traces with the proper average phase ramp indicates that these sensors are at 
least partially responding to the actual acoustic signal. Additionally. the channets 
exhibiting low signal levels tend to be closely grouped in depth. as seen here. 
Finally, the particular groups of low SNR sensors are frequency dependent: that ts. 
the particular hydrophones exhibiting low signal levels change when the frequency 
bemg transmitted changes. As an example, consider the log magnitude traces 
exhibited in Figure 3-11 and the associated phase traces illustrated in Figure 3-12. 
These are identical in format to the data displayed ın Figures 3-9 and 3-10, 
meluding the channels displayed, but are for the 17.75 Hz signal from the vertical 
array data set. Here it is the very lowest sensors in the water column, the two at 
870 and 960 m, which display the lack of signal to noise ratio. One must then 
conclude that these low SNR channels actually represent nulls m= the vertical 
structure of the signal field, caused. erther by. nulls in some dominant mode shape or 


bv the coherent mterference of multiple modes that might be present. 


The latter hypothesis is strongly supported by another observation that can be 
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Figure 3-11: Preprocessing Output Time Series- Vertical Array 
Sensor Magnitudes at 17.75 Hz (part 2 of 2) 
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made from this raw data. If a single dominant mode were present and the array 
were truly vertical, then all the phase traces would exhibit either identical absolute 
phase, or would be shifted by z after going through a null. An easy way to judge 
the relationship of the absolute phase between sensors is to note the point at which 
phase wrap around occurs; for two channels to have identical absolute phase, the 
wrap around must occur at the same point in the trace. The phase traces displayed 
in Figure 3-10 most certainly do not exhibit either identical absolute phase or 7 
phase shifts; moreover, it is very difficult to explain the phase shifts exhibited in 
terms of array tilt away from the vertical. In order to explain the shifts of just the 
top three channels, (those at 30, 60 and 90 m) in terms of array tilt, it Is necessary 
to assume that the local slope of the array near the surface is about 7.5 in the 
plane of propagation. about double the value shown in Figure 3-5 and well outside 
the range of reasonable array tilts. Therefore, ıt ıs нь to assume that the 


data displayed in Figures 3-9 and 3-10 include a number of reasonably coherent 


interfering modes. 


For comparison, consider the data displayed in Figures 3-11 and 3-12, which. 
as mentioned earlier, is the output of the preprocessing for the vertical array data 
taken at 17.75 Hz. Although the data at first glance is similar to the 47 Hz data 
just presented, a closer inspection reveals some very interesting differences. First, 
notice the presence of two strong events, one Just before the start of the signal and 
the other about 600 seconds later. These are known to be blasts from an airgun 
that was deployed near the FRAM ice camp. Note that the signal phase tracks 
through the second blast in a reasonable fashion. In general. the signal to поље 
ratio is slightly lower than that encountered earlier. and the short term variations, 
particularly those of the phase traces, reflect the increased relative noise level. The 


frequency mismateh with the demodulator 1s somewhat higher (about 1.1 mfz) but 
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is still well within both the bandwidth of the filter and the range of reasonable 
Doppler shifts for source-receiver motion. Additionally, the Doppler shift appears to 


be more constant over the signal duration than that observed previously. 


The most important difference between the two sets of data, however, is the 
way the absolute phases align in Figure 3-12. This is a good example of the 
absolute phase alignment discussed above, leading to the conclusion that a single 
mode dominates the vertical structure of the signal at this frequency. Here, the 
dominant mode must be the first, since no z phase shifts (corresponding to sign 
changes in the mode shape) are observed down the array. Anticipating the results 
of Chapter 4, the relative levels of the log magnitude plots as a function of depth 
also indicate a dominant first mode, particularly the low SNR levels of the bottom 
sensors in the water column (the 870 m and 960 m traces); the significant extent of 
the first mode at 17.75 Hz is only about 800 m. Note that the traces for the 
hydrophones at 180 and 210 m are offset by 6 dB from the remaining traces, due to 


their different sensitivities (as indicated in Table 3-III). 


The two sets of traces just presented are representative of the outputs of the 
preprocessing for both Ше horizontal and vertical array data sets. For all the 
signals detected, signal to noise ratios fall in the 15 to 40 dB range for the final 
filter bandwidth of 12.5 mllz, and both stable log maguitude and phase traces are 
observed. The short term variations (those on the order of minutes) on both type of 
traces appear to be related to background noise corruption, since their size appears 
to vary as the SNR observed at the hydrophone varies. Some longer term signal 
variations on scales of tens of minutes are observed: while their cause is not well 
understood, it is possible that array movement is partially responsible, although it is 
unlikely that this is the sole cause. Plis leads to the conclusion that the vertical 


array may be considered stationary for periods of up to 10 or 20 minutes: at can 
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probably not be considered constant over periods of an hour or longer, though. 


3.13 Summary 


The ЕКАМ IV Arctic Acoustic Experiment has been summarized, including 
descriptions of the general locale and arrangement, the scientific objectives, the 
equipment used, and the data set of interest in this thesis. The data presented here 
consists of a rich set of tonals in the frequency range from 5 to 71 Hz. These were 
broadcast from the low frequency source deployed at the TRISTEN ice camp. and 
then propagated roughly 250 km to the main FRAM ice camp. Here, they were 
received on both of the available sensor arrays, and were recorded in a digital 


format on a 24 channel acquisition system. 


The nature of the two arrays has been discussed in some detail. providing 
some insight from an engineering viewpoint into the non-ideal aspects of vertical 
arrays versus those of horizontal arrays. In general, it may be said that the process 
of designing and implementing a good vertical array is much more demanding than 
the equivalent development of a horizontal array. In investigating some of the 
possible difficulties, useful estimates of their various magnitudes for the FRAM 
arrays have been made. In particular, maximum sensor offsets of about 10 m and 
50 m can be expected for the FRAM IV horizontal and vertical arrays, respectively. 


For the vertical array, this corresponds to tilt angles of somewhat less than 5 


A detailed description of the preprocessing scheme used to compress the data 
set from more than 1500 Mbytes to a workable size has been provided. This scheme 
involves two stages, each consisting of quadrature demodulation, filtering with a 
very narrow FIR low pass filter, and decimation. The output time series resulting 
from the preprocessing scheme have been demodulated to within I mllz of the 


nominal signal frequency, filtered in à bandwidth of 12.5 mllz. and decimated to a 
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sample period of 1.6 sec. The last filter in the process provides a correlation length 
of about 20 post-decimated points, leaving approximately 100 degrees of freedom in 


a typical 55 minute signal. 


Two examples of the preprocessing output for the vertical array have been 
presented and analyzed. Signal to noise ratios of from 15 1o 40 dB are tvpical in 
the data set, due to the very narrow filtering of the data. The high SNR leads to 
nearly ideal log magnitude and phase traces. The demonstrated ability to 
coherently process the data in the very narrow bandwidth employed without loss of 
signal power or phase continuity indicates that no significant variations exist with 
time scales of less than about one minute. Short term fluctuations of one minute 
Order seen in these traces are clearly a function of SNR, indicating that thev are 
probably noise related. Longer term signal variations with periods of tens of 
minutes are observed; their cause is not well understood. Some, although not all, 
are likely due to array motion. A more important cause of these variations Is 
probably changes in source-receiver range over time due to relative drift between 
ice at the two locations. All the observed phase ramps are within the range of 


Doppler shifts that reasonable values of source-recelver motion would produce. 


The absohite phase alignment down the sensors of the vertical array provide a 
great deal of insight into the underlying modal structure. For the 47 Hz data of 
Figures 3-9 and 3-10, this phase alignment implies a. number of significant modes 
probably interfering coherently, since nulls in the vertical structure of the field are 
visible. For the 17.75 Hz data of Figure 3-11 and 3-12, however, a single inode, the 
first, appears to dominate the vertical structure. Ñ study of the pattern of recen ed 


апр] паде versus depth for the 17.75 Hz data supports the same conclusion. 





Сћарђег 4 


The Modal Structure of the Arctic 
Sound Channel 


The purpose of this chapter is to characterize the aspects of the central Arctic 
environment that impact the propagation of sound; to utilize that characterization 
to study the modal structure of the channel in which the FRAM IV experiment was 
conducted; and to examine the interaction between the modes and the source and 
receiver geometries present during the experiment. These results are directly 
needed for the work that follows, but also allow some more general observations to 
be made. Typically, such a study would be emploved only in the interpretation and 
modeling phases of an experiment. In the present effort, however, the modal 
decomposition techniques to be developed all require both the shapes and the 
associated horizontal phase speeds of the significant modes to be known a priori. 
Prediction of the channel's modes thus assnmes a more important role here, that of 


a necessary preliminary to the beamforming implementation. 


The requirement for prior knowledge of the local modal structure has 
important theoretical and practical consequences. Among the most significant 15 
that a whole range of modeling and prediction Issues now have a direct impact on 
the performance of the beamformer. From a theoretical viewpommt, questions of 
processing performance become dependent upon the assumptions made about the 
environment in which the array is deployed. Mere practical issues include the 
adequate ¿n situ measurement of the local sound speed profile and methods for 
estimating other needed environmental parameters. Techniques for computing the 


mode shapes and horizontal phase speeds must also be addressed. Since this 
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involves solution of the one-dimensional Helmholtz equation (2.7), usually by 
numerical or other approximate methods, the applicabilty, stability, and acenracy of 
the method chosen must be considered. The question of the channels gross 
temporal stability becomes important since dynamic updates of the modal data mav 
be necessary. These issues all affect beamforming performance to the extent that 
they generate mismatches between the assumed and actual modal structures: the 
mismatches then skew the results of the modal decomposition in proportion to their 


severity. 


The starting point for this investigation is a compilation of the environmental 
data necessary to describe the acoustic channel. The most important element is a 
knowledge of the average water column sound speed profile at the receiver: 
additional useful information includes descriptions of the bottom and overhead ice 
cover as well as an understanding of the spatial and temporal variations of the 
channel profile. The method of mode generation to be employed Is then discussed 
and the results of its application to the central Arctic sound channel summarized. 
Finally, these results are analyzed so that the impact of varions environmental 


features on the modal structure can be assessed. 


Coachman and Aagard [19] provides a good general reference for Arctic 
oceanography. Two works by Chen [15], [16] provide more detailed discussions of 
Arctic sound speed profiles directly applicable to the data set of interest. Some 
work on the modal structure of the central Arctic has been done previously: for 
comparison with the results shown here, the reader may consult Duckworth [29] or 
Polcari [67]. No significant variation of results is found among the three. Another 
good reference is Kutschale [48}, while Mellen [51] and Yang and Giellis [81]. both 


provide additional discussion. 
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Figure 4-1: Sound Speed Profile of the Water Column 
at the Receiving Arrays 


1500 


1500 





-92- 


4.1 A Representative Sound Speed Profile from the Central Arctic Ocean 


Figure 4-1 indicates the sound speed profile on which the following modal 
calculations are based. The left hand plot shows the detail of the upper 1000 ma, 
while the right hand plot shows the complete water column. This profile represents 
the best average estimate of sound speed at the FRAM receiver camp during mid- 


April, 1982, and is identical to the profile presented in Polcari [67]. 


As is found throughout the Arctic, the sound speed is essentially a monotonic 
increasing function of depth. Within the uppermost regions (down to about 150 m). 
the controlling mechanism is not the water temperature, as Is typically the case for 
most open ocean areas, but rather the salinity gradient generated by the presence of 
the overhead ice canopy [15]. The variation in salinity provides a change of about 5 
m/sec in sound speed across the first 100 m. By comparison, temperature variations 
account for about 1 m/sec of change in this zone. Below the 150 m point, the 
temperature takes over as the dominant factor, resulting in an almost isovelocity 
profile to roughly 700 m. The deep regions (below 800 m) are primarily pressure 


controlled, as is the case with deep ocean profiles in more temperate climates, 


The net effect of the combination of salinity and temperature variations Is to 
produce a very strong surface duct encompassing about the upper 200 m of the 
water column. As will be shown, this surface duct plays a significant role in 
determining both the modal structure of the Arctic channel and the nature of the 
propagation there. An average slope of the deeper portion of the central Arctie 
profile (that below 200 m) can be taken as 0.013 (m/sec)/m, not far different from 
the sound speed pressure dependence of 0.016 (m/sec)/m. By contrast, the mean 


eradient in the surface duct is 0.13 (m/sec)/im. about ten times larger. 


l'ieure. 1-2 details the sound speed profile used to model the Aretic bottom, 
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Figure 4-2: Assumed Bottom Characteristics at the Recening -Arrays 


И 
The main features here have been chosen based on the composite profile presented 
in Duckworth and Baggeroer [30]. The bottom is modeled as four layers. A 400 11 
sediment layer of density 2.2 g/cm’ extends to 4200 m from the assumed bottom 
depth of 3800 m, reaching a sound speed of 2.3 km/sec. Two layers are then used 
to transition at a depth of 5900 m to a soft basement capable of propagating sound. 
The sound speed of this basement is assumed to be 6.0 km/sec, while its density is 


s 3 
assumed to be 2.9 g/cm”. 
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Figure 4-3: Historical Bathymetry along the TRISTEN/FRAM 
Propagation Path (taken from [11]) 


[Note different range and depth scales] 
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4.2 Temporal and Spatial Sound Speed Variations 


The sound speed profile of the central Arctic has been shown to be unusually 
stable in both time and space [16]. This is due to a number of factors, the most 
significant being the geographical nature of the area and the vear-round continuous 
ice canopy. The isolation of the central Arctic limits variations due to external 
influences, while the ice cover insulates the ocean below from the solar and wave 
effects found elsewhere. Figure 4-1 can therefore be taken as representative of the 
profile that would be encountered anywhere along the transmission path throughout 
the length of the experiment. Some range variation mn the sound speed profile can 


-be expected due to large scale oceanographic features such as eddy systems, though. 


Although the water column is spatially well behaved, the same EET 
cannot be drawn of the bottom. Figure 4-3 indicates the bathymetry along tlie 
propagation path, based on charted data (in situ measurements were not available). 
The TRISTEN source camp was located almost over top of the Mid-Arctic Ridge. 
Which is easily identified in the figure from the central trench structure at a range 
of 40 km. The series of peaks on either side of this trench extend to the 2500 m 
depth range. The severe range dependence of the bottom significantly affects the 
nature of the propagation, as has been graphically illustrated by the comparison of 
shot data with that taken from earlier experiments (3). Any mode that interacts 
significantly with the boitom in the region of the ridge can be expected to be 
attenuated to such an extent that its contribution to the total sound field at the 
receiver would be smal; thus, modes having turning pomts much below 2000 ın are 


probably negligible for the TRISTEN/FRAM propagation path. 


(Da 
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4.3 The Acoustic Impact of the Ice Cover 


Much of the Arctic’s uniqueness from an acoustic viewpoint is due to its 
relatively continuous year-round ice canopy. Both ice camps were located well 
within the zone of continuous ice coverage. At times (especially in the late 
summer}, significant regions of open water can occur even in the area nominally 
designated as continuously covered. However, this effect is not overly significant 
during the early spring. The assumption of continuous ice cover over the full extent 


of the propagation path is therefore reasonable. 


The influence of the ice canopy on the propagation path is a significant effect 
which is not yet well understood. An excellent reference that reflects the current 
thought on the subject is a recent memorandum by DiNapoli and Mellen [25]. The 
primary effect of the ice cover at low frequencies is to increase the effective 
attenuation coefficient by about two orders of magnitude over valnes found in open 
water. Prediction of such a loss falls into the general category of rough surface 
scattering. A number of different models of varving complexity may be applicable. 
ranging from a simple rough free surface to a full solid laver model. including shear 
wave effects. To date, none of tlie models investigated shows acceptable agreeinent 
with experimental data. The most important parameters needed to describe the ice 
canopy are the mean ice draft and the rms roughness, correlation length. and power 
spectrum of the random surface that forms the ice-water mterface. Typical values 
reported by DiNapoli and Mellen for these parameters are a mean draft of t m. with 
an rms roughness of 2 m and a correlation length of 1£8 m. These results are in 
reasonable agreement. with earlier work published by Wadhams [83]. DiNapoli and 
Mellen also provide an analytical fit. to. the. experimentally measured ice-water 


interface power spectrum. They give the two-dimenstonal power speetrum as 
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assuming an isotropic random surface. Here x is the two-dimensional wavenumber 


vector, c ls the rms roughness of the surface, and L is the correlation length. 


In this thesis, the ice canopy is modeled by a rough free surface at the ice- 
water interface, more for simplicity than any other reason. While the roughness has 
an important direct effeet on modal attenuation, its Impact on the gross properties 
of the modes (the mode shapes and their horizontal phase speeds) is determined by 
its effect on the mean channel profile. Within the limits of the model, the position 
of the mean free surface is displaced by the amount of the mean ice draft (1-1). 
which is probably insignificant when compared to other sources of depth error in 
the profile. Ice effects have therefore been ignored in the sound speed profile used 


to generate the mode parameters. 


4.4 Solutions of the Helmholtz Equation for the Mode Shapes 


To convert the environmental information presented above into mode shape 
and phase speed data, one must solve the one-dimensional Helmholtz equation 
eigenvalue problem (2.7) for the specified sound speed profile at each desired 
frequency. The form of equation. (2.7) assumes that all lavers, meluding the 
propagating basement, are fluids; shear wave propagation in the bottom Is mplicitly 
ignored. In the sequel, it is shown that neither the mode shapes nor their associated 
horizontal phase speeds are significantly influenced by the bottom characteristics. 
Thus, the neglection of bottom shear effects Is a reasonable approximation. АН 
attenuation effects are also ignored here, as they do not affect the shapes and phase 


speeds of the modes to any great extent. 
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One method for finding the mode shapes and phase speeds is to employ the 
WKB approximation in solving the Helmholtz equation. This approach is 
questionable at the low frequencies of Interest here. Another technique involves the 
use of propagator matrices, but this method Is more applicable to a discretely 
layered channel model than it is to the continously varying profile of Figure 1-1. 
The most accurate method for the present circumstances Is a shooting technique, 
involving direct numerical integration of equation (2.7). The integration is 
initialized to meet one of the boundary conditions, and the assumed eigenvalue is 
varied until the other boundary condition is met. Lavers of varying density can be 
accommodated through the proper application. of. cross-boundary— continuity 
conditions (given in Figure 2-1 as the propagating bottom boundary conditions). 
The mode's horizontal phase speed is recovered from the eigenvalue, while. the 
eigenfunction, when properly normalized to meet equation (2.9), vields the mode 


shape. 


The main drawback of this approach is the instability of. the. numerical 
integration in the evanescent regions of the mode shape. Since equation (2.7) is 
second order, it supports two possible solutions, one exponentially increasing and 
the other exponentially decreasing in the direction of integration. Numerical errors 
that excite the increasing solution can quickly donate the desired solution, even 
though they are insignificant when introduced. The Arctic profile, with = its 
monotonic Increasing sound speed, is particularly susceptible to this effect, since the 


mode shapes in this instance typically exhibit extensive evanescent regions, 


To circumvent tlus problem, a modified integration method developed by 
Baggeroer [1] has been employed in generating the mode shapes presented here. 
The Helmholtz equation ts first cast into a state variable form of second order. 


Treating the state variables as a rectangular Cartestan pair, the problem is then 
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transformed to a magnitude/phase pair of equivalent polar coordinate variables. 
The differential equation for the magnitude variable decouples from the equation for 
the phase variable and the boundary conditions take on rather simple forms, so that 
only the first order differential equation for the phase need be solved via shooting to 
obtain the eigenvalue. Then. given the phase solution, the mode shape Is recovered 
by integrating for the magnitude and converting back to the original variables. The 
development by Baggeroer includes a series of conclusions about the proper 
directions of integration for numerical stability. This approach has also been 


successfully emploved by Duckworth [29]. 


To ensure the greatest possible accuracy, the Baggeroer approach has been 
implemented here with a series of fourth order Adams-Bashforth predictor-corrector 
integrators [12]. А depth grid of 5000 points is employed for the integrations, 
providing a sample interval of 1.2 m. The resulting mode shapes have been 
normalized for the pressure ‚field (rather than the velocity potential field) in 
accordance with equation (2.9). As a test of the internal consistency of the 
algorithm, the orthogonality integral of equation (2.8) has been computed for 
different mode pair combinations at several frequencies of interest. Typical 
integration results are of the order of 10°*, with the worst errors generally occurring 
either when the pair consists of neighboring modes or when the first mode is one of 
the pair. However, even in these conditions results as large as 107% are obtained 


only rarelv. 
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Figure 4-4: Shape of the First Mode at Varions Frequencies 


4.5 An Analysis of the FRAM IV Modal Structure 


All modes with turning points above the nominal bottom depth of 3800 m 
have been computed at the frequencies available in the vertical array data set (the 
nominal sound speed at 3800 m is 1508 m/sec). Table t-I indicates the number of 
modes found at each frequency, as well as the number with turning points above 
2500 m: as discussed above, the modes contributing significantly to the sound. field 
at the receiver are all probably tncluded in this last group. Figures 1-1 through 1-7 
display the shapes of the first four modes at some of the frequencies for whieh data 


Is available. The frequencies chosen for display were selected so as to span those 
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Figure 4-5: Shape of the Second Mode at Various Frequencies 


frequencies of interest that are not shown. 


4.5.1 Mode Shape Analysis 

The extent of a particular mode is inversely related to frequency. Higher 
order modes look similar to the fourth, but extend further into the water column 
with a greater number of sign changes. Their general shape is not far different from 
that of Airy functions. In all cases, the absolute maximum of the mode shape 
occurs in the half cycle immediately above the mode turning point. 

As both the source and the horizontal reeeiving array were deployed to a 


depth of 91 m. the values of the various mode shapes there are of particular 
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Figure 4-6: Shape of the Third Mode at Various Frequencies 


interest. Figure 4-8 indicates the levels of the mode shapes in dB at 91 m for the 
first 10 modes at 15, 30, 47, and 71 Hz. As indicated in Chapter 2. it is these levels 
which control the distribution of source energy in the various modes. Figure 4-8 
shows that, at its 91 m deployment depth, the TRISTEN source drives the the first 
mode much better than any other mode. The higher order modes, while excited less 
efficiently than the first, are all driven about equally; moreover, the level of this 
plateau appears to be insensitive to frequency. The #7 Hz curve is something of an 
anomaly, in that at this frequency modes 2 through 6 all exhibit nulls vers close to 
the source depth (see Figures t-f through 4-7): this accounts for the unusually. low 


levels of modes 2. 3 and 4. However, the higher order modes at 10 Hz still appear 
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Figure 4-7: Shape of the Fourth Mode at Various Frequencies 


to plateau at about the same levels as those of other frequencies. 


Figure 4-8 may also be interpreted in a different context. Because the 
horizontal array at the FRAM receiving camp was deployed to the same 91 m depth 
as the source, the discussion above apphes equally well to the relative sensitivity of 


the horizontal array elements to energy carried by the various modes, 


From Figure 1-8, it can be seen that tlie advantage in excitation that the first 
mode possesses over higher modes is a funetion of frequency. Winle this difference 
Is quite large in the two middle frequencies, tt i8 much smaller at tlie edees of tlie 
frequency. baud of interest. To highlight this effect; Figure 4-9 has been included. 


This figure displays the relative advantage of the first mode over its neighboring 
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Frequency 3800 m 2500 m 
15.00 Hz 14 T 
7105.47 17 9 
20.00 Hz 19 10 
23.50 Hz 22 2 
30.00 Hz 28 15 
29 25 H2 33 17 
35.50 Hz 34 18 
47.00 Hz 45 2 
5325 H-2 51 26 
55.00 Hz 93 о 
71.00 Hz 65 92 


Table 4-1: Number of Modes with Turning Points above 
3800 m and 2500 m 


modes as a function of frequency. The relative advantage formula given in the 
figure is a crude method of quantifying how much the first mode extends above 
neighboring mode shape levels in Figure 4-8. Note the strong 30 dB peak at about 
45 Hz in Figure 4-9. About 10 dB of the plotted value near the peak ts due to the 
uncharacteristicly low levels of modes 2 through 5 (the modes used in the average 
indicated ш the figure): however. a reduction of the central peak by even 10 dB 
does not change the fundamental nature of the plot. The explanation of the roll off 
at the extremes of the frequency range is straightforward. At the lower frequencies, 
it is the mode shape level for the first mode that is decreasing. rather than those of 
higher modes increasing. At the high frequency end of the spectrum (above 70 lz). 
the mode shape levels for modes 2 and 3 finally begin to approach that of the first 
mode, again causing a decrease in the first mode relative advantage. Both of these 
effects can also be seen in Figure 4-8. 

Two important conclusions must be drawn from Figures 4-8 and 4-9. First. 
the deployment of the TRISTEN source heavily favored first mode excitation. by as 


mueh as 15 dB or more in the middle frequencies of the range. H (8 not until 71 И. 
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Figure 4-9: Relative Advantage of the First Mode over the Next 
Four Modes at 91 m versus Frequency 


the highest frequency available in the vertical array data set, that the second and 
third modes begin to reach excitation levels equivalent with that of the first mode. 
Second, and equally important, because the horizontal array at FRAM was deployed 
to the same depth as the source, it was also heavily biased towards the first mode. 
so much so that data taken from the horizontal array hydrophones must ђе 
considered to be almost completely dominated by first mode characteristics, While 
this effect was anticipated, its magnitude is rather surprising, especially for sound 


propagation from the TRISTEN source to the horizontal array. where it comes into 
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play at both the source and receiver. 


4.5.2 Phase Speed Analysis 


While an analysis of the mode shapes provides information about relative 
levels to which the source excites the various modes, an analysis of the horizontal 
phase speeds provides much more insight into the effects of the various parts of the 
sound speed profile on the modal structure. This is because the phase speeds are 
much more sensitive to variations in the profile than the mode shapes. Figure 4-10 
provides plots of the phase speeds for the first three modes of the central Arctic 
channel in the 0 to 80 Hz frequency range. In general, these plots consist of three 
regions: a high frequency region, Where the modes are effectively trapped їп the 
surface duct; a transition region; and a low frequency region where the modal 
characteristics are dominated by the deeper portions of the profile. To emphasize 
this effect, modal phase speed results obtained from the two sound speed profiles of 
Figure 4-11 are also plotted. In the left hand profile, the effect of the surface duct 
has been eliminated by extending the profile from the 200 m point to the surface at 
the gradient existing just below that point. In the right hand profile, the surface 
duct has been retained and extended all the way to the bottom of the water columna 
in a similar fashion. The gradient used here to extend the profile is the slope 


existing just above the 200 m point. 


Given the results of the previous section that indicate the dominant role it 
plays in the transmission of sound from the TRISTTEN source, a study of the first 
mode is of special interest here. As can be seen in Figure 4-10, below 19 Hz the 
effect of the surface duct on first mode propagation i small... From 10 t0 30 Elz. the 
mode is in transition into the surface duct, while above 30 Elz it can be considered 


to be trapped im the surface duct. As the mode falls into. the surface duct. a 
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significant reduction of phase speed (from about 1460 m/sec to about 1150 m/sec) is 
evident. Examination of the mode shape as a function of frequency (Figure 4-1) 
supports this conclusion. By comparison, the second mode does not begin its 
transition until about 35 Hz and is not completely trapped in the surface duct until 
a frequency of almost 60 Hz is reached. The center of the third mode transition 
region is about 65 Hz, but it does not become completely trapped until well above 
80 Hz. As the various modes fall into the surface auct, the spread in horizontal 


phase speed from those not yet trapped becomes much more pronounced. 


Phase speed variations have also been used to Investigate the sensitivity of the 
results to the bottom characteristics. To accomplish this. the propagating bottom of 
Figure 4-2 is replaced by a hard bottom and the modes recomputed at 23.5. 35.25. 
and 47 Hz. .Measurable differences between the resultant phase speeds and those 
computed earlier are first detected at modes 19, 30. and 41, respectively. All the 
modes exhibiting sensitivity to the bottom structure have corresponding horizontal 
phase speeds in excess of 1500 m/sec, or, equivalently, possess turning points within 
500 m of the bottom. Thus, they are excluded from the set of modes thought to 
most significantly contribute to propagation in the channel (those with turning 


points above 2500m). 


4.5.3 General Observations 


Comparing the data from Figure 4-9 with that of Figure 4-10. ome observes 
that the frequency region in which the first mode is most dominant (between about 
20 and 60 Hz) coincides with it being the only mode trapped in the surface duet. [t 
can then be concluded that any source or receiver deployed іш the duct i very well 
tuned to tlie modes trapped there. but only poorly matched to those that extend 


outside 1t. In such a case. the environmental effects that most greath affect the 
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sound propagation are those that are located at the top of the water column. This 
is Why parameters such as basement characteristics and bottom topography have 


relatively little impact on the problem at hand. 


In the upper layer of the Arctic, the only two significant environmental 
considerations are the ice cover and range variations in the sound speed profile. 
The former is known to have a pronounced impact on the sound propagation. as 
mentioned earlier. At the low frequencies of Interest here, the effect of the latter i 
probably small by comparison for a number of reasons. Surface duct range 
variations influence propagation in the duct in three ways. First. they provide 
density and refractive index contrasts that. scatter. tlie. incident. sound. energy, 
effectively increasing the attenuation. However, this effect is negligible when 
compared to the scattering caused by the ice, which provides a much better acoustic 
contrast. Second, they affect the way that the source and receiver couple into the 
channel by altering the shape of the modes that are trapped in the duct. Finally, 
thev affect the distribution of modal energy by generating mode coupling. A review 
of the coupling coefficients (equations (2.18) and (2.19)), though. reveals that. their 
size Is primarily a function of the net change encountered in. the. mode shapes. 
Thus, both of these last two effects are controlled by the amount of mode shape 
Variation that can be generated. In general. it can be said that the shapes of the 
lowest order modes (the ones which the source primarily excites) are very Insensitive 
to small variations in the souud speed profile. These mode shapes tend to be more 
dependent on the larger scale trends in the profile than on its details. A change of 
as much as several meters per second in sound speed at any point in the duct can 
therefore be expected to have only a very small affect on the шо фе shape. 
Баг А. variations in duct depth of as much as 20 or 30 m are not overly 


significant. Since variations m sound speed much beyond these values are unlikely, 
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it is hard to argue that range variations in the sound speed profile have an 


important influence on low frequency Arctic surface duct propagation. 


The above conclusions are very dependent on the location of the source in 
depth. If the source is positioned below the surface duct, a much different situation 
is encountered. Under these conditions, a large spectrum of higher order modes are 
excited at the expense of the modes trapped in the surface duct. The bottom 
characteristics can then be expected to play a much more important role in 


controlling the overall propagation picture. 
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Figure 4-12: Deepest RSR Rav for the TIISTEN/FRAM Transmission Path 


[Note different range and depth scales] 
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4.6 The Deepest RSR Ray Paths of the TRISTEN/FRAM Channel 


As concluded in Polcari [67], it is at times easier to characterize the 
propagation in terms of a ray representation than it is to describe them in terms of 
modes. This is especially true of the deepest diving RSR гауз, which correspond to 
coherent sums of large numbers of higher order modes. all of which are closely 
spaced in horizontal phase speed. Figure 4-12 is adapted from [67] to provide a 
general feel for the characteristics of the deepest RSR ray paths for the 
TRISTEN/FRAM transmission path. These rays are seen to exhibit: horizontal 


phase speeds in the 1480 to 1495 m/sec range. 


4.7 Summary 


The environmental aspects of the central Arctic sound channel and their effect 
on the modal structure encountered during the FRAM IV Experiment have been 
investigated. The water column sound speed profile is presented in Figure 4-1, and 
can be assumed to be reasonably stable over both time and space. The bottom 
topography creates major deterministic range dependence in the channel. The 
overhead ice canopy most greatly effects the attenuation of modal energy. but has 
little influence on the gross aspects of modal structure. the mode shapes and their 
horizontal phase speeds. The exact mechanism for the attenuation Increase ік по! 


well understood. 


The sound speed profile consists of three parts: a strong surface duct 
extending to about 200 m, which is unique in that it is generated by vartations In 
salinity rather than temperature; a monotonte increasing gradient below the surface 
diet that extends to the bottom at 38500 m, and which is generally pressure 


controlled: and a basement for which some expertmental data on the geophysteal 
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structure is available. The profile is typical of that generally encountered in the 
basins of the eastern region of the central Arctic. Such a profile seriously increases 
the possibility of stability problems when solving the Helmholtz equation 
numerically in order to obtain the mode shapes; special efforts have been made here 


to avoid this problem. 


An examination of the mode shapes reveals that. because of the source and 
receiver geometries involved, first mode propagation plays a dominant role in the 
FRAM IV data set. It is exhibited in the relative excitation levels of the various 
modes by the TRISTEN source, where the first mode is often driven at levels more 
than 15 dB above those of other modes. Also, the sensitivity of the horizontal 
receiving array to energy propagated’ in the first mode exceeds its sensitivity to 
other modes by equal amounts. In general, the relative mode shape levels for higher 
order modes appear to reach a plateau whose level Is not very sensitive to the 


frequency being considered. The shapes themselves are similar to Airy functions. 


By studying variations In the associated modal horizontal phase speeds as a 
function of frequency, the effect of the various parts of the central Aretic sound 
speed profile on the modal structure has been deduced. Each mode can be seen to 
transition from the deeper part of the profile to the surface duct. with. increasing 
frequency; for the first mode, this transition starts at roughly 10 Hz and rm 
completed by about 30 Hz. The dominance of the first mode is intimately linked to 
this process. The bottom characteristics do not appear to influence the modes 
thought to contribute most heavily to sound propagation through the channel. 
Finally, the deepest propagating RSR ray paths occur at horizontal phase speeds in 
the 1180 to 1195 m/sec range: these represent coherent sums of large numbers of 


In general. sensors deploved in the Arctie surface duet couple very well to 
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those modes that are trapped in the duct, and only poorly to any mode extending 
beyond it. This makes environmental effects that exist in the uppermost laver of 
the Arctic much more significant than effects located at other depths for any 
situation involving surface duct propagation. For this case, the impact of variations 
in the sound speed profile appear to be small compared to that of the ice canopy, 
since the mode shapes of importance are relatively insensitive to small changes in 
the profile. Sources deployed below the duct can be expected to exhibit 


significantly different propagation regimes. 
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Сћарђег 5 


Modal Beamforming: 


Theory and Sımulation 


In the three previous chapters, the preliminaries necessary for "modal 
beamforming" have been reviewed. The source model has been discussed: the 
available receiving arrays have been studied in detail; aud the modal structure of 
the sound channel has been carefully analyzed. This chapter draws npon elements 
of all of the three previous chapters to study the various mathematical algorithms 
that might be employed in the modal decomposition of a narrowband sound field at 
a receiving array. Such decompositions are close analogs of the methods used in 
traditional plane wave beamforming; from a theoretical viewpoint. the fundamental 
difference is the choice of the basis. set on. which the observed sound field ts 
decomposed. To emphasize the similarity, the term “modal beamforming" is used 


here to refer to the modal decomposition process. 


In plane wave beamforming, one attempts to measure the spatial structure of 
a narrowband signal field by estimating its three-dimensional spatial Fourier 
transform, or wavenumber spectrum. In modal beamforniiug. the estimate Is of a 
trausforin tliat is only partially Fourier in nature; a modal transform ts Instead used 
m the vertical dimension. From a practical viewpoint, each of these choices 
possesses distinct advantages and disadvantages. The pure plane wave 
representation is very simple. yet very general, requiring only limited. à prior: 
knowledge of the characteristics of the propagation channel. This makes it a vers 
robust. chaunel-independent procedure. Additionally. the extensive bods of 


experience from other disciplines on the use (and misuse? of the Fourier transform is 
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directly applicable in the study of performance issues. Finally, the existence of 
efficient computational algorithms for evaluation of the transform (1.e.. the FET), 


enhance the desirability of such an approach. 


Unlike Fourier decomposition, modal beamforming requires intimate prior 
knowledge of the propagation channel, so that the modal structure needed to 
implement the algorithms can be determined. In many respects, it may be thought 
of as a beamforming approach that is tuned to the propagation channel in which it 
is implemented. As discussed in Chapter 6, modal beamforming is also less forgiving 
of errors in sensor position, particularly when these errors are generated by tilt in a 
vertical array. The added sensitivity and decreased robustness of a modal approach 
are balanced by a number of benefits, however. The most important of these is the 
very efficient representation of the received signal field provided by the modal 
decomposition. This efficiency has both mathematical and physical aspects. 
Mathematically speaking, the modal representation provides a countably infinite 
basis set as opposed to the infinitely. dense set of plane waves: further. this modal 
basis set can often be limited to a small number of significantly contributing modes 
with relatively simple arguments (such as those presented іш Chapter 2). From a 
physical viewpoint, the modal approach is in some sense the natural decomposition 
to employ, because It represents a direct eigenfunction representation of the sound 


field. 


The chapter starts with the derivation of modal beamformers in a general 
sense from the principles of least squares estimation theory. The modal equivalents 
of two popular plane wave beamforming techniques are developed: spatial matched 
filtering: aud maximun liklhood (MEA!) algorithm. Each of these methods can be 
implemented in either a single beam or a multiple beam variant. providing a total of 


four different appronehs. A comparison of the four is made to highlight therr 
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theoretical differences. Two issues then remain. The first is the question of relative 
performance, which, because of its complexity, ıs studied in some detail. The other. 
which is closely related, involves the proper selection of a modal source model. 
Both of these issues are investigated using simulated data in addition to theoretical 


results. 


The theoretical developments presented here are covered in the literature to 
varying degrees. The fundamental techniques of least squares estimation theory 
have a rich history, and any number of good standard references are available: one 
such example is Lawson and Hanson [19]. References which deal with direct modal 
amplitude estimation techniques are limited. Hinich [12] first proposed the multiple 
beam MLM algorithm derived here, but did not applied it to actual data. Shang, 
Clay, and Wang [75] develop the multiple beam least squares algorithm and apply it 
to data generated in a laboratory setting. Signal processing aspects of thé algorithm 
(particularly performance issues) are not discussed, however. While Ferris [33] uses 
a single beam least squares algorithm to make modal amplitude estimates from field 
data, his experiment is purposely structured to avoid the full modal decomposition 
issue. A similar approach is taken true of the tow tank experiment of Hobnek. 


Tindle. and Muir [4 1]. 


More research has been completed in two related areas. One of these areas Is 
the study of the plane wave equivalents to the algorithms Investigated here. 
Baggeroer [2] provides a comprehensive reference on plane wave beamforming, 
including the MLM approach, which was originally introduced in the classic paper 
by Capon [10]. Cox [21] analyzes. the performance of MLM plane wave 
beamformers in detail. Schweppe [71] studies the concept of multiple beam least 
squares beamforming. Finally. a recent paper by Hayes. Ripa., and Mangum 


[38] provides further insight into the relationship between the least squares and 
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MLM approaches. 


Related work, much of it recent, can also be found in the area that has come 
to be known as “matched field processing". The techniques of interest in this field 
are not directly involved with modal decomposition; rather, they use prior 
predictions of the total signal field to generate direct estimates of source 
parameters, such as range and depth. However, their performance is fundamentally 
linked to the modal characteristics of the sound channel in which they are 
employed, making them of much interest in the present effort. Clay [17] first 
proposed the use of matched field processing; this work also forms the basis for 
Chapter 7 of Tolstoy and Clay [79] and Appendix .\5 of Clay and Medwin [18]. Ina 
more recent paper, Bucker [9] proposes a similar but more robust scheme. Hinich 
develops an MLM algorithm for source depth estimation [12] and studies the 
problem of the optimal array locations for making such estimates [13]. Heitmever. 
Moseley, and Fizell [39] have completed detailed simmlations of matehed field 
processing in a Pekeris waveguide, while Porter, Dicus, and Fizell [70] have done 
similar work for a deep Pacific channel.  Fizell and Wales [31] have had some 
preliminary success applving the techniqne to actual field data from the Aretic 
Ocean. Finally, Shang, Clay, and Wang [75] have proposed a related (but 


somewhat more simplistic) approach for passive source ranging. 


5.1 Development of the Least Squares Modal Beamformer 


The theoretical development of the least squares modal beamformer is 
straightforward. Let the № х I complex vector P. denote the complex amplitude 
observed at each of the V sensors in the array after completion of the two steps of 
quadrature demodulation and filtering described in Section 3.11. For the moment, 


assume that Pois a constant: the case of temporally varying complex amplitudes ds 
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taken up in the sequel. With the exception of the assumption of time invariance, P 
is then identical to p(f) in equation (1.5). Because P is a complex vector. its 


elements include information about both the magnitude and phase of the signal. 


The least squares beamformer computes the set of complex modal amplitudes 
which minimizes the total square error between the received pressure field and the 
modal source description developed ın Chapter 2. This error may be quantified by 


defining an error vector 


e=P-EÁ. (5.1) 


The second term in this equation represents the assumed source model summarized 
in equations (2.33) through (2.36). А is an .M X 1 time invariant complex. eolumn 
vector representing the resulting magnitude and phase estimates for the amplitude 
of the M modes that are being modeled. Note that the choice of M ix a modeling 
issue Of some complexity Whose resolution is left to the user. The resultant 
performance of the processing is quite sensitive to the proper selection of this 


parameter, as Investigated in Section 5.10. The N X M steering matrix E is the 


same às that described in equation (1.4). 
The least squares requirement is that the real scalar 
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be minimized över all possible choices of the A vector. The minimization Is easily 
performed by variational caleulus. Setting the first variation of equation (5.2) to 


zero. one obtains 





-121- 


ll 
[> 
| 
үт | 
= 
іт 


= (АТЕТЕ—Р' Е)ОА + дАЗ (Е' 


which requires that 


ТЕА=Е*Р. (5.4) 


[En | 


Assuming the Af X M matrix ETE to be non-singular. equation (5.4) may be 


solved directly to find the least squares modal estimator 


~ 


A (5.5) 


Such an assumption requires that M not be larger than N. This requirement may 
be circumvented by identifying the linear transformation in equation (5.5) as the 
standard matrix pseudo-inverse; the techniques for extending this operator to the 
case Where Af is larger than № аге quite standard [19]. Such an approach has a 
number of drawbacks, though; the most important being that the resulting modal 
amplitude estimates are no longer linearly independent. Therefore. it is not 
considered here. From a practical standpoint, the assumption of non-singularity 
also raises the issue of the numerical precision needed to accurately compute the 
Inverse, a question intimately linked to the proper selection of M. 

An investigation of the second variation demonstrates that the solution results 
in a minimum total square error if the E* E matrix is positive definite (which it 
must be if it is non-singular), 


АЕ ЕСА > 0. (5.6) 


The residual total square error for the beamformer can be obtained by substituting 


equation (5.5) back into equation (5.2) and simplifying. which results in 
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? min = 


Based on equation (1.6), an estimate of the M X M modal cross-coherence 


matrix may be computed as 


Sy =E [A 4*] (5.3) 
— (E* EJ! E* S; E(E* EJ! 
where 
Se Ee | (5.9) 


is the more traditional N X N sensor cross-coherence matrix. 


If the vector Р is actually time invariant, then the modal cross-coherence 
estimator Is redundant, since it provides no more information than the amplitude 
estimate given in equation (5.5). On a practical basis, though, this vector fluctuates 
randomly over time. The application of equation (5.5) to such a problem can be 
justified only on an instantaneous basis. One сап envision using a time slice of data 
to form P, which is then transformed into the complex modal amplitude estimates 
for that instant of time. Under such conditions, equation (5.5) provides only an 
estimate of the random complex amplitude time series of each mode, since it 
includes no averaging. By contrast, equation (5.8) vields information that ts both 
more useful and more stable, since it instead estimates the second order spatial 
statistics of the vector random process. Since the expected value operator Is almost 
always implemented by temporal averaging rather than ensemble averaging, the 
modal cross-coherence estunator may be thought of as averagme a number of 


estimates of the matrix А А“. Where A is determined on an instantaneous basis 
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from equation (5.5). Alternatively, one can. view the averaging taking place in the 
computation of the sensor cross-coherence matrix, with equation (5.8) providing the 
transformation to the equivalent modal cross-coherence matrix. The two views are 


identical because the processor is linear. 


5.2 Single Beam versus Multiple Beam Beamformers 


As developed, equations (5.5), (5.7), and (5.8) define the multiple beam least 
squares modal beamformer. This beamformer ts exactly analogous to Schweppe's 
multiple beam decoupled processor [7-4] for plane wave beamforming. with the sole 
exception that the construction of steering matrix E is different. While it Is possible 
to implement multiple beam beamformers directly, a different approach is more 
commonly used, especially in azimuthal plane wave beamforming. A single plane 
wave source model is chosen, so that M = 1. The resulting algorithm i5 then 
evaluated repeatedly for various steering vectors, each of which corresponds to a 
plane wave propagating in a different direction. One can develop an equivalent 
approach for modal beamforming by again selecting M to be l, and then making 
successive independent estimates by indexing the steering vector over the mode set 
of interest. This implementation is called the single beam least squares beamformer 


in this thesis. 


Such an approach has several advantages. Most importantly, it eliminates the 
need for the user to select a valne for M a priori. The resulting processor is simple 
and intuitive, forming the modal equivalent of a conventional beamformer. 
Alternatively, mt may be considered to be a spatial implementation of the matched 
filter concept. Finally. the requirement that A be less than NA is relaxed. so that 


one is not constrained in the number of different modes that can be estimated. 


As long us the mode set of interest is smaller than the number of sensors 
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available in the array, the processing can be implemented by either approach. 
Therefore, it is important to understand the relationship of the two methods. This 
can be readily demonstrated with a simple example. Consider a signal which 
consists only of a single mode, with no background noise. Let the steering vector 
associated with the mode forming the signal be denoted by E, and let a, be the 


mode's complex amplitude. In this situation, the sensor cross-coherence matrix Is 


"E Et. (5.10) 


Assume that the single beam least squares beamformer Is steered to a mode 
different than the one present in the signal. Let the steering vector associated with 
this mode be E,. From equation (5.8), the resulting single beam least squares 


estimate for the squared magnitude of this mode is easilv found to be 


€) 


| [E E |" 
a,” = (а, 


ILE. (5.11) 
Е |" 

Note that when the beamformer is steered to the mode actually present m the 

signal, the estimate in equation (5.11) Is exact. When steered to any other mode, 

though, what is generally a small amount of energy (but need not necessarily be) is 

allowed to leak through, contaminating the estimate to a greater or lesser degree, 

Such behavior ts equivalent to the sidelobe patterns developed In plane wave 


beamforming. 


Now consider a multiple beam least squares beamformer with M = 2. Let 
one of tlie two modes included in the beamformer be the mode actually present iin 
the signal: let the other be the same mode to which the single beam beamformer ts 
steered. The steering matrix for the multiple beam: algorithm: is then 


ЕЕ! ENEI 
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While algebraicly tedious, it is not difficult to directly evaluate equation (5.3) for 
the sensor cross-coherence matrix of equation (5.10) and the steering matrix given in 
equation (5.12). An estimate equivalent to that of equation (5.11) may be extracted 


from the lower right hand corner element of S, When this is done, it is found that 


во Ze 
а |“ = 0 (5.13) 


The relationship between the single beam and multiple beain least: squares 
beamformers is now evident. The single beam algorithm allows energy from one 
mode to leak into the estimates of other modes, while the multiple beam algorithm 
does not allow such leakage between modes that are included in the beamformer. 
The effect of the multiple beam approach is best described by saving that the 
effective beam pattern for each mode to be estimated is forced to have a null in the 
direction of each of the other modes also included in the beamformer. This 
elimination of cross-talk and sidelobe leakage between the various modal estimates 
represents the fundamental advantage of multiple beam algorithms over their single 


beam equivalents. 


5.3 Generalizations of the Least Squares Modal Beamformer and the 


MLM Modal Beamformer 


The results of Section 5.1 can be generalized through the introduction of a 
weighting matrix W in the error calculation of equation (5.2). W IS E К 
assumed to be Hermetian and positive definite. It controls the relative importance 
of the various terms that contribute to the error expression. The generalized form 


of equation (5.2) 15 
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=P'WP-P'WEÂ-ÂtEtWP+ÂtEŻWEÀ. 
for which the minimum total square error 15 generated by 

А = (ЕМЕ) ЕУР. (5.15) 
The residual error result equivalent to equation (5.7) 1s 

O nin = PY [W- WE(E*t WE) 'E* WIP. (5.16) 
and the modal cross-coherence matrix estimator becomes 

EN E SWR RS ER WwE) (5.17) 


The choice of W is somewhat arbitrary. The results reduce to those of the 
least squares beamformer if W is chosen to be an Y X .V identity matrix. Other 
choices are possible; the discussion here is limited to the selection often made 


high resolution beamforming 


p 
ZN 


W= (3.18) 


The motivation behind choosing the mverse of the sensor cross-coherence 
matrix to be the weighting matrix is quite simple, but the implementation of such a 
choice is rather complex. In theory. one desires to. weight. observations made on 
sensors having low noie more heavily than. those encountenng high none levels. 
Obviously, if the noise is uncorrelated between sensors and has the same power 


(variance) at all the hydrophones. then selecting W to be the identity matris ts 





reasonable. For spatially uncorrelated noise that is of varying intensity at different 
sensors, it is appropriate to choose W as a diagonal matrix with each non-zero 
element being proportional to the inverse of the noise power seen at the respective 
sensor. Such a selection most certainly weights the sensors receiving low noise levels 


more heavily than those receiving high noise levels. 


For spatially coherent noise, the concept is naturally extended by identifying 
the weighting matrix with the inverse of the full sensor cross-coherence matrix in a 
condition where no signal is present, so that only noise statistics are measured. 
Finally, it is shown in Appendix A that the results of equations (5.15) and (5.17) are 
invariant to the inclusion of signal components in the weighting matrix, so that the 
full sensor cross-coherence matrix may be used in place of just its noise portion. 
This invariance is very important, since it allows rn situ adaptive estimates of the 
weighting matrix to be made without requiring prior removal of the signal from the 
noise. To be valid, it does require the assumption of statistical independence 
between signal and noise, though, which then leads to the issue of coherent 


interference. This very significant problem is addressed at length in the sequel. 


The choice of W indicated in equation (5.13) allows equations (5.15) and (5.17) 
to be recast in more recognizable forms. After simplification, one obtains 


SE ES. (5.19) 
and 
S, = (Et S! E)! (5.20) 


These results are the modal analog of what is known as the MLM algorithm: the 
name vrrnmman energy method is also commonly used. equations (5.19) and (5.20) 


mav be derived in a number of different ways, each of which provides insight into 
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the properties of the solution. The name "MLM algorithm" comes from the fact 
that the solution is the maximum liklihood estimator for a known signal in zero 
mean, spatially coherent Gaussian norme. [t can also be shown that the solution 
provides the maximum SNR processor for a signal in statistically independent 


coherent noise. 


The MLM algorithm just derived represents the modal analog of the full 
multiple beam MLM beamformer. It agrees with the result of Hinich [42], who 
approached the problem from the perspective of a maximum liklihood estimator for 
a signal in Gaussian noise. The algorithm can also be implemented as an analog of 
the single beam least squares approach, where M is chosen as l and the resulting 
processor is repeatedly evaluated by indexing the steering vector E over the mode 
set of interest. The relationship between the single beam and multiple beam 
variants of the MLM algorithm is identical to that of the two least) squares 
approaches, in that the effective beam pattern for each mode included іш tlre 
multiple beam algorithm is forced to have nulls in the direction of all the other 
modes also included: such nulls are not generated in the single beam MLM 
algorithm. This null generation has some interesting consequences for MLN] 


algorithm: performance. These effects are discussed later mn the chapter. 


5.4 A Comparison of Modal Beamforming with Plane Wave 


Beamforming 


As ls apparent from the previous discussion, the most significant difference 
between a modal beamforming implementation and more traditional plane wave 
methods is in the construction of tlie steermg matrix E. If one considers the ~шоје 
beam case for simplicity, the steering veetor required for the detection of a plane 


wave traveling with wavenumber vector x Is 
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where F, is the three-dimensional position vector of the i sensor with respect to 
the receiver reference point. For the plane wave case, 4 has a magnitude equal to 
the wavenumber, and a direction that is parallel to the direction of propagation. It 
also is three-dimensional, so that it may include a component in x. Because a plane 
wave Is homogeneous In amplitude across any phase front, the terms of the steering 
vector all have unit magnitude. The phase of each term indicates the compensation 
necessary to correct for the phase accumulation of the wave in its travel from the 


receiver reference point to the sensor. 


г 


Bv contrast, the steering vector describing the m mode is 


MUT s 
L от) е " 


Where ¢ (=) denotes the shape of the mode, and £,, is the horizontal wavenumber 
vector associated with the mode. к, has a length equal to the mode's horizontal 
wavenumber and a direction parallel to the horizontal direction of propagation. In 
contrast to the plane wave case, both the wavenumber vector * and the position 
vectors 7; are now restricted to just the two horizontal dimensions. Also, the modal 
steering vector elements are non-uniformly weighted by the size of the mode shape 
at each sensors depth. These two differences are indicative of the fact that the 


transform used in the z direction is modal rather than Fourter. 


It is of interest to specialize these results to the case of a horizontal array. For 
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such a situation, the depths of all the sensors can be assumed equal. This makes all 
the 6, (2;) in the modal steering vector identical. For the plane wave steering 
vector, it also ensures that the phase contribution of the x component of the 
wavenumber vector is constant across all the array elements. For a horizontal 
array, therefore, each modal steering vector can be considered to be related to an 
equivalent plane wave steering vector. The only differences between these two 
steering vectors are a scalar phase factor and a scalar magnitude factor. One can 
easily compensate for the effect of the magnitude scalar, since it is equal to the size 
of the mode shape at the horizontal array deployment depth, while the phase scalar 
has no effect on the modal cross-coherence estimate, which provides only relative 
phase information. Thus. for a horizontal array, plane wave beamforming in the 
vertical direction provides information which is equivalent to the modal amplitude 
estimates of direct modal beamformers. “This observation becomes important in 
Chapter 7, where data from the horizontal array (which is processed with plane 
wave beamformers) is compared with direet modal beamforming results from the 


vertical array. 


5.5 Estimation of the Sensor Cross Coherence Matrix 


The beamformers used for the generation of the data presented in this thesis 
are fashioned around estimates of the modal cross-coliereuce matrix (equations (5.8} 
and (5.20)) ratlier than the estimates of tlie direct. complex. modal. amplitudes 
(equations (5.5) and (5.19)). There are several reasons for selecting this approach. 
As discussed in Section 5.1]. the eross-coherence matrix is related to the fundamental 
statistics of the random amplitude processes, Whereas direct — amplitude 
measurements provide data only about a single realization of these processes. Thus. 


the former can be considered. to. be. somewhat. more. stable than. the. latter. 
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Additionally, the cross-coherence matrix estimator is slightly more efficient 
computationally, because its implicit averaging reduces the number of data points 
that must be processed through the beamformer. Finally, the matrix estimate 
provides direct measurements of the coherence for the various modes in its off- 


diagonal terms. The only information lost in this approach is the absolute phase of 


the received signal, which is of no interest here. 


As written, equations (5.8) and (5.20) depend upon the sensor cross-coherence 
matrix S,,, which is an unknown statistic of the random propagation problem. In 
practice, S,. is first estimated over a finite window of data, and then this estimate Is 


substituted into the equations in place of the true sensor cross-cohereuce matrix. 


^ 
= 


Typically, Sọ would be computed by integrating over the band of frequencies 
containing the signal. For discrete time signals that have been complex 
demodulated, the calculation may be written as [10] 


B 
Sy, = y X (m) Х.(т) ; (5.23) 


т----В 


where ЗВ + I is the bandwidth of the signal and A (n) is the FFT of x n). tlie 


complex preprocessed time series received at the i sensor. 


A somewhat simpler procedure may be used in the present situation. Because 
of its narrowband nature, the signal may be assumed to reside m only a single 
frequency bin. The preprocessed time sertes from each sensor may then be viewed 
as successive estimates of that bin’s complex amplitude, since the digital filtering 
employed is very narrow and cuts off quite sharply. For these conditions. a natural 


estimate of the semsor cross-colrerence matrix is the direct tine average 
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This method has the advantage of not requiring FFT computations in the 
estimation procedure. The two approaches are easily shown to be equivalent. If the 
sum in equation (5.23) is extended to cover all frequency bins generated by the FFT, 
then the results of equations (5.23) and (5.24) must be equal bv Parseval's relation. 
Since the narrowband nature of the preprocessing ensures that any energy 
remaining in frequency bins outside of the 12.5 mHz analysis bandwidth is 


negligible, the resulting estimates must then be the same. 


The replacement of S by an estimate raises the issue of bias generation in 
the estimates of Sy This is not a problem with the least squares algoritlims, since 
only linear transformations are made on the unbiased estimate of the sensor cross- 
coherence matrix. Such is not the case for the MLM algorithms, though, because. 
they involve non-linear inversion operations. Both Capon and Goodman [11] and 
Scheer [73] have investigated the resulting bias in some detail. While this is 
generally a matter of some concern when using MLM > algorithms. the coherent 
_ interference problem (to be discussed presently) causes such large errors in this 
application that the bias issue may be considered moot here. | None the less, it 
should be noted that, in general, abs: ute levels returned bv MEM algorithms must 


be corrected in order to guarantee quantitative validity. 


5.6 Performance of the Various Modal Beamformers 


The previous discussion defines four different possible approaches to the modal 
beamforming problem. These approaches are: 


l. Single beam least squares (conventional) modal beamforming. 
2. Multiple beam least squares modal beamforming. 


3. single beam MEM modat beamforming. 
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4. Multiple beam MLM modal beamforming. 


Judicious selection from among these various methods requires that their 
relative performances be assessed under conditions likely to be encountered during 
the actual data reduction. The next several sections address this task. The length 
and complexity of the discussion make a short overview of the eventual conclusions 


ın order. 


The theoretical performance of the single beam least squares beamformer Is 
investigated in next section. This simple approach has two attractions: it is easy to 
implement; and it is quite robust. However. even fundamental theoretical 
calculations demonstrate that the FRAM IV vertical array allows only the first 
mode to be rehably estimated at the lowest) frequencies of interest. Such 
performance is clearly inadequate, so that the other algorithms must then be 


considered. 


The two MLM algorithms are investigated in the following section. Mode 
cohereuce is shown to be a major factor in the performance assessment of both these 
methods. The single beam algorithm is shown to be theoretically superior to the 
multiple beam algorithm if the individual acoustic modes that make up the signal 
are incoherent (phase random) with respect to each other. This method would be 
the processor of choice for such a situation, since it also offers higher resolution than 
either of the least squares approaches. One of the most siguificant conclusions of 
Chapter 7 is that tlie received mode field is highly coherent in tlie Arctic, though. 
Indeed, some preliminary indications of this conclusion have already been presented 
i] Figures 3-9 through 3-12. For è signal field consisting of coherent modes, both 
the single and multiple beam MLM algorithms perform Inadequatels when applied 
to the vertical array data. although they may still be successfully applied. to the 


horizontal array data. This issue is discussed in some length. and и ts demonstrated 
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that the fundamental cause of the difficulties is mutual coherent interference 


between the various modes. 


In the final section, the multiple beam least squares beamformer is studied. 
This approach is shown to possess somewhat better resolution than the single beam 
least squares processor and is not affected by the coherent interference issues that 
plague the two MLM algorithms. Therefore, the multiple beam least squares 


method is selected for use in reducing the vertical array data. 


5.7 Performance of the Single Beam Least Squares Beamformer 


The fundamental issue involved in any assessment of the performance of the 
single beam least squares Poner is its ability (or inability) to distinguish 
between signal energy propagating in different modes. To investigate this question, 
Figure 5-1 has been generated to provide a measure of the theoretical modal 
resolution that can be expected from the algorithm. This figure consists of a 
contour plot of a surface, the height of which represents the output of the 
beamformer for an array consisting of the 18 recorded elements of the FRAM IV 
vertical array. The vertical axis shows the mode which the beamformer ts to detect 
(or, in a more traditional parlance, the mode to which the beamformer is steered), 
While the horizontal axis indicates the one mode actually present m the artificially 
created signal field which the beamformer ts processing. The theoretical output of 
the beamformer for each possible combination of the mode actually present and the 
mode steered to is used to define the surface, which is then contoured on a dl? scale 
(2.5 dB per contour). The first 10 modes of the 30 Hz mode set lrave been used in 


the generation of Figure 5-1. 


In analogy to plane wave beamforming in bearing. each horizontal cross- 


section. of Figure 5-1 may be considered to be the modal "beam. pattern" that is 
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Figure 5-1: Theoretical Modal Sidelobe Pattern - 30 Hz Single Beam 
Least Squares Beamformer for 10 Modes (all 18 elements) 
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obtained when the array is steered to that particular mode. However, tt must be 
stressed that the mode domain differs from the bearing domain ђу being 
fuudamentally discrete in nature. Thus, for this and all other plots having one or 
more axes indexed in mode number, the results shown at non-integer values of mode 
number are only linear interpolations between the true data, which exists only for 
integer mode numbers. This technique of displaying results is used so that large 
quantities of information may be displaved in a compact manner. To emphasize the 
discrete nature of the data, a dark grid has been added to the plot at integer values 


of the mode number. 


The calculations displaved in Figure 5-1 are theoretical in that thev have been 
computed with perfectly known mode shapes, sensor positions. and sensor cross 
coherence matrices. In particular, complications caused by the substitution of an 
estimate for the actual Sy, matrix have not been considered. The cross-coherence 
matrix used in the computation of the figure consists of one mode of unit amplitude 
in background noise that 1s spatiallv white, 1.e 


Se EE EL. ЕИ 


where 


Е EE = (5.26) 


Here E. is the steering vector of the mode indicated by the horizontal axis. That is. 
it is the particular column of the steering matrix E defined in equation (1.4) that 
describes the mode assumed to be present in the signal field. à is the inverse of the 
effective modal signal-to-noise ratio: a value of 107 * has been used in the generation 
of Figure 5-1. 

X beamformer of perfect resolution would create a ridge of 0 dB outputs alone 


the major diagonal of the plot. This ridge would be of narrow cross-section, with 
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Figure 5-2: Theoretical Modal Sidelobe Pattern - 47 Hz Single Beam 
Least Squares Beamformer for 45 Modes (all 18 elements) 
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Figure 5-3: Theoretical Modal Sidelobe Pattern - 47 Hz Single Beam 
Least Squares Beamformer for 45 Modes (top 9 elements only) 
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Figure 5-3: Theoretical Modal Sidelobe Pattern - 17 Hz Single Beam 
Least Squares Beamformer for 45 Modes (top 9 elements only) 
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Figure 5-4: Theoretical Modal Sidelobe Pattern - 47 Hz Single Beam 
Least Squares Beamformer for 45 Modes (every other element only) 
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the output ımmediately fallıng to a floor of about -30 dB (the level of the white 
noise included in the cross-coherence matrix) for any off-diagonal points in the grid. 
Figure 5-1 shows quite graphically that the single beam least squares modal 
beamformer is far from ideal at 30 Hz when applied to the FRAM IV vertical array. 
The first mode ts clearly resolved; about 20 dB of separation 1s observed between it 
and the second mode. To lesser degrees, modes 2 and 3 can also be separated. with 
about -10 dB of leakage between modes 2 and 3. and about -7 dB between modes 3 
and 4. However, all modes beyond the third can be seen to be fundamentally 
unresolvable with the available array. For example, tle sidelobe between modes 4 
and 5 is, at best, down by only 3 dB; hence, energy from eitlier mode for mode 5 Is 
capable of generating large output levels when the beamformer is steered to mode 1. 
It can be concluded that the single beam least squares modal beamformer is 
theoretically capable of resolving no more than 3 modes at 30 Hz. The realities of 


dealing with actual data can be expected to degrade performance still further. - 


Similar investigations at other frequencies indicate that the algorithm can 
resolve only a single mode at 15 Hz, 5 modes at 47 Hz, and only about & modes even 
at 71 112. Figure 5-2 displays the results for the 18 elements of the FRAM IV 
vertical array at 47 ilz. Note that the range of modes displayed in Figure 5-2 
extends to 45 versus the 10 displayed in Figure 5-1. This is done to demonstrate the 
different forms of aliasing that are possible. The major modal sidelobes observed in 
Figure 5-2 can be divided into two types. One type Is a general broadening of the 
central ridge with increasing mode number. The other type is a general banding 
that occurs at right angles to the central ridge, and ìs prominent at higher mode 


numbers. 


The broadening of the central peak sidelobe structure in Figure 5-2 can be 


attributed to the problem of madequate array length. as shown m Figure 5-3. This 
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figure has been generated from the 47 Hz mode set in a fashion identical to that 
used in the creation of Figure 5-2, except that the bottom 9 elements in the array 
have been excluded from the beamformer, effectively halving the array aperture. 


The resultant large increase in the width of the central peak is readily observed. 


In contrast, Figure 5-4 has been recreated from Figure 5-2 by eliminating 
every other recorded sensor, rather than the bottom nine. The number of sensors 
has again been halved, but now in such a way as to retain the total array aperture. 
As expected, the width of the central peak in Figure 5-4 is comparable with that of 
Figure 5-2. However, the banding observed at the higher modes has now moved 
down in mode number. This ts clearly due to an inadequate sensor population in 
the available array aperture. Note that the pattern does not occur until one or both 
of the injected mode number and the beamformed mode number is greater than the 


number of sensors in the vertical array (13 in Figure 5-2 and 9 in Figure 5-4). 


To better understand the causes of the two types of allasing, one must first 
consider how the beamformer works in a physical sense. In essence. the single beam 
least squares algorithm uses the elements of the vertical array to effect a finite sum 
approximation of the mode orthogonality integral (equation (2.8)). The amount of 
aliasing is an inverse measure of the effectiveness of this sum in approximating the 
integral. Too great an mter-sensor spacing results in undersampling the aperture of 
the array, which, in turn, causes aliasing stmilar to that generated when a time 
series Is sampled below the Nyquist rate. This type of aliasing generally occurs 
between modes of significantly different mode number, as Is demonstrated by the 
banding effect it generates. On the other land, a vertical array aperture of 
inadequate length can cause significant portions of some of the mode shapes 
mvolved to be excluded from the integration. also causing anomalous results. 1n 


this case, though, itis generally modes that are close in mode number that look 
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most similar; hence, this type of aliasing primarily contributes by widening the 


ғ 


central ridge o? i 


3 


igures 5-1 through 5-4. A study of the mode shapes presented in 
Figures 4-4 through 4-7 verifies this conclusion. At every frequency, the first mode 
not resolved by the single beam least squares algorithm is also the first to have a 
significant portion of its mode shape extend below the 960 m length of the array. In 
some respects, the two types of aliasing are also akin to the local versus global error 


problems encountered in target detection and estimation theory [82]. 


From the above discussion, it becomes obvious that, in the frequency range of 
interest, the modal resolution available from the FRAM IV vertical array is 
fundamentally limited by its length. At the lowest frequencies available in the data 
set, not more than a single mode can’ be resolved by the conventional spatial 
matched-filter; even at the highest frequencies in the data set, only about 10° of 
the modes that might contribute significantly to the received. field can be resolved. 
Such performance is clearly inadequate, requiring consideration of higher resolution 


modal beamforming methods. 


5.8 Performance of the MLM Beamformers 


9.8.1 Performance of the Single Beam MLM Beamformer for Incoherent 


Modes 


When applied to plane wave beamforming, the MLM algorithm is normaly 
considered to be a high resolution alternative to least squares (conventional) 
beamforming. From 3 theoretical viewpoint, a similar interpretation can be made 
for modal beamforming. This is illustrated in Figure 5-5. The contour plot 
displaved in Uns figure is identical to that of Figure 5-1, except for the algorithm 
employed in computing the side lobe patterns. Here, it is the theoretical output of 


the single beam MLN algorithm at 30 Hz that ts shown, rather than the output of 
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Figure 5-5: Theoretical Modal Sidelobe Pattern - 30 Hz Single Beam 
MLM Beamformer for 28 Modes (all 18 elements) 
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Figure 5-5: Theoretical Modal Sidelobe Pattern - 50 H7 5inegle Deam 
MEM Beamformer for 28 Alodes (all IS elements} 





ET 
the least squares algorithm displayed earlier. The same SNR as used earlier (30 dB) 


has been assumed in the generation of this figure. 


The result displayed in Figure 5-5 is very close to that of an ideal modal 
beamformer, and provides the motivation for investigating the MLM algorithm in 
modal decomposition applications. However, practical issues involving the 
implementation of the processor typically limit this ideal theoretical performance. 
In particular, Figure 5-5, like Figures 5-1 through 5-4, has been computed using a 
perfectly known sensor cross-coherence matrix. “To measure the the effect of ising 
only an estimate of the I matrix, sets of synthetic data have been created. aud 
processed using the single beam MLM algorithm. The processed output for one of 
these data sets is displayed in Figure 5-6. This figure is once again a contour plot 
with the surface height representing modal amplitude on a dB scale. Here. though, 
the horizontal axis represents time, while the vertical axis represents the mode 
number to which the beamformer is steered. The plot consists of the time series 
over 20 minutes of the amplitude estimate for each of the 15 modes displaved. 
These estimates have been melded together to form the contoured surface by linear 
interpolation between adjacent modes. As was the case earlier, the only real data in 
the vertical occurs at integer values of the mode number. The contoured format ts 


used only because it allows information to be displaved ata high visual density. 


To interpret the results displayed in Figure 5-6, one must first understand the 
details involved in the creation of the data set. The sound field has been 
synthesized from the 30 Hz mode shapes presented in Chapter £, with each of the 
odd modes between the first and the fifteenth contributing to the sum with unit 
amplitude. Since only the odd modes are included in the sound field, the even mode 
amplitude ontputs indicate the amount of energy that the algorithm allows to leak 


into their estimates; To better differentiate between Various modes in the output, 
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Figure 5-6: Output of Single Beam MLM Beamformer versus Time 
for Incoherent Synthetic Data at 30 Hz 
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Figure 5-6: Output of Single Beam MEM Beamformer versus Time 
for Incoherent Synthetic Data at 30 Hz 





- 146- 
the unit amplitudes have been modulated by sinusoids with periods of 20 minutes: 
the full time extent of the plot corresponds to exactly one period of the modulation. 
The starting point of the modulating sinusoid has been adjusted by 00 for each 
succeeding odd mode, so that the modulation for mode 1 looks like a cosine 
function, while the modulation for mode 3 appears as a sine function, etc. Each of 
the various modal contributions has been provided with an independent random 
phase component, so that they are phase random with respect to each other. The 
resulting field can then be considered to be a sum of incoherent modes. The 
random phase fluctuations have been created with a temporal correlation length of 
32 sec (20 points), so as to approximate the effect of the final preprocessing filter. А 
small amount of spatially uncorrelated white noise (at about 30 dB SNR) is included 


in each hvdrophone's time series. 


The synthetic data set has been created in a format identical to that of tlie 
preprocessed real data; namely, a complex time series sampled every 1.6 seconds for 
each recorded sensor. The beamforming itself has been computed using an 
averaging window of 2 minutes for S. representing about 4 degrees of freedom. 
with updates once a minute (for a 50*6 overlap factor). The resulting estimate of 


the sensor cross-coherence matrix is subsequently stabilized for inversion by the 


addition of a small positive constant to the diagonal terms. 


In Figure 5-6, the first mode is clearly resolved. The amphtude modulation is 
apparent, demonstrating both the proper period and the correct initial phase. The 
output for the second mode indicates that the processing provides more than 20 dB 
of rejection between it and modes modes 1 and 3. Similarly, all modes through the 
seventh are well resolved, and modes 8 through 11 сап be at least partially 
distinguished. This performance certainly exceeds that of the single beam least 


squares algorithm, even without considering the effects of an using only an estimate 





BIDS 


of S. on the latter. 


5.8.2 The Relationship Between the Single and Multiple Beam Variants 
of the MLM Beamformer for Incoherent Modes 


Because Hinich [42] has proposed a multiple beam MLM algorithm for modal 
beamforming, it is important to consider the relationship of the two algorithms in 
the case of incoherently summed modes. To see the nature of the relationship. one 
must review briefly the technique by which the MLM algorithm assures maximum 
signal-to-noise ratios at the beamformer output. From a linear algebra perspective, 
the problem consists of determining V complex unknowns, each representing the 
complex gain applied to the observed field at a given sensor. The problem then 
contains N degrees of freedom. For the single beam case, the MLM result is 
obtained by effecting a constrained minimization of the noise power. The EE 
guarantees undistorted passage of the desired signal through the beamforming 
process, and uses one degree of freedom. The remaining N — 1 degrees are then 
employed to minimize the power of the noise in which the signal is embedded. This 
minimization may be visualized as an optimal placement of № = I nullis in the 


modal beam pattern. 


By comparison, the multiple beam variant can be considered to have M № 
degrees of freedom, representing the unknown complex gams to be applied at any 
given sensor for each of Af beams. Cursory consideration of the problem shows that 
at least Af distortionless constraints (one for each beam) are required. A more 
detailled investigation shows that M? constraints are actually placed upon the 
problem. The extra constraints represent decoupling requirements between the 
different beams. These are needed to guarantee accurate energy accounting: 


without them, energy detected on one beam could also leak throveh to other beams, 





SE 
and thus be accounted for multiple times. Thus, there are a total of M - (.N — M) 
degrees of freedom available for minimizing the noise on M beams, or only N — Af 
degrees of freedom for each beam. The multiple beam MLM algorithm can tlıen be 
viewed as the equivalent of single beam result with some extra decoupling 
constraints included. "These extra constraints have the effect of locating a priori 


some of the nulls that would otherwise be free for noise minimization. 


The upshot of this rather abstract discussion is that one always obtains better 
SNR performance from the single beam form of the MLM algorithm than from the 
multiple beam variant when processing an incoherently summed mode field. This is 
obvious, since the multiple beam approach requires prior placement of mulls that 
would otherwise be located to optimally reject noise; this prior null placement can 
never improve the total noise rejection. A full and mathematically rigorous 
argument to the same effect is presented in Appendix B. One simple indication of 
this behavior is that when M = N, both the least squares and MLM algorithms 


reduce to the same solution 


| 2 > 


_p-le (pt)! — 
MZE SN (E | , ОО] 
where the steering matrix E is now NY X N and may be assumed to be non-singular. 
Such a result is expected, since under these conditions all the available degrees of 


freedom are utilized as constraints, leaving none for the noise minimization that Is 


the source of the MLM algorithm's superior performance. 


The conclusion to be drawn from Figures 5-1 through 5-6 and the previous 
discussion is that the single beam MLM modal beamformer is the algorithni of 
choice for vertical array data modal decomposition when the sound field consists of 
incoherent modes. This finding must be tempered by the results of the following 
section, which show that it is very sensitive to how well the incoherent mode 


assumption is met ir practice. 
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Figure 5-7: Output of Single Beam MLM Beamformer versus Time 
for Incoherent Synthetic Data at 47 Hz 
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Figure 5-8: Output of Single Beam MLM Beamformer versus Time 
for Coherent Synthetic Data at 47 Hz 
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Figure 5-8: Output of Single Beam MLM Beamformer versus Time 
for Coherent Synthetic Data at 17 Hz 
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5.8.3 Performance of the Single Beam MLM Algorithm for Vertical 


Arrays in the Presence of Coherent Interference 


As has already been indicated in Chapter З апа is conclusively demonstrated 
in Chapter 7, the various complex mode amplitudes are highly coherent in the 
FRAM IV data set. This mode coherence has a significant effect on the 
performance both of the MLM algorithms, making neither the single beam processor 
nor the multiple beam processor applicable to the present problem. It is revealing 
to study this problem in some detail, so that the effects of mode coherence on the 


MLM algorithm can be better understood. 


Pigure 5-7 displays the output of the single bean MEM modal beamformer for 
a synthetic sound field at 47 Hz. In all other particulars except frequency, both the 
synthetic data generation and the processing is identical to that of Figure 5-6. Most 
importantly, the synthetic field once again consists of the odd modes between I and 
15 combined in a temporally phase random fashion, so that the various modal 
contributions can be considered to be mcoherent. Each mode is once more assigned 
an amplitude of unity and sinusoidally modulated. At 47 Hz, the first 9 modes are 
easily resolved, and all of the modes through mode 15 (the last one present tn the 


signal field) are at least partially resolved. 


Like Figure 5-7, Figure 5-8 also displays the output of the single beam MEM 
modal beamformer for a synthetic sound field at 47 Hz. Both the syntlietic data 
generation and the processing is identical to that of both Figure 5-6 and Figure 5-7, 
with one important difference; while the modal contributions to the sound field in 
Figures 5-6 and 5-7 are summed incoherently by injecting random phase 
fluctuations that vary mdependently from mode to mode. the modal contributions 


for Figure 5-8 are summed in a coherent fashion. Vhis ts accomplished by injecting 


3 phase fluctirations that are constant across all niodes. so that а constant phase 





-152- 
relationship is maintained over time between the various modal amplitudes. In 
striking contrast to the results of Figure 5-7, the resulting output from the single 


beam MLM beamformer for this case is so poor as to be essentially meaningless. 


Having demonstrated numerically that the issue of mode coherence is of 
extremely important in determining the performance of MLM modal о а 
it is next important to develop a more intuitive understanding of the problem. 
Consider a sound field consisting of a single mode ın spatially white noise. Let the 
the complex amplitude of the mode that is present be а, and the steering vector 
associated. with it be denoted by E, This modal steering vector 1s the particular 
column of the steering matrix E defined in equation (1.4) that describes tlie mode 
assumed to be present in the signal field. The sensor cross-coherence matrix for 


such a situation 15 


9 9 oe 
Є = |а |" Е E? +1]la |“ [E |-1 2.28 
where X is the inverse of the effective modal signal-to-noise ratio and 
9 
Е |" = ЕТЕ . (5.29) 
Ep e 


The single beam MLM beamformer estimates the squared magnitude of the 
; : A. 2 1? 
complex amplitude of the mode to which it is steered. [a |”, as 


1 


= (5.30) 
(ES E) 


. 19 
ШЕ = 
In contrast to Е. Which is the steering vector of the mode actually present m the 
signal field, E, is the steering vector of the mode to which the beamformer is 
steered. Using the identity given in equation (1.3), 1t Is possible t0 compute the 
Inverse of Se as given by equation (5.28). This result may then be substituted into 


equation (5.30) to vield 
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Here p, is the normalized inner product between E, and E 
bp =? =p 


ЕЕ 
жс ee (5.32) 


Pp = . 
^ |EJIEJ 
The magnitude of Php always lies in a range from 0.0 to 1.0. 


Figure 5-5 is nothing more than equation (5.31) evaluated at various 
combinations of E, and E, It is easy to see how the MLM beamformer obtains its 
high resolution. When the beamformer is steered to the mode that is present in the 


field, Poy takes on a value of 1.0, leading to an output of 


9 


а | = ја |“ Це ЈЕ (5.33) 


р 
On the other hand, when the mode present in the signal field is different from the 
mode to which the beamformer is steered, the value of 77 is small if the two modes 
are well resolved by the array. In fact, for the limiting case of a vertical array of 
continuous aperture extending across the full sound channel, 2» Is guaranteed to be 
zero for all possible mode pairs, since the numerator of equation (5.32) approaches 
the orthogonality integral of equation (2.8). Thus, the squared magnitude of 7» IS 
generally negligible for well resolved modes, leading to 

ja ЈЕ |“ 


ај == | = ] . (5.31) 
IE, 


2 
› 


Since the value of A is small for high SNR signals, (in Figure 5-5, a value of 10” 1s 
used), excellent rejection is obtained when the mode present is not the ome to which 


the beamformer Is steered. 


In the case where the sound field consists of two modes whose contributions 
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are incoherent with respect to each other, the sensor cross-coherence matrix can be 
written as 
5 


== 2 + 2 + 2 27 и 
Sy — м Е T EL E + |а [E^ I. (3.35) 


pl 


where E, and Бо are the steering vectors of the two modes contributing to the 
signal field. Here, for simplicity, the white noise scaling has been kept identical to 
that of equation (5.28). Equation (5.35) can also be inverted and then substituted 
into (5.30) to obtain a general result. In this case, it is of the most mterest to study 
the effect of the presence of the second mode on the beamformer estimate for the 


first mode, which can be obtained by setting E, equal to Ба This leads to 


| : A (05 +1 
po = lanl E | (5.36} 
where 
e) e) 
Eos 
b= — M — (Бал 
a PIE, 


Again, И 15 easy to see how high resolution is obtained; if the two modes 


present are well resolved by the array, then the squared magnitude of p,, , is once 


plp 
more negligible, giving an output of 


Ж 2 жесе 2 И N 
|а E Lm (1 + МЕ (9.38) 
This result is identical to that of equation (5.33), which means that the single beam 
MLM beamformer almost completely elimmates the effect of the presence of the 
second incoherent mode in the signal field. Such a conclusion is certainty in 


agreement with the results of Figures 5-6 and 5-7. 


'Fhe ability. of the single beam. MEM algorithm t0. discriminate against. the 


presence of the second mode is inherently tied to the ability of the array to resolve 
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vw) р 
the two modes. If the two modes are well resolved, then ЈУНИ is small. and 
therefore justifiably negligible. On the other hand, if the two modes are 
fundamentally unresolvable, then ignoring this term is a poor approximation at best. 


In either case, though, the performance is at least as good as that of the single beam 


least squares beamformer. 


To highlight the difference between case where the *he signal modes combine 
incoherently and the case where they remain phase locked with respect to each 
other, the preceding problem is now repeated for a sound field consisting of two 


modes that sum coherently. In this situation, the sensor cross-coherence matrix is 


- \+ - 
Sy = (а E ға, Е.) (а Ба p E ој c (251) 


э 22027 
\ |а, ЕТ. 


Again, the white noise scaling has been kept identical to that of equation (5.28). 
Comparing equations (5.35) and (5.39), it can be seen that the effect of the modes 
being coherent ts to force the retention of cross-terms that would average to zero in 
the incoherent case. Inversion of equation (5.39) and subsequent substitution into 
(5.30) is easily accomplished. Upon steering the beamformer to the first mode, one 


obtains 


ја [a (bx +c 
JA (5.10) 


. SR Ж 
um (b) +1) — |р 


01021 
where b is again defined as in equation (5.37) and c is 


la Е ae [s 


p2 =p? 


5 5 (5.41) 
а о! Е ol 


= 


5 *) 
Assuming as before that the modes are well resolved by the array, а (ШІ 


ve 


again be neglected, so that for typical white noise levels the output is approximately 
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The resulting output is now proportional to X. If both the modal steering vectors 
are assumed to be of roughly equal norm, and if the modal amplitudes are 
approximately equal, then c can take on values of between 0.0 and 4.0, depending 
on the relative phase between the two modes and how well they are resolved by the 
array. For well resolved modes, c approaches a value of 2.0, since the steering 
vectors associated with the modes are close to orthogonal. On the other land. 
equation (5.34) demonstrates that the background noise floor that leaks mto any 
mode is approximately A 47 under the assumption of equal norms for different 
modal steering vectors. Thus, the presence in the signal field of two coherent, well 
resolved modes of roughly equal strength causes the single beam MLM algorithm to 
generate output estimates of these modes that are only 3 dB above the noise floor, 


no matter how strong the actual modal signals. This ts exactly the phenomenon 


seen in Figure 5-8. 


Further investigation of equation (5.40) vields two significant conclusions. 
First, even very small levels of a second coherent mode can trigger tlie interference 
effect. Again assuming that the modes are well resolved, consider tlie limit. as ШЕ 
vanishes. For this case, c approaches the value of 6, and equation (5.10) then 


siniplifies to 


|a ——. (5.43) 


: D А : 
The result approaches the correct value of Ja 1 only if 5X is nnıch greater than 1.0, 
indicating that coherent contributions from other modes are negligible to the MEM 
processor only if they are at or below the background notse level. Second, the effect 


of colrerent interference oi the single beam MILAM algorithm: is worse when the two 





modal steering vectors are well resolved than it is when they are poorly resolved. lu 


u. ә : 
the limit as Е о approaches Е р |p |" takes on a value of unity. and the output 


plp2 


becomes a very reasonable 


+a 


i p (5.44) 
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This is exactly the opposite of what might be intuitively expected. Beamforming 
algorithms typically experience greater difficulty when multiple sources look very 
similar than when they are easily distinguished. But the nature of equation (5.10) 1s 
such that more realistic results are achieved from the single beam MIM algorithm 
when the value of [ОО is near 1.0 (indicating that the interfering modes are 
poorly resolved) than when it is close to 0.0 (meaning that the interfering modes are 


well resolved ). 


For the ideal case of a continuous vertical aperture extending across the full 
channel, the single beam MLM approach is completely inappropriate when the 
modes contributing to the signal field remain coherent. This Is because mode 
orthogonality guarantees that рр? Is always zero, a situation which leads to the 
worst possible performance of the MILM processor. As will be investigated m the 


sequel, its suitability for use with horizontal arrays is somewhat better. 


The fundamental reason that coherently interfering modes cause the single 
beam MLM algorithm difficulty is that they violate the assumption of statistical 
independence between signal and noise. In the single beam approach, the signal, by 
definition, is the mode selected for detection; any other modes present аге 
considered to be part of the noise against which the processor Is to discriminate. If 
the modes are phase random with respect to each other, then the selection of one 
mode as tlie signal while tlie others are included in the noise causes no violation of 


this fundamental assumption. On the other hand. the definition of a single mode as 
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the signal destroys the assumption of independence between signal and noise when 
the modes remain coherent. The effect of the resulting coherent mterference Is to 
corrupt the spatial structure of the signal mode. The beamformer then rejects the 
energy carried in the signal mode because the mode is no longer recognizable. This 
is why equation (5.42) indicates that the output is roughly the same as background 


noise levels. 


5.8.4 Performance of the Multiple Beam MLM Algorithm for Vertical 


Arrays in the Presence of Coherent Interference 


Since the single beam MUM algorithm is clearly inappropriate for the task at 
hand, the effect of coherent modal interference on the multiple beam variant is now 
investigated. From the discussion in Sections 5.2 and 5.3, it can be seen that the 
fundamental advantage of the multiple beam approach is that it expands the 
number of modes which are simultaneously designated as signal from 1 to M, the 
number of beams included in the multiple beam beamformer. If a multiple beam 
MLM algorithm that includes all modes making significant coherent contributions to 
the total signal field rs implemented, then the statistical independence of signal and 
noise can be reestablished, and the resulting algorithm should, in theory, produce 
acceptable results. .À number of practical considerations limit tlie applicability. of 


this approach, though. 


The issue of greatest importance is the proper selection of Af, the number of 
modes to include m the beamformer. As mentioned earher, this choirce is basically a 
modeling decision, since the inclusion of extra modes in the beamformer is 
tantamount to improving the initially assumed signal model. — Including. an 
insufficient number of modes leaves one with the same problem that plagues the 


single beam approach; the algorithm: performs poorly because the signal and. tlie 
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noise, as implicitly defined by the model selection, are not statistically independent. 
On the other hand, making Af too large is also detrimental, for two reasons. The 
original derivation limits the number of modes that can be handled simultaneously 
to the number of sensors available. Otherwise, Et S E is of less than full rank, 


and therefore uninvertible. 


For the FRAM IV vertical array (and many other vertical arrays), tlie fact 
that many of the higher order modes cannot be adequately resolved places an even 
more restrictive limit on M. Since two modes become unresolvable when tlie 
steering vectors describing them become too similar, the inclusion of multiple 
unresolvable modes in the steering matrix quickly drives it towards rank deficiency. 
so that inversion of the quadratic product creates serious difficulties from the 
viewpoint of numerical stability. The impact of this issue on the multiple beam 
least squares algorithm is studied in some detail in Section 5.9. Its effect on the 
multiple beam MLM algorithm is similar, hmiting the effective number of beams 
that can be handled simultaneously to as few as 3 or 4 below 20 Hz, and no more 


than perhaps a dozen even at 71 Hz. 


It has been shown previously that any mode contributing coherently to the 
signal field in an amount that exceeds the background noise level must De 
considered significant (see equation (5.43)). Thus, one reaches the conclusion that it 
is not possible to implement a multiple beam MLM processor capable of dealing 
with the coherent interference problem for the FRAM IV vertical array. However, 
given a sufficiently improved array, such an approach could be possible. The 
primary remaining issue under these conditions would be one’s ability to invert 
large dimension arrays (perhaps 30 x 30 or 50 X 50) im a numerically stable 
fashion. The array necessary to support such processing would need to extend 


much further in depth. so that ıt could adequately resolve significantly greater 
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numbers of modes. Only a relatively modest increase in the number of hydrophones 
would be needed. The addition of further sensors beyond this amount would be 
beneficial, though, since the performance improvement that MLM algorithms over 
equivalent least squares algorithms is fundamentally linked to the size of N — M 


(the excess of sensor count over beam count). 


5.8.5 Alternative Implementations of the MLM Algorithm in the 


Presence of Coherent Interference 


Alternative implementations of the MLM = algorithm in the presence of 
coherent interference are possible if one is not interested in actually measuring 
individual modal characteristics. While these approaches are generally not 
applicable to the present problem, it is of interest to understand their relationship 
to it. 

Consider again the sensor cross-coherence matrix for two coherently 
interfering modes (equation (5.39)). A consideration of the eigenvalues and 
eigenvectors of this matrix quickly leads to the conclusion that the output of the 
single beam MLM beamformer is maximized when a steering vector of the form 


Е, = К (a, E, +a Е.) (5.15) 


p-— 
Is used [20j.. Here & is an arbitrary scaling factor that would typically be used for 
normalization of the steering vector. This observation forms the basis for all these 
alternative implementations. If the relative mode amplitudes and phases were 
known «a priori, then an optimal detector could be created by beamforming to a 
steering vector which was the properly weighted sum of the individual modal 
steering vectors. In this thesis, of course, these are exactly the parameters to be 
measured, but there may be ways of predicting them accurately enough to allow 


significant improvements im detection processing. 
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One alternative approach ıs to index the beamforming across a set of steering 
vectors generated by various weighted sums of the modal steering vectors; arbitrary 
variations of the modal amplitudes and phases become impractical beyond a very 
few modes, however. A more promising variant Is to index the beamforming across 
the two parameters of source range and depth, using one or another of the different 
propagation theories currently available to predict the mode amplitudes and phases 
needed to compute the sum. The advantage of such à method is that source range 
and depth are immediately available once detection is made. Several preliminary 
efforts in this direction show promise. The work of Fizell and Wales [34] is of 


special note, since It applies the approach to actual field data. 


5.8.6 Performance of MLM Algorithms for Horizontal Arrays in the 


Presence of Coherent Interference 


After much effort, two Important conclusions have been reached concerning 
the use of the MLM algorithm for reduction of the FRAM IV vertical array data set. 
The first is that the MLM algorithm works very well when the modes making up 
the sound field are incoherent. For this situation. the single beam variant is 
guaranteed to provide better performance than any multiple beam variant. The 
second is that the MLM algorithm has a great deal of difficulty when the modes 
comprising the sound field remain coherent with respect to each other. In this case, 
the performance is worst when the coherently interfering modes are well resolved by 
the array rather than poorly resolved, as might otherwise be expected. The 
difficulties that arise when coherent interference Is encountered are severe enough to 
make both the single and multiple beam versions unattractive for use with the 
vertical array data. On the other hand, the MLM approach can still be used on the 
horizontal array data if adequate care and caution are exercised. И is revealing to 


Investigate why better performance might be expected. of the single beam MLM 
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Е: 


Figure 5-9: Components of a Cohereut Sound Field for a Notional 
Horizontal Array 


beamformer mn a coherent mode field when it is applied to horizontal array data 


than when it is used with vertical array data. 


Consider a situation where a horizontal array Is receiviug a large number of 
coherent modal arrivals. Let the steering vector associated with the KU node be 


Е 


= The sensor cross-coherence matrix 


„and its complex amplitude be denoted by a. 
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is then of the form 


Sy = E, E} + o” I, (5.46) 


= 


where E, is a weighted sum the steering vectors associated with the various modes 
that are present 

ER Y ЕК (5.17) 
This sum ıs shou Bietörially in Figure 5-9. As discussed in Section 5.4, the modal 
steering vectors for a horizontal array are identical to the steering vectors for plane 
waves, with the exception of a scalar equal to the size of the mode shape at the 
array's deployment depth. The equivalent plane waves are all traveling in tlie same 
horizontal direction but with slightly different horizontal phase speeds. Under these 
conditions, the modal steering vectors can be expected to be almost parallel to both 
one another and to the total sum. Figure 5-9 has been drawn to reflect this 
situation. 


h mode. Its 


Assume that the single beam MLM beamformer is steered to the ;' 
output can then be evaluated directly from equation (5.31) by setting m 1. 
E, — E. and E, =E. Since E and E, are almost parallel, it is not unreasonable 


to represent their inner product as 


-= ļ] а (5. 481 





where e,, 15 small with respect to unity but still large when compared to à. The 
case where є. 15 the same or smaller than ^ is not of interest, since this generally 
corresponds to operational situations where the array is fundamentally unable to 
resolve the modal structure (such as à horizontal hne array operated near 


broadside). The amplitude estimate that the beamformer makes for the II anode is 


then 
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a = | = | —. (5.19) 
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Е 
St 15 
Although this output is still proportional to the noise level ^, the effect is greatly 


mitigated by the small denominator ¢.., which is typically on the order of 107! or 


1077. 


The result presented in equation (5.49) raises the very important question of 
why the resultant output should bear any resemblence to the actual modal 
amplitude distribution. It 1s not at all obvious why the estimate should be large 
when the actual mode amplitude 15 large and why it should be small when the mode 
amplitude is small. Of the four terms on the right hand side of equation (5.49). only 
the two in the denominator depend upon the mode to which the beamformer is 
steered. The variation of IE | with mode number is only a scalar related to the size 
of the mode's shape at the array depth, and, hence, of no practical consequence. 
Indeed, if the steering vector of the equivalent plane wave is instead considered, 
even this simple variation with mode number is eliminated, and the term can be 
considered to be constant. Thus, the change in output level with mode number is 
principally caused by the e, term, so that the estimator output tracks the actual 
modal amplitude distribution only if e. is generally small for modes of large 
amplitude and big for modes of small amplitude. While this is never guaranteed to 
be the case, a graphical argument can be made to support the conclusion that it ts 


often so. 


Consider once more the sum of equation (5.47) and its graphical depiction 
shown in Figure 5-9. E can be expected to lie in a direction that is closest to the 
direction of the modes providing the largest contributions. For these modes, then, 


c, ls generally smaller than average, providing the desired peak mm the modal 
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amplitude distribution. The effect is most reliable when a single modal contribution 
dominates all others, as is often the case in the Arctic. This situation is shown in 
Figure 5-9. In creating this figure, the lengths of the various modal steering vectors 
have been assumed to be roughly equal, and the amplitude of mode 1 has been 
assumed to be much greater than the amplitude of any other mode. For this 
situation, it is apparent that the angle between E, and E, is smaller than the angle 


E, makes with other modal steering vectors. 


It is the natural sensitivity of equation (5.49) to variations in e, that provides 
the MLM algorithm with better resolution than the conventional beamformer For 
comparison, the output of the single beam least) squares beamformer is 


approximately 


9 





талы” IE, == 
|a dre = ЕР (1 = 6.) . (5.50) 
= 


While this result provides a much better absolute level than equation (5.49), its form 
is such as to make it significantly less sensitive to variations in є. than the MLM 


result. 


Thus, the single beam MLM beamformer ends up providing the desired result 
of higher resolution than the conventional beamformer when. applied. to. modal 
separation with horizontal arrays, although not for the reasons typically given to 


explain its operation. Additionally, the linkage through е. is highly non-linear, 


Ls 
making both absolute levels and relative peak levels highly suspect, although 
positions of the peaks in mode number (or, equivalently, grazing angle or phase 
speed) can generally be accepted as accurate. In previous cases where the MEM 
algorithm has been used in the presence of multiple coherent sources. or with a 


single coherent) source in a stable multi-path environment, these effects have 


probably been incorrectly attributed to the effect of the MEM bias discussed. earlier. 
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The distinction is important, because in many situations one can accurately 
compensate for the bias effect, whereas the extraction of reliable levels in the 


presence of coherent interference is Impossible. 


In conclusion, the MLM algorithm still provides higher resolution in the 
presence of coherent interference than conventional processing does when applied to 
horizontal array data, despite the fact that the various arrivals are highly coherent. 
This higher resolution is achieved at the expense of accurate peak levels and а 
marked decrease in the robustness of the processing, though. The results are most 
trustworthy in the case of one strong modal arrival dominating all other coherent 
contributions to the signal field. All of these effects can be seen in the resultant 


outputs of the MLM processor that are presented in Chapter 7. 


5.9 Performance of the Multiple Beam Least Squares Beamformer 


Since neither form of the MLM algorithm can be reliably. applied. to. the 
FRAM IV vertical array data in a modal beamforming context, the multiple beam 
least squares beamformer must then be considered as a high resolution alternative, 
Such an approach offers several advantages. First, the multiple beam algorithm can 
be expected to provide at least marginally better performance than the single beam 
variant if the number of modes included in the beamformer is not excessive. 
Because the multiple beam algorithm essentially attempts to fit the observed 
pressure field with a model containing a greater number of parameters, closer 
agreement with the data generally results. Additionally, since the derivation of the 
multiple beam least squares algorithm requires no assumption about the statistical 
relationship between the signal and the noise, its performance Is not degraded by 
coherent interference. It also provides some relief from the inter-modal aliasing 


problems seen in the single beam least squares algorithm, although the improvement 
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Figure 5-10: Output of Multiple Beam Least Squares Beamformer 
versus Time for Coherent Synthetic Data at 47 Hz 
(7 modes included in beamformer) 
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Figure 5-11: Output of Single Beam Least Squares Beamformer 
versus Time for Coherent Synthetic Data at 47 Hz 
(7 modes shown) 
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is not as great as that which the single beam MLM algorithm would provide if it 
could be applied. Finally, the multiple beam version i$ somewhat more 


computationally efficient than the single beam variant. 


A typical set of modal amplitude estimates generated by the multiple beam 
least squares approach is presented in Figure 5-10. For comparison, the equivalent 
set of outputs from the single beam variant is presented in Figure 5-11. Similar to 
Figures 5-6 through 5-8, the modal amplitudes are depicted as contours over a 
surface consisting of time on the horizontal axis and mode number on the vertical 
axis. As mentioned earlier, the only real data points occur for integer values of the 
mode number; results plotted at non-integer mode values are merely interpolations 


between adjacent data points. 


The particulars of the synthetic data sets used to generate both of these | 
figures are identical to those of the earlier figures. Again, the synthetic signal field 
consists of a coherent sum of the odd modes between | апа 15, each with unit 
amplitude. The various modal amplitudes have once more been alternatively 
modulated with cosine and sine waveforms with periods of 20 minutes. Both plots 
span the first seven modes at 47 Hz. All seven are resolvable in Figure 5-10. About 
20 dB of rejection 1s observed in the second mode and about 17 dB is seen in thre 
fourth mode. These results show some Improvement over those of Figure 5-11, 
where only the lower order modes can effectively be resolved, and where rejection in 


the second and fourth mode troughs is somewhat less. 


While the additional resolution provided by the multiple beam least squares 
algorithm Is not great, it is significant. This ts especially true at frequencies below 
20 Hz. where the single beam approach has difficulty resolving more than one mode. 
It is for this reason that the multiple beam least squares beamformer has been 


selected as the most attractive alternative for the modal decomposition of the 
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FRAM IV vertical array data set. 


5.10 Determination of the Number of Modes to Model 


The final topic to be addressed in this chapter is the proper choice of M, the 
number of modes to be included in the multiple beam algorithm. Ideally, the choice 
would include any mode likely to provide a significant contribution to the total 
sound field seen at the receiver. However, the inclusion of an excessive number of 
modes in the beamformer can lead to problems involving numerical stability, as can 
be seen in Figure 5-12. In this figure, the same 47 Hz synthetic data used to 
generate Figures 5-10 and 5-11 has been reprocessed with a multiple beam least 
squares algorithm that includes the first erght modes instead of the seven used 
above. By increasing Af from seven to eight, the effective resolution has been 
decreased rather than increased; only four modes are now actually resolved. while 


the rest are saturated by processing-induced noise. 


To understand this effect, consider the output of the multiple beam least 
squares processor for an input consisting only of spatially white noise. For such a 


case, the sensor cross-coherence matrix Is 


€) 


Sy 2 (2.21) 
From equation (5.8), the output of the processor can then be evaluated as 

а DEEST er 

Su Ee (ЕБ Е) (3.92) 


It is obvious that as the M Xx M matrix E* E. becomes. singular, one or more 
outputs of the processor grows without bound. Numerical. nistability tm the 
Inversion can then make the algorithm extremely sensitive to even small amounts of 


background noise. Therefore, a proper choice of M must be based on a compromise 
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Figure 5-12: Output of Multiple Beam Least Squares Beamformer 
versus Time for Coherent Synthetic Data at 47 Hz 
(8 modes included in beamformer) 
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between the increased resolution and the decreased numerical stability caused by 


the inclusion of additional modes. 


The relative singularity of the above matrix can be quantified through the use 


of the singularity coefficient g, which is defined as 


Det [Et E] i 
а == Биш с. 1239 
E IE I” 


pz] 


еге, (фе ре | | notation Is used to indicate the deterininant ol the Рм у а је 





E, through E, are the various column vectors comprising the steering matrix E. 


The denominator consists of the product of the norms of these vectors. 


The singularity coefficient has several very useful properties. First. it is 
always unity when Af = 1. Second, for any other value of M, it is unity only if all 
the component column vectors are all mutually orthogonal. Finally, its value is 
always in the range from 0.0 to 1.0. Since the coefficient is proportional to the 
value of the matrix determinant, it can approach zero only as the matrix becomes 
singular. Thus, it provides a good indication of the relative singularity of the 
matrix Et Е. 

Figure 5-13 provides a dB scale plot of this singularity coefficient for steering 
matrices consisting of up to the first nme modes at frequencies of 15 Hz, 30 112. and 
47 Hz. In general, the coefficient remains close to 1.0 when small numbers of modes 
are included in the steering matrix, gradually decreasing as the mode count 
increases. Beyond a certain critical number of modes, though, inclusion of further 
modes in the steering matrix then causes the coefficient to fall rapidly towards zero. 
Various methods may be used to define precisely this critical number of modes. 
Here, a trial and error procedure has been used to find a reasonable breakpoint. 
The singularity coefficient has been compared against the processing output for 
synthetic data similar to that found in Figures 5-10 and 5-12 at various 
combinations of frequency and mode count. This comparison shows that optimum 
beamformer performance occurs for values of the singularity coefficient that are 
slightly greater than 0.4. The critical mode count is then the largest number of 
modes that can be included in the steering matrix without the coefficient exceeding 


this threshold. 


This definition produces optimal results at all the frequencies tested in the 


MO and error procedure. llowever, it may ђе influenced by the arras 
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configuration, the modal structure of the particular channel investigated. and the 
choice of modes included in the synthetic data. The sensitivity of the definition to 


these various effects has not been investigated. 


Using the above definition, the number of modes that can be included in the 
multiple beam algorithm at any frequency is readily found. This information is 
presented in Table 5-I for the frequencies available in the FRAM IV vertical array 
data set. Since the multiple beam processor computes the output of all its beams 
simultaneously, the value also represents the number of modes that can be resolved 
at each frequency. In all cases, at least three modes are resolvable; the number 
increases with frequency, reaching a value of 9 modes at 71 Mz. A comparison with 
the mode shapes of ee 4 reveals that the last resolvable mode 15 the one in 


which the lowest oscillation of the mode shape is still sampled by at least one sensor. 


Frequency . Mode Count 
15.00 Hz 3 
17.75 Hz 3 
20.00 Hz 3 
23. 50 17 4 
30.00 Hz 5 
33 795917 5 


35.50 Hz 5 
47.00 Hz T 
OZ 7 
55.00 Hz 7 
71.00 Hz 9 


Table 5-I: \odes Resolvable by Multiple Beam Least Squares 
Modal Beamformer 
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5.11 Summary 


Several different methods for the decomposition into its modal components of 
the sound field observed at an array have been presented and analyzed. Modal 
decomposition is intimately related to the problem of beamforming, which has a rich 
history in the literature. The primary theoretical difference between the two results 
is the choice of the basis set of functions on which the received field ts decomposed: 
traditional beamforming utilizes the set of plane waves, while the present effort 
employs the normal mode functions associated with the sound channel in which the 
receiving array resides. When the two approaches are compared, the major 
advantages of plane wave beamforming are the independence of the technique from 
knowledge of the sound channel and its equivalence with Fourier. transform 
techniques, which allows a wide body of knowledge to also be tapped. The principal 
advantages of modal beamforming include the mathematically efficient nature of 


the representation and the physical relevance of the resultant outputs. 


Two different approaches to modal beamforming have been addressed: the 
least squares method and the MLM algorithm. The two are differentiated by their 
choice of weighting matrix used im computing the total square error to be 
minimized. The least squares approach werghts all errors equally, attempting to 
match the full received field. The MLM algorithm, on the other hand, tries to 
compensate for norse by selecting the inverse of the sensor cross-coherence matrix as 
the weighting matrix. ‘Tins selection has the effect of weighting noisy sensors less 
heavily than sensors with clean signals for purposes of calculating the total square 
error. Both methods can be implemented іп either à. single or. multiple beam 
format, generating four different approaches, all of which have been considered 


here. 
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As might be expected, the resulting solutions for these different approaches to 
modal beamforming all take on the forms identical to their plane wave processing 
algorithm analogs. The only difference occurs in the definition of the steering 
matrix, which is the parameter describing the signal to be detected. In plane wave 
beamforming, this steering matrix (or steering vector, when dealing with single 
beam variants) consists solely of the complex phase accumulated ш the travel of the 
plane wave from the receiver reference point to the individual sensor. For modal 
beamforming, the direction of travel is restricted to the horizontal plane. but the 
size of the mode shape at each sensor depth is included in addition to the phase 


accumulation due to travel. 


The relative performance of the each of the four variants proposed for modal 
beamforming has been investigated for the case of the FRAM IV vertical array 
when deployed in the sound channel described in Chapter 4. The single beam least 
squares beamformer, which is the equivalent of the conventional beamformer, has 
been shown to exhibit inadequate modal resolution, being able to resolve only a 
single mode at 15 Hz, and no more than 8 modes even at 71 Hz. At intermediate 
frequencies, typical resolution is 3 modes at 30 Hz and 5 modes at 47 lz. In all 
cases, the: number of modes that can be resolved is but a small percentage of those 
that might possibly provide significant contributions to the received sound field. 
The resolution has been shown to be limited by the length of the array rather than 


by sensor spacing within the array aperture. 


The single beam MLM beamformer can provide significantly better resolution 
than the single beam least squares beamformer if the received sound field is 
comprised of modes which sum incoherently. И has also been demonstrated 
theoretically that the single beam variant of the MLM algorithm always ехјиђих 


better performance than the multiple beam varrant under these conditions. Thus. 





the single beam MLM algorithm is the method of choice for accomplishing modal 
decomposition with the FRAM IV vertical array if the modes at the receiver 


combine incoherently. 


The performance of the MLM algorithm has been shown to be very sensitive 
to the assumption of mode incoherence. The single beam variant of the algorithm 
operates only very poorly in the presence of coherently summed modes. This 
problem has been examined in detail, both theoretically and through simulation: the 
cause has been shown to be a violation of the assumed independence of signal and 
noise. This violation is produced because one of the coherent modes ts implicitly 
classified as the signal and the others as noise in the single beam MLM variant. lt 
has been shown that this problem can be triggered by the presence of even small” 
amounts of a second coherent mode, and is most severe for the case of orthogonal 
steering vectors, as is the case for modes well resolved by a vertical array. The 
effect has been demonstrated to be less significant when the steering vectors of the 
coherently summing modes are approximately parallel to each other, such as occurs 
when considering high resolution horizontal array beamforming. Finally, it can be 
demonstrated that, in the presence of coherent modes, the single beam MILAM 
beamformer still provides higher resolution than the single beam least) squares 
beamformer when apphed to the horizontal array. In this case, the peak locations 


normally remain accurate; the resultant output levels are not reliable, though. 


The multiple beam variant of the MLM algorithm is theoretically capable of 
eliminating the difficulties encountered by the single beam in processing a field 
consisting of coherent modes. This is accomplished only at the expense of some 
reduction in performance, and is implemented by simultaneously ineluding in the 
multiple beam MEM beamformer all modes with stgmiffeant coherent contributions. 


MS a practical investigation indicates that the PRAM IV. vertical array. lacks 
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the fundamental resolution necessary to guarantee the numerical stability of such an 
approach. Thus, neither variant of the MLM algorithm is considered to be adequate 
for the reduction of FRAM IV vertical array data, since (as is shown in Chapter 7) 


the modal contributions turn out to be highly coherent. 


The multiple beam least squares algorithm has been shown to be an adequate 
high resolution alternative to the single beam least squares algorithm for vertical 
array data, because it is insensitive to the issue of mode coherence, and because it 
provides the multiple mode resolution at lower frequencies that is lacking in the 
single beam variant. The major issue involved in implementing any such algorithm 
is the selection of the number of modes to model (or. equivalentlv, the number of 
modal beams to estimate simultaneously). Inclusion of t06 many modes that are 
fundamentally unresolvable by the array causes numerical stability problems in the 
matrix inversions required to implement the processing. A quantitative method for 
selecting the number of modes to be include in the multiple beam least squares 
beamformer is presented. This method is based on the relative singularity of the 
steermg matrix when different numbers of modes are included. Using the method, 
which has been verified through simulation, Table 5-I has been compiled: it 
indicates the number of modes that can be simultaneously resolved at different 
frequencies when the multiple beam least squares processor Is apphed to data taken 
from the FRAM IV vertical array. The results of Table 5-I indicate that. multiple 


modes can be resolved for all frequencies of interest. 
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Chapter 6 


The Array Tilt Problem 


The signal processing aspects of modal beamforming have been discussed m 
the previous chapter. For the purposes of that discussion, the array has been 
assumed to be ideal in all respects. In particular, the locations of the sensors that 
comprised the array have been assimed to be known exactly. In practice, such 
sensor positions are rarely known precisely; this is certainly true for the FRAM IV 
arrays, as Is discussed in Sections 3.5 and 3.7. The purpose of the present chapter Is 
to investigate the sensitivity of the modal decomposition process to the types of 


sensor position errors likely to be encountered in practice. 


The chapter begins with the selection of the hnear tilt model for the shape of 
the FRAM IV vertical array. The model is very convenient, since only one 
parameter, the effective array tilt angle, is needed to completely specify It. even 
when the three-dimensional aspects of the problem are considered. Such an 
approach can be justified in two ways. First, it accounts for the majority of the 
sensor offsets found in the notional shapes presented in Section 3.7, especially at 
larger tilt angles, where both the size and the effect of sensor position errors are 
greatest. Second, the physics of the problem are consistent with the conjecture that 
higher order array shape variations have less of an Impact on the modal 


decomposition process than does the linear tilt. 


Next, the tilt problem is studied so that one can qualitatively understand its 
effect on vertical array modal beamforming. This study is supported by results 
generated from synthetic data which depret graphically the ult angle sensitivity of 


modal beamformers. The physical msight generated also leads to a simple method 


? 
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for estimating the accuracy cf tilt angle measurements needed for valid modal 
separation to be guaranteed. This method is then applied to the FRAM IV vertical 
array and sound channel, and the results are both reported and compared with the 


synthetically generated outputs. 


The discussion Just outlined contributes in three ways. First, it emphasizes 
the general importance of vertical array tilt in the modal decomposition problem. 
Second, it provides at least a rudimentary method of estimating the quality of tilt 
angle measurements that are needed. Finally, it proves that the array tilt problem 
must be addressed for successful application of modal beamforming to the FRAM IV 
vertical array data. However, no tilt angle measurements are directly available 
from the experimental data package. In light of this situation, one must then 


attempt to recover estimates of the effective array tilt from the acoustic data. 


The remainder of the chapter presents a simple method for deducing the 
effective array tilt from the acoustic data set. This method mvolves the 
minimization of the residual total square error over a range of reasonable tilt angles. 
The residual total square error 1s a measure of how well the sound pressure field 
observed at the various array sensors is matched bv a pressure field. reconstructed 
from the complex modal amplitude outputs of the beamformer. The accuracy of 
the method is investigated with synthetic data, and some conclusions are drawn 


about the conditions required to obtain reliable results. 


6.1 A Simplified Model for the Shape of the Vertical Array 


Two estimates of the true shape of the vertical array in the presence of 
reasonably severe relative currents have been presented in Section 3.7. "The 
methods used to generate these estimates can be considered mathematical models of 


the vertical array shape. Both of the models are relatively crude, even to the point 
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of being somewhat unrealistic. From a sıgnal processing point of view. though. still 
further simplification is desirable, so that the number of parameters needed to 
define the shape of the array can be kept to a minimum. In the remainder of this 
thesis, the array ıs modeled as a simple straight line that makes some angle with the 
true vertical. In terms of the signal processing, this is the simplest model possible, 
since it requires only one parameter, the effective array tilt angle. to completely 


determine the model array shape. 


The selection of the simple linear tilt model for the vertical array shape can be 
justified by examining how well it fits the actual shape of the array under different 
conditions. Since no actual array shape data is available, this is difficult. to. do 
directly. Instead, the comparison is made with the NRL array model on the 
assumption that the model is at least a reasonable reproduction of the general shape 
of the array. The fit can be quantified by breaking the NRL result presented in 
Figure 3-5 into polynomials involving increasing powers of z. The natural basis set 
for such a decomposition is the set of normalized Legendre polynomials, since they 
are the orthogonal function set involving increasing powers in z. The coefficient of 
the first order term in the expansion represents the assumed tilt angle that the 
arrav makes with the vertical, while those of higher terms indicate the relative 
contributions that more complex curves make to the total array shape. The more 
commonly used Tavlor series expansion is misleading in this situation. since the 
Tavlor polynomials are not orthogonal. The coefficient of the first order Taylor 
term does not reflect the full tilt of the array. since odd Tavlor polynomials of 
higher order (A >, etc.) all contain additional implied array tilt. 


th 


From orthogonal function. theory [20]. the coeffietent of the 7° term of the 


series can be computed as 


ah 
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where z' is a normalized depth coordinate in the range from -1 to I. being related to 
- № (€ € - - | € - « сә ~ - ` a t 


the unnormalized depth x and the total length of the array // by 


лз A ү. (6.2) 


The notation r'(z' is used to indicate the actual array offset as measured at the 
equivalent unnormalized depth. Pi) is the г“ Legendre polynomial: the lowest 
order functions of this set are provided in Table 6-I. In addition, the sanie theory 


guarantees that 


rS = e 
[=> / СЕЕ > U, (6.3) 
іш- 0 


so that the percentage contribution of any term to the full array shape can be 


computed by dividing the square of the appropriate coefficient by Z. 


Table 6-I provides the results of the Legendre polynomial decomposition for 
the first few terms of the series. It is easy to see that the series can be reasonably 
well approximated with as few as the first two terms. Since the zero order term 
represents nothing more than a translation of the horizontal axis, tt has no practical 
effect on the beamforming, and can be completely eliminated by proper redefinition 
of the coordinate system. However, even if the zero order term Is ignored, the first 


order linear tilt term still accounts for about 90°C of the remaining array offset. 


Thus, the conclusion is reached that the array shape displaved m Figure 3-5 
can be adequately represented by an array that remains a straight Ime but which is 
tilted at some angle to the vertical. Whether or not the array can normally be 
represented by a linear tilt is more problematie. A quantitative discussion of this 
point m clearly impossible. due to the lies of experinrental array shape 


measurements. However, two general arguments can be made in favor of such a 
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Relative Contribution 





ре ове: 3/2 К | 
Polynomial Coefficient (m*/ with dy without d 
В) = \/0.5 9.82 x 107 80.02% 
Ри“) = у15 = 4.62 X 10° 17.72%% 88.68% 
Pe) = \/2.5 (1.5 z ^ — 0.5) -1.65 x 10° 2.26% 11.32% 


Р. (=!) = 4/3.5 (2.5 218 – 1.5 и) 3.35 x 10° : Е 


I— 1305 x 109 m3 


Table 6-I: Results of Legendre Polynomial Analysis Results for 
the NRL Model of the FRAM IV Vertical Array Shape 


conclusion. First, the structural mechanics of the vertical array suggests that the 
lowest order shape components should dominate. For the length and weight 
involved (1000 m and roughly 2 tons). it is difficult to visualize a condition where 
much more than a very few low order terms of the expansion contribute 
significantly. This supposition agrees with the results presented in Table 6-1. where 
anything beyond the parabolic term is obviously negligible. Second, it must be 
realized that the accuracy of the linear tilt model can be expected to generally 
increase at higher currents, Where the array shape effects are larger and their 
impact on the problem is more significant. As the array becomes strung out at 
these currents, the proportion of the array shape provided by the tilt term increases. 
Thus, the model is most accurate in the worst case condition, where tt is needed 
most. Finally, it is posstble to argue that it is necessary to retain the second order 
term for accuracy, but such a decision must be balanced against the increased 


complexity of the resulting model. In the present case, the parabolic term appears 
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to be small enough that more ıs gained by dropping ıt than by retaining it. 


The true array shape ıs a three-dimensional function, as shown in Figure 3-5. 
However, it is only one projection of the full shape (that in the vertical plane of 
sound propagation) that affects the modal beamforming. For the purposes of the 
Legendre polynomial analysis just conducted, this vertical plane was arbitrarily 
assumed to be coincident with the vertical plane of current flow at the surface, 
which will not generally be the case. In a more general situation, the Imear tilt 
model can be thought of as the first order Legendre polynomial fit to the projection 
of the array shape into the vertical plane of propagation. The tilt angle associated 
with this linear tilt is then not a true array tilt, but only an effective one measured 


in the direction of propagation. 


6.2 The Sensitivity of Modal Beamforming to Array Tilt 


Having demonstrated that the first order effect of sensor offsets in vertical 
arrays Is to provide an effective tilt to the array, the sensitivity of the modal 
decomposition process to this tilt must next be considered. Consider the simple 
channel shown in Figure 6-1, which consists of a free surface at = = 0, a hard 
bottom at a depth of z = №, апа а sound speed c, that is constant throughout the 
channel. For simplicity, assume an array of continuous aperture that extends to the 
full depth of the channel, and a situation where only the i" mode is present. For 
such a case, the multiple beam least squares result of Chapter 5 reduces to the 
modal equivalent of the conventional beamformer. The resulting processor may also 
be interpreted as a traditional matched filter, but in a spatial sense. Thus. to detect 


h 


4 . . . 
the / ^ mode, the output of the array at any depth should be werghted by the size of 


the mode shape at that depth. which for thts channet ts 
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Figure 6-1: Modal Beamforming for an Untilted Array in a Simple Channel 
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Here v; is the vertical wavenumber of the mode, given by 


| (i - 1/2) 2 e 
n = re (6.5) 


The second expression in equation (6.4) provides another interpretation of 


modal beamforming. The weighting of the array aperture by the mode shape № 
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Figure 6-2: Modal Beamforming for a Tilted Array ina Simple Channel 


equivalent to simultaneously forming two plane wave beams, one at angle 7, to the 
horizontal and the other at angle —3,, where 


pue С 
МЕ = () m 
B, = sin е N. (6.6) 


2 2 
274, a 


The outputs of these two beams are then subtracted and properly scaled to form the 
modal amplitude estimate. In equation (6.6) c, is the horizontal phase speed 


associated with the jth mode, Which is related to the vertical wavenumber by 


Ae "eq 
С, == E = Er. : ) (6.7) 
ЛГ 24, 
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This interpretation is not surprising, since the equivalence of each mode with a pair 
of plane waves is well established for this particular channel. 


Because the modal beam is no more than the sum of a pair of plane wave 


h 


beams for this case, the ¿'* modal beam can be expected to respond to energy 


traveling in either the upward or downward direction at a grazing angle of 5; 
further, this is true whether or not the energy is really being carried by the ү! 
mode. Now, consider the effect of imposing a linear tilt on the array, which ts 
illustrated in Figure 6-2. The tilt causes the plane wave beam pairs to be aimed 
away from the directions in which the modal energy is actually arriving. In fact. if 
the array tilts too much, then one or both of the equivalent plane wave beams can 
end up pointed in a direction where it detects the arriving energy of a neighboring 
mode instead of the energy carried by its own mode. Since it is obvious that the 
modal beamformer responds to this energy, one can readily anticipate severe 


distortions in the modal decomposition if the effective array tilt angle becomes too 


large. 


It is of interest to consider the effect of higher order shape terms from a 
similar viewpoint. Whereas the linear tilt term tends to redirect the equivalent 
plane wave beams, the presence of higher order shape terms causes decreasing 
amounts of mis-steering In the beam formation and increasing amounts of beam 
defocusing. Thus, the hnear tilt approximation is probably better suited to the 
array shape modeling problem than expected, since it retains the portion of the 


actual array shape that most significantly affects modal decomposition. 
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Figure 6-3: Sensitivity of Vertical Array Modal Beamforming to 
Array Tilt (Synthetic Data at 30 Hz - First Realization) 
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Figure 6-3: Sensitivity of Vertical Array Modal Beamforming to 
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Figure 6-4: Sensitivity of Vertical Array Modal Beamforming to 
Array Tilt (Synthetic Data at 30 Hz - Second Realization) 
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Figure 6-4: Sensitivity of Vertical Array Modal Beamforming to 
Array Tilt (Synthetic Data at 30 Hz - Second Realization) 
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Figure 6-5: Sensitivity of Vertical Array Modal Beamforming to 
Array Tilt (Synthetic Data at 47 Hz) 
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6.3 Examples of the Array Tilt Sensitivity in the Processing 


To highlight the sensitivity of modal beamforming to the assumed tilt angle, 
several synthetic data sets have been created and processed with the multiple beam 
least squares algorithm described. in Chapter 5. The synthesized sound fields 
examined are similar to those used in the previous chapter, since each signal field 
has been generated by coherently summing contributions from the odd modes 
between mode 1 and mode 15. Each mode has been assigned a unit amplitude and 
a random initial phase. Unlike the synthetic data presented in Chapter 5, though, 
the signals analyzed here have been created with constant, unmodulated modal 
amplitudes. The array has been assumed to be precisely vertical for the full leneth 
of each synthetic signal. The multiple beam least squares algorithm used to 
accomplish the beamforming is identical in all respects to that emploved in the 
previous chapter. The window length and processing interval selections are also the 


same, 


Figures 6-3 and 6-4 illustrate two different examples of the results for 
synthetic data sets at 30 Hz, while Figure 6-5 provides a similar example at 4% Pz. 
To develop these figures, the modal amplitudes for a single processing window of 
data have been computed at a series of different assumed vertical array tilt angles. 
The resultant modal amplitude outputs have then been plotted as a function of the 
effective tilt angle assumed in generating them. As in Chapter 5, the vertical axis 
represents the mode number to which the amplitude data applies. Since the mode 
number axis really ranges over a discrete set, actual amplitude data occurs only. at 
the integer mode numbers, where the black horizontal hnes have been drawn. The 
contours display the variation of amphtudes on the mode-tilt plane for one 
particular window of data. Results are provided for all 5 of the modes resolvable by 


the vertical array at 30 [2 and for the 7 modes that can be distinguished at 47 Iz. 
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The horizontał axis indicates the effective tilt angle assumed in the beamforming. 
Thus, the vertical column of data at 0° indicates the true distribution of energy in 
the mode set (since the synthetic field is generated for a truly vertical array). On 
the other hand, the column of data at 1° indicates the modal distribution generated 
by the beamformer if 1° of mismatch is present between the actual tilt and the tilt 


assumed by the beamformer. 


All three plots vividly illustrate the variation in modal beamformer outputs 
that can occur when even slightly different vertical array tilts are assumed. As 
expected, the O° outputs accurately reflect the actual modal distribution present in 
the synthetic field. It is easy to observe the expected alternating arrangement of 
deep nulls (at even modes) and strong peaks (at odd modes) at this tilt. Although 
one is unlikely to encounter the full range of assumed tilt angles (+10 ) in practice, 
it is not unreasonable to expect excursions of as much as £5 . The large variations 
of the modal amplitude estimates even within this smaller range of tilts show that 


knowledge of the actual tilt angle is necessary for valid modal decomposition. 


The synthetic data sets used to generate Figures 6-3 and 6-4 differ in only one 
respect. that being the inttial phase relationships assumed between the varıous 
modes making up the signal field. Both figures are included to demonstrate that 
the pattern of peaks and nulls generated is highly dependent on these phase 
relationships. A comparison between either of these figures and Figure 6-5 shows 
that, as might be expected, the processing becomes more sensitive to tilt angle at 
higher frequencies. Amplitude estimates for higher modes also exhibit) greater 
sensitivity than those for lower modes at the same frequeney. If one wishes to 
assure reasonable accuracy in first mode amplitude estimates. tt appears that. the 
array tilt must be known to roughly the nearest I. Estimates at higher. modes 


require much more precise tilt data to be valid. 
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6.4 The Array Tilt Accuracy Required for Valid Modal Estimates 


Given the sensitivity of the vertical array processing to assumed array tilt, 
estimation of the accuracy required in the tilt angle measurement to guarantee valid 
modal decomposition becomes important. One simple approach to this problem is 
outlined here. From Figure 6-2, it is obvious that, for a simple uniform channel, the 
modal separation is certainly invalid by the time the array has been tilted far 
enough to cause one of the pair of equivalent plane wave beams to point in the 
direction of a neighboring mode’s arriving energy. Adopting a criterion for valid 
results of one half of this tilt leads to the expression 


C с 

ES ( S 
В о | | (6.3) 
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where the second half of equation (6.6) has been used to evaluate each of the 
grazing angles. The minimization must be conducted over all modes to be properly 


estimated. 


Equation (6.8) 1s strictly correct only for a hard bottom, constant sound speed 
channel, where the equivalence of the mode to a pair of plane waves ls exact, 
Applying it to more general channels requires further justification. Typically., each 
mode is equivalent to a full spectrum of plane waves across vertical wavenumber 
space, rather than just the discrete pair considered here. This spectrum can be 
calculated by taking the Fourier transform in depth of the mode shape. However, 
mode shapes often exhibit a sort of characteristic vertical wavenumber, as can 
readily be observed by studying Figures 1-6 and 1-7. In such a case, the resulting 
wavenumber spectrum сап be expected to show a sharp peak near this 
characteristic value. It is not unreasonable to approximate such a spectrum by а 
pair of delta functions; which then completes the identification with an equivalent 


plane wave pair. Yang and Giellis [84] provide several examples of modal spectra 
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that confirm the presence of these strong wavenumber peaks in the central Arctic 


mode set. 


There remains the issue converting the characteristic vertical wavenumber 
value for each mode into an equivalent grazing angle. This is trivial for a channel 
of uniform sound speed, since it is easy to define a horizontal reference; any mode 
propagating in exactly the horizontal direction must possess a horizontal phase 
speed of cy. The choice of a proper definition for the horizontal reference in a 
channel of arbitrary sound speed profile is much more subtle. Care in this matter i5 


important, since the form of equation (6.3) makes it sensitive to the choice of v. 


The problem is often dealt with by making a "relative" conversion: that is, 
the difference in grazing angle between two modes is computed by arbitrarily 
assigning the role of cy to the mode having the lower horizontal phase speed. A. 
simple numerical example is enough to demonstrate the flawed nature of such an 
approach. Assume that the phase "e of the first mode is 1450 m/sec at a given 
frequency, while that of the second mode is 1455 m/sec. Also assume that the 
reference phase speed associated with horizontal propagation is 1410 m/sec. Then 
the actual difference between the equivalent grazing angles of the two modes ts 
1.5°, whereas a value of 4.8° is given by the relative conversion method. While 
this level of accuracy may be acceptable for problems to which the relative 


conversion is typically applied, it is obviously inadequate for the present purposes. 


A definition of c (the horizontal phase speed associated with propagation at a 
grazing angle of 0°) is needed that is general enough to be applicable to a channel 
having an arbitrary sound speed profile. Such a definition can be developed by 
identifying horizontal propagation with the mode having the smallest. horizontal 
phase speed. The mode possessing this property is always the "DC" mode. Le., the 


lowest order mode for the desired channel when both уун surface and rigid 
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Figure 6-6: Horizontal Phase Speed for the DC Mode versus Frequency 
bottom boundary conditions are assumed. It is named the DC mode here because 
the resulting mode shape is constant with depth for a channel of uniform sound 
speed, and, thus, represents the DC component of the resultant Fourier series. Note 
that this definition is consistent with the exact result derived above, since the phase 
speed of the DC mode for a channel of uniform sound speed is exactly ege The 


yhase speed of this mode for tlre sound speed profile of Figure 4-1 is plotted versus 
(=> 


frequency m Figure 6-6. 


Table 6-II has been computed from equation (6.8) and Figure 6-6. [t indicates 
the accuracy required in the array tilt measurements 1f. valid. modal amplitude 
estimates are to be assured for the modes resolvable in the FRAN IN vertical array 
data set. It can be seen that accuracies of better than 0.25 are often needed. 


These results agree well with the tilt angle sensitivities observed in Preures 6-2 
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through 6-5, where significant changes in the modal amplitude estimates can ђе 
observed over tilt changes of much less than 1°. The higher tolerance of the first 
mode estimate to tilt angle variations that was noted earlier appears to be due to 
the greater separation in phase speed between it and the other modes; this causes an 
equivalent increase in the tilt angle needed to alias the resultant amplitude estimate. 
A similar phenomenon can begin to be seen for the second mode at higher 
frequencies. Note that the accuracy requirements generally increase at both higlrer 
frequencies and higher mode numbers. This is in keeping with the observations 


made from Figures 6-3 through 6-5. 


Frequenev 1 Моде 2+ Modes 

15.00 Hz 0.45. 0.40 

ІШ Boule 0.50 | 0.35 

20.00 Hz 0.65 0.30 

ОБОН? 0.90" . 025° 

30.00 Hz 1.35 0.20 ° 

Dolo LZ 1.60” uS 

50 Н? 1.60` 0.15 
Frequency ] Mode 2 Modes 3+ Modes 

47.00 117 1.20 0.35 0.10 

29 795 H7 1.80 0.50 0.10 

55.00 Hz 1.20 0.55 0.10 
Frequency 1 Моде 2 Modes 3 Modes 44 Modes 
11.00 HHz 1.60 0.90 0.20 0.07 


Table 6-II: Tilt Angle Accuracy Requirements for the FRAM IV 
Vertical Array 


From Table 6-If, it can be seen that the tilt accuracies required to assure 
proper modal decomposition of the received field are significantly smaller than the 
range of tilt angles likely to be encountered im the field. Thus. very accurate arrays 
tilt measurements must be made either directly or indirectly if a valid modal 


beamforming experntent 1s to be conducted. The FRAÄM IN vertical array tilt 
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accuracy requirements are probably typical of those found elsewhere, while in other 
situations somewhat larger ranges of tilt angle will be encountered. Therefore, it is 
reasonable to conclude that vertical array tilt monitoring must be considered as an 


important integral element of any reliable modal decomposition scheme. 


6.5 The Effect of Tilt on Horizontal Arrays 


The sensitivity of the vertical array to array tilt is generated bv a confluence 
of factors; these factors tend to work in just the opposite fashion for horizontal 
arrays, making them relatively insensitive to array tilt issues. The first order effect 
of ocean currents on a vertical array is the generation of a tilt across the face of the 
array. By comparison, the first order effect for a horizontal array ıs a uniform 
offset of sensors across the array, which has no influence on array performance for 
far field beamforming. Some variation of individual sensor offsets is expected. but 
this is generally small when compared to the mean offset. Such variations may 
cause some degradation in horizontal array beamformer performance at higher 
frequencies, Where exact sensor position information Is more critical. Even the 
fluctuations around the mean offset generate no real array tilt, since they tend to be 
randomly distributed across the array. Finally. because a horizontal array is 
essentially being operated at endfire when used for modal decomposition. the 
algorithm itself can be expected to be relatively insensitive to tilt. However, it 
should be realized that it is exactly this endfire operation that limits the modal 
resolution of horizontal arrays in the first place. In general, then, array tilt does not 
appear to be as significant an issue for the horizontal arrays as it is for vertical 


NS. 
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6.6 Recovery of Vertical Array Tilts from the Acoustic Data 


The conclusion that even small amounts of array tilt are important im the 
modal decomposition problem has major consequences for vertical array design. In 
a different way, it also is significant for the FRAM IV vertical array data to be 
analyzed here. Since no field measurements of verticai array tilt were taken during 
the experiment, one must devise a scheme to extract the effective array tilt froni the 
acoustic data if further progress is to be made. The present section presents such а 
scheme for recovering the effective tilt angle. Such an approach is obviously not 
ideal; its ability to provide even approximate results is mtimately linked to thre 
underlying quality of the data set. In the present instance, the high signal to noise 
ratio developed by the preprocessing (typically 20 to 30 dB) and the exceptional 
stability of the received signals provide some hope that some reasonable tilt angle 


estimates can be made. 


To estimate the effective array tilt from the acoustic data, the residual error 
from the multiple beam least squares beamformer is evaluated as a function of the 


tilt angle assumed in the beamforming process. Equation (5.7) may be rewritten as 


NORA ( | 


“тіп 


Pe TENES EN EM). | 
<p Pp  -— (6.9) 


where the dependence ou the assumed tilt angle 4 has now been explicitly indicated. 
Additionally, the error in equation (6.9) has been normalized by the total square 
pressure present in the observed field, P* P, so that it now represents a fractional 
goodness of fit measurement. When Qu V) is evaluated across tlie range of 
physically reasouable tilt angles, a consistent global mimmmum can typically be 
expected. The tilt angle that generates ihis minimum can then be used for the 


estimated tilt of the array. The underlying assumption is that the beamformer ts 
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better able to match the observed sound pressure field when using the correct array 


tilt angle than when working with any incorrect tilt value. 


It is important to review what the procedure just outlined accomplishes from a 
processing standpoint. The multiple beam least squares modal beamformer 
computes the set of complex modal amplitudes that produce the best fit to the 
observed pressure field for one assumed array tilt. The residual error Is a measure 
of how well this "best fit" actually matches the observed field. Note that the 
residual error is always smaller than the total square error generated by any other 
set of modal amplitudes at the assumed tilt angle. If the data is then processed 
using the tilt angle which minimizes the residual error, the resulting set of 
amplitudes is guaranteed to produce a better match to the observed field than any 
other possible combination of modal amplitudes and array tilts. Thuas. the 
minimized error philosophy, which forms the basis for the original beamforming 


algorithm, has been extended further to cover the array tilt angle. 


Synthetic data has once more been employed to study the performance of the 
scheme. One such result is displayed in Figure 6-7. Here the normalized residual 
error has been plotted as a function of assumed tilt angle and time. The normalized 
error has been plotted on a dB scale, so that the -10 dB line marks a 90% fit and 
the -20 dB line indicates a 99*6 fit, ete. As a practical matter, the normalized 
residual error has been computed as an average value over a two minute window of 
data, rather than instantaneously (as indicated in equation (6.9)). The computation 
has been made once each minute, giving a 50% overlap factor. This windowing 
procedure Is identical to that carried out for synthetic data throughout the thesis. 

The synthetie data set from which Figure 6-7 has been ereated is similar in 
most ways to those presented previously. For the present figure. modes 1. 3, 5. and 


1 of the 17 Hz mode set have been coherently summed to generate the signal field: 
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Figure 6-7: Normalized Residual Error as a Function of Tilt Angle 
for Synthetic Data at 47 Hz (No Excluded Modes Present) 
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Figure 6-7: Normalized Residual Error as a Function of Tilt Angle 
for Svnthetic Data at 47 112 (Хо Excluded Modes Present) 
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Figure 6-8: Normalized Residual Error as a Function of Tilt Angle 
for Synthetic Data at 47 Hz (Excluded Modes Present) 
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Figure 6-8: Normalized Residual Error as a Function of Fih Angle 
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02 
each mode has been assigned an amplitude of unity. Unlike the synthetic.data sets 
used in Chapter 5, the modal amplitudes shown here are constant and unmodulated 
over time. Rather, in Figure 6-7 it is the effective array tilt that has been 
sinusoidally varied over time. This variation occurs over a range from vertical to 
5° and back, and has an oscillation period of 3 hours (10,800 sec). Figure 6-7, 


therefore, covers the first half cycle of this tilt oscillation. 


Figure 6-7 demonstrates the ability of the tilt angle estimation scheme 
described above to properly track vertical array tilt over time. The minimum 
errors observed are in the -15 dB range, representing a 97°¢ fit of the synthetically 
generated data. Note that the minimum error tracks the actıral tilt of the array 


quite well as it swings away from the vertical and then back. 


Figure 6-7 1s a highly idealized situation, for two reasons: a tilted line array 
possessing no higher order shape components has been assumed in generating the 
synthetic field; and only the seven modes included in the beamformer have been 
injected into the signal. While the effect of the presence of higher order shape 
components on the array tilt estimation process is difficult to quantify, the 
fundamental Issue underlying it remains the adequacy of the linear tlt model for 
array shape, which has already been discussed at some length. Figure 6-8 
demonstrates the effect of having significant amounts of signal energy carried бу 
modes not included in the beamformer. Here, both the synthetic data and the 
processing used are identical to those of Figure 6-7, with the lone exception that 
modes 9, 11, 13, and 15 have additionally been inserted into the signal field at unit 
amplitude. The presence of this extra energy m modes excluded) from the 
beamformer makes rt more difficult for the algorithm to attam a good match with 
the observed field at any tilt angle. which is reflected bv. the increased: minimum 


restdual error 0 ІП Figure G-S Versus Qu In Figure 6-4). It also ar РОМ | bias 
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of about 2° in the resulting tilt angle estimate. While both the amount and 
direction of the bias are difficult to predict, its cause is not hard to explain. Even 
though some energy in the included modes 15 not accounted for, the biased tilt angle 
generates less residual error than the true tilt because even more of the energy in 
the excluded modes is allowed to leak into the result. The amount of the bias ts 
directly reiated to the relative proportions of signal energy present in the modes 
included and excluded by the beamformer. The actual minimum error level 
achieved can be used as a measure of this effect. Processing that results in 
minimum normalized errors of much above -10 dB should be regarded with 


suspicion. 


6.7 Summary 


The least squares modal beamforming algorithm developed in Chapter 5 has 
been shown to be quite sensitive to sensor position errors, particularly when it is 
applied to vertical arrays. The effective linear tilt of the array away from true 
vertical appears to be the component of array shape that contributes most 
significantly to this sensitivity. Two reasons account account for the effect. The 
shapes of long, heavy vertical arrays (such as the one deploved during the FRAM IV 
Arctic Experiment) are composed, in a mean square sense. primarily of low order 
components, such as constant translations, linear tilts, and parabolic shapes. The 
zero order constant translation term has no effect on the beamforming problem., and 
may be completely removed by proper redefinition of the coordinate system. The 
first order linear tilt term can be expected to account for the majority of the total 
sensor position error, especially in strong currents where the sensor offsets are 
greatest. To demonstrate this concept. a model-based notional shape generated for 


the FRAM IV vertical array has been analyzed using Legendre polynomial. Also, 
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. the physics of the problem support the argument that the impact of the linear tilt 
component of the array shape on the modal decomposition process 15 тоге 
significant than that of higher shape components. This is because the array tilt 
causes each modal beam to effectively "point away" from the vertical directions in 
which that mode’s energy arrives. By contrast, the primary effect of higher order 


shape components is only to defocus the modal beams. 


The above considerations have lead to the selection of a simple tilted line 
model for the true array shape. Such a model is also attractive because it. allows 
the assumed shape to be completely specified with a single parameter, the effective 
array tilt angle. This simplicity holds even when the three-dimensional nature of 


the problem is considered. 


The sensitivity of vertical array modal beamforming to array tilt has been 
studied both theoretically and through the use of synthetic data. A simple method 
has been developed to estimate the accuracy to which one must know the tilt angle 
if valid modal decompositions are to be assured. This method implicitly assumes 
that a rough equivalence exists between each mode and a corresponding pair of 
plane waves. The applicability of the resulting accuracy estimates is dependent 
upon the degree to which this equivalence holds in practice. The equivalence may 
be tested qualitatively by examining the Fourier transform of each mode shape: if a 
single strong peak exists over some narrow range of vertical wavenumbers, then a 
reasonable equivalence exists. This appears to be the case for the central Arctic 
mode set of Chapter 4. Processing results from synthetic data have been utilized to 
actually demonstrate the sensitivity of the processing to array tilt. 

Application of this tilt accuracy estimation method to the FRAME IV vertical 
array shows that the array tit must be known to about +£0.25 to guarantee 


meaningful modal separation. These results are in good agreement with tilt angle 





sensitivities demonstrated by synthetic data. They also correspond well to the tilt 
sensitivities observed in real data taken from the vertical array, which are presented 


in Section 7.3. 


Since no field measurements of vertical array tilt were made during the 
FRAM IV Experiment, a scheme has been developed to recover tilt values from the 
acoustic data itself. This scheme involves evaluation of the residual beamforming 
error as a function of the tilt angle assumed in the beamforming compntatitons. The 
restdual error is a measure of the difference between the the sound pressure field 
observed at the sensors and the field synthetically generated from the complex 
amplitude outputs of the modal beamformer. As such, ıt indicates how well the 
outputs of the beamformer match the observed pressure field. Thus. the tilt that 


generates the smallest residual error is the one which best fits the observations, and 


may be assumed to be the actual vertical array tit angle. 


Through the use of synthetic data, the tilt angle estimation scheme has been 


shown to work well when applied to high SNR signals. In making this 


5 
demonstration, an ideal linear shape has been assumed for the array. The low SNR 
situation has not been studied, but it is expected that performance would be 
severcly degraded under these conditions. The relative amount of signal energy 
carried by modes not included in the multiple beam least squares algorithm ах 
been shown to affect the accuracy of the estimate. When the contribution of these 
modes is negligible, the method can be expected to provide accurate estimates: as 
the contribution becomes more important, a bias m injected into the ем паге, 
leading to progressively poorer performance. When significant signal energy Is 
resideut ш modes not included іш the beamformer, higher values of the normalized 
residual error are also generated, even at the minimum. This is because the 


beamformer has increasing difficulty mm matching the total freld with just the modes 
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available to it. Thus, the value of the normalized residual error at the minimum 
can be used as a guide to the accuracy of the tilt angle estimate. Fits of much less 


than 90€6 (greater than -10 dB normalized error) should be considered to be 


suspect. 
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Chapter 7 


Data Presentation and Analysis 


The purpose of this chapter is to present the results obtained from applving 
the processing techniques outlined in Chapters 5 and 6 to the data set described in 
Chapter 3, based on the modal structure of Chapter 4. In so doing, a number of 
different hypotheses can be tested. First, the vahdity of the assumed modal 
structure can be established. Second, the results of processing the vertical array 
data can be compared to those generated from the more traditional horizontal array 
processing, providing a qualitative measure of the vertical array processing 
performance. These efforts generate confidence in the computed mode propertics 


and the processing performance, after which a full analysis of the implications of 


both the horizontal and vertical array results can be undertaken. 


The processing output includes three different types of useful data. Тһе 
absolute mode amplitudes are directly available. Their relative sizes provide 
immediate insight into the nature of low frequeney Arctic propagation, and they can 
be used to compute observed values of the modal attenuation. coefficients. 
Comparison of these coefficients with. predicted. values sheds light on the channel 
parameters that most influence Arctic sound transmission. The time series of the 
amplitudes also provide some general indications of the temporal stability of the 


central Arctic channel. 


Measurements of mode coherence provide data of a different type. and 
indicate. whether the stgnals component modes are phase-random or phase-locked 
with respect to each other. This data can be used to test the uncorrelated muki- 


path assumption that i$ widely used in theoretical developments. 
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Finally, if the different modes do maintain a deterministic phase relationship, 
then it is meaningful to measure the relative phases between various pairs of modes. 
These relative phase values provide a measure of the range dependence of the 
channel, since they represent an integration of horizontal phase speed fluctuations 
that the modes have encountered during propagation (see equation (2.30)). Thus, if 
the measured values disagree with values predicted from the phase speed data of 
Figure 4-10, then the channel must be assumed to show non-negligible variations i 
range; if the agreement is good, the simpler range independent channel model will 


suffice. 


This chapter is organized into five major sections, some containing niultiple 
subsections. Following some preliminary comments, the horizontal array data 
processing results are presented and summarized. The section is divided into four 
subsections: one deals with azimuthal DIS and another with vertical 
steering. The third summarizes the results across the full data set. In the last 
subsection. the predicted first mode horizontal phase speeds are compared with 


those observed from the horizontal array data, so that the accuracy. of the modal 


characteristics computed in Chapter 4 can be assessed. 


A similar section summarizes the results for the vertical array. Two 
subsections discuss the array tilt problem and the tilt angle estimates made with the 
procedure discussed tm Section 6.6. The third presents some typical modal 
amplitude levels versus time, while the fourth discusses observations of possible 


array shape changes over time that can be observed in the data. 


Following this are three sections dealing with data analysis. The first of these 
sections Investigates the modal amplitude data. One subsection. compares. tlie 
| 


horizontal and vertical array results. The second analyzes the observed modal 


attenuation coefficients versus theoretical predictions, and the last subsection 
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contrasts both sets of results with some shot data available for the same channel. 
The second section presents and discusses the mode coherence estimates made. The 
final section presents the relative phases observed between various modes. A 


summary of major points completes the chapter. 


7.1 General Comments 


Before commencing the actual presentation of processing results, some general 
comments concerning the formats used and the processing parameters selected are 
in order. First, all amplitude levels shown in this chapter are displayed on a decibel 
scale referenced to units of absolute pressure (1 ¿Pa rms); conversion from the 
voltage levels measured has been done using the hydrophone sensitivities quoted in 
Figure 3-6 and Table 3-H. A logarithmic scale is also used for displaying residual 
error data; the values shown are 10 times the logarithm of the normalized residual 
error, as defined in Section 6.6. An error level of -10 dB represents a 10%. error, or, 
equivalently, a 90% fit to the observed data, while a level of -20 dB represents a 
99% fit. Phase data are presented on a linear scale from -180 to 180 : phase 
wrap Is often observed, as only the principal value of the phase is displayed. 
Coherence measurements are presented in normalized form, so that the values range 


from 0.0 to 1.0. 


On all plots, amplitudes or residual error values are marked with DB. Моде 
numbers (such as mode 1, mode 2, ete.) are denoted by MD, while horizontal phase 
speeds (in m/sec) are indicated by PV. DG is used for azimuth angles on plots of 
Bonızontalarray data and for array tilt angles on plots of vertical array data. AI 
azimuth angles are measured with respect to the convention discussed in Section 
3.3: the tilt angles displaved are effective angles measured in the vertical plane of 


propagation (see Section 6.1). Where axes are indexed by mode number, it is 





-210- 
important to remember that data only exists at the integer values: any markings 
between these discrete points, such as connecting lines on one-dimensional plots or 
contours on two-dimensional plots, have no physical significance and are included 
only to help in judging relative trends by eye. Where contoured data indexed by 
mode have been presented, overlying black lines have been added to the plot to 
emphasize its discrete nature. All levels plotted versus mode number are modal 
amplitudes. For the horizontal array, levels plotted versus horizontal phase speed 
are plane wave amplitudes, and differ from modal amplitude measurements by a 
factor equal to the size of the mode shape at the 91 m depth of the array. The 


necessary conversion factors can be found in Figure 4-3. 


Processing windows of 28.8 sec (18 points, providing a time-bandwidth 
product of 0.7) have been used for estimating the sensor cross-coherence matrices 
for horizontal array data. One set of beams is computed every 15 seconds. 
providing roughly a 50*6 overlap factor. With some noted exceptions, vertical array 
data have been processed in 120 sec windows (75 points, giving a time-bandwidth 
product of 3.0), with a computation interval of 60 sec, again giving a 50° overlap. 
The verv stable and coherent nature of the Arctic sound channel makes the 


difference negligible. No windowing functions (such as Bartlett) or Hanning 


windows) have been emploved in the processing. 


The 17.75 and 47 Hz frequency data from both the horizontal and vertical 
arrays have been chosen for display throughout the chapter, since they are 
reasonably representative of the results at other frequencies. This is done to allow 


direct comparison of the various processing results while limiting the chapter to an 


acceptable length. Results from other frequencies. are presented when thes 
demonstrate unique features of interest. Тһе most important results are 


summarnzed by tables or plots that span the full range of available frequencies. 





e 


7.2 Presentation of the Horizontal Array Processing Results 


The data taken from the horizontal array have been processed using both the 
single beam MLM beamformer and the conventional (single beam least squares) 
beamformer. In general, the MLM algorithm provides much higher resolution than 
the conventional approach, but does not produce” reliable quantitative 
measurements, because of the presence of coherent interference. Quantitative 
amplitude estimates of observed MLM peaks have been made by reprocessing the 
data using the conventional beamformer. This technique works best when a strong 
single peak is present, and is considerably less reliable when closely spaced multiple 


peaks of roughly equal height are present.. 


7.2.1 Azimuthal Beamforming 


Figure 7-1 displays the output of the single beam MEM processor versus time 
for the 27 Hz tone from the horizontal array data set. For an assumed sound speed 
of 1450 m/sec, the beamformer successively scans through 3860 beams spaced | 
apart in horizontal azimuth angle. Note that although 360 beam outputs are beme 
evaluated at each time step, every beam is stu} being computed individually: a 
multiple beam algorithm would compute a number of beams all with a single set of 


matrix operations. 


For the azimuthal beamforming done here, the MEM processor is preferable to 
the conventional beamformer for two reasons: it provides much higher resolution: 
and it does a significantly better Job of suppressing side lobes. Quantitative level 
measurements are not an issue here. In Figure 7-1, the tonal signal is readily 
apparent at an azimuth angle of 86.14 . A low level sidelobe that is down about by 
about S dB is visible at 270 . This iss a back lobe generated by the north-south leg 


of the array, Which is essentially broadside to the Incoming signal. 
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Figure 7-1: Output of Single Beam MLM Beamformer Scanned 
in Azimuth for 27.00 Hz Tone (Horizontal Array) 
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Figure 7-1: Output of Single Beam MLM Beamformer Seanned 
in Azimuth for 27.00 12 Tone (Horizontal Array) 
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By processing the data with beams spaced 0.1” apart, one can resolve the 
source azimuth angle to about +0.1 , an accuracy not achievable using а 
conventional beamformer. For comparison, the navigation data provide an azimuth 
angle estimate of 86.1° (see Table 3-I). The difference is well within the limits of 
accuracy of the array rotation measurements, which were made daily from the sun. 


Similar agreement is obtained for other horizontal array tones. 


With one exception, the tone is extremely stable over its full thirty minute 
duration. That exception is the fade that can be seen at the 120 sec mark. The 
fade appears to be caused by an instananeous shift ш the phase of the received 
signal, which in turn causes the signal frequency to walk out of the filter passband 
for a short time. In many respects, it is similar to the fade observed by 
Mikhalevsky [57]. Neither the azimuthal beamforming outputs nor the vertical 
steering results (which are shown in Figure 7-2) display any noticeable difference 
between the spatial structure before and after the event. This observation argues 
against the assumption that the fade represents a shift ın transmission path 
structure; some type of interruption of the source seems more likely. Background 
noise variations can be seen to occur in the absence of the signal and during the 
fade. These are processing artifacts caused by variations in the amount of white 


noise used to stabilize the cross-collerence matrix estimate. 


7.2.2 Vertical Beamforming 


Figure 7-2 shows the output of the single beam MEM processor versus time for 
the 27 Hz tone. Unlike Figure 7-1, however, the beamformer now scans in the 
vertical at the previously determined signal azimuth. This is accomphshed by 
stepping the beam’s assumed horizontal phase speed from 1100 m/sec to 1600 m/sec 


i] to m/sec increments. with the higher phase speeds corresponding to larger 
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Figure 7-2: Output of Single Beam MLM Beamformer Scanned 
in Phase Speed for 27.00 Hz Tone (Horizontal Array) 
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Figure 7-2: Output of Single Beam MEM Beamformer Seanned 
in Phase Speed for 27.00 Hz Tone (Horizontal Array) 
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Figure 7-3: Output of Single Beam Least Squares Beamformer Scanned 
in Phase Speed for 27.00 Hz Tone (Horizontal Array) 
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equivalent grazing angles in the beams formed. The signal appears to be composed 
almost entirely of first mode propagation arriving at about 1460 m/sec. Such a 
result is expected, given the analysis of Section 4.5. If the attenuation of the 
various modes is assumed to be equal, the received level of the first mode at AI m 


should dominate its neighbors by 27 dB (twice the value plotted in Figure 4-9). 


p 


Figure 7-3 shows the same 27 Hz data processed in the vertical at the signal 
azimuth, but now using the conventional beamformer. By comparing this figure 
with Figure 7-2, one can obtain a feel for the resolution gains possible with the 
MLM algorithm when it is applied in a carefully controlled fashion. Figure 7-3 is 
useful for measuring the the first mode's average received signal level (about 91 AB). 
and its fluctuation range (which appears to be no more than I dB). The complete 
saturation of the beamformer output by the first mode arrival makes estimates of 


higher mode levels impossible. 


In Figures 7-4 and 7-5, the output of the single beam MLM processor is agai 
presented versus time. Here the 17.75 Hz tone recorded on the horizontal array has 
been processed, and the beamformer has been steered in the vertical at the 
appropriate azimuth angle. The plots indicate some of the practical difficulties 
involved m properly assessing channel stability. Figure 7-4 was generated first, by 
including all 24 hydrophones available from the horizontal array m the 
beamforming. In an effort to assess the natnre of the fluctuation seen about half 
way mto the signal, the preprocessed time series data were carefully reviewed on a 
sensor by sensor basis. Three different sensors (not adjacent to one another) were 
fonnd to exhibit abnormally high signal levels for short periods near the time of the 
disturbance, and have been eliminated as inputs to the beamformer in the results 
displayed in Figure 7-5. The previous fluctuation has now been completely 


Ва КЕ те сличе об «ће Иоса оп the three sensors. i. not. understood. 
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Figure 7-4: Output of Single Beam MLM Beamformer Scanned 
in Phase Speed for 17.75 Hz Tone (Horizontal Array) 
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Figure 7-5: Output of Single Beam MLM Beamformer Scanned 
ın Phase Speed for 17.75 Hz Tone 
(Horizontal Array with 3 Sensors Excluded) 
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However, the signal level variations point more to a mechanical or electrical 
phenomenon than an acoustic cause. Such variations have also been observed in 


other data taken from the horizontal array. 


The results of Figure 7-5 are similar to those obtained at 27 Hz, and once 
again show dominant first mode propagation. This is in qualitative agreement with 
the preliminary analysis done in Section 3.12 for the 17.75 Hz tone recorded on the 
vertical array. The widening of the first mode peak in Figure 7-5 compared with 


Figure 7-2 is attributed to the loss of array directionality at the lower frequency. 


The final signal displayed from the horizontal array data set Is the 47 117 топе 
shown in Figure 7-6. This figure has again been generated from the output of the 
MLM beamformer when steered in phase speed at the signal azimuth. Note the 
marked lack of a strong first mode arrival, which should appear at a phase speed of 
about 1449 m/sec. Instead, the most intense arrival occurs near 1480 m/sec, and 
can be associated with the deepest diving RSR rays of the TRISTEN/FRAM 
channel (see Figure 4-12). Even this arrival is only marginally stronger than those 
seen at other phase speeds, though. The received levels appear to indicate a rough 
equipartition of energy im the various modes; however, the array lacks sufficient 
resolution in the vertical to make a definitive statement about the distribution of 
energy In the higher modes. In any case, the strong first mode arrival seen at lower 


frequencies has now disappeared. 


This result agrees closely with the initial analysis of the 47 Hz tone from the 
vertical array done in Section 3.12, and is quite remarkable when one considers that 
47 Hz lies almost at the peak of the relative advantage curve shown in figure 4-9. 
The theory predicts that the first mode should dominate its neighbors by almost 60 
dB И the attenuation is constant across all modes, Which is in direct opposition to 


the analysis completed both here and in Section 3.42, where it is concluded that the 








Figure 7-6: Output of Single Beam MLM Beamformer Scanned 
in Phase Speed for 47.00 Hz Tone (Horizontal Array) 
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variations of the received signal with depth can only be explained by the presence of 
a number of coherently interfering modes. The 53.25 Hz tone from the horizontal 


array data set indicates similar low levels for the first mode. 


7.2.3 Summary of the Horizontal Array Processing Results 


Several important general observations can be made from the processed 
horizontal array data. An extremely stable channel has been observed at all 
frequencies in the data set below 70 Hz. Fluctuations in beamformer outputs over 
periods as long as 30 minutes rarely exceed 1 dB. The data appear to be separable 
into two or, possibly, three regimes. The data below about 30 Hz show the 
expected dominant first mode propagation. Despite the tremendous relative 
К ntage possessed by first mode's coupling to both the source and the horizontal 
array, the data above 40 Hz show only very weak first mode propagation. In this 
frequency regime, the deep RSR rays of the channel appear to contribute marginally 
more to the total propagation than does any other E While a strong first mode 
arrival is also seen in the region between 30 and 40 Hz, it is perhaps not as 
dominant as might be expected from the results of Section 45. This frequency 
region appears to be a transition zone between the other two propagation regimes, 
The primary source of this conclusion is the data from 35.25 Hz signal, which has 


not been displayed here. 


The vertical resolution of the horizontal array is Imadequate to allow 
separation of the individual modes, even when the higher resolution MEM algorithm 
ts employed. At frequencies below 40 Hz, good quantitative estimates of the 
received level of the first mode can be made because of its relative strength 
compared to other modal arrivals: this same dominance prohibits estimates for 


higher modes at these frequencies. At frequencies above 10 Hz. all mode level 





measurements are of questionable accuracy because of the lack of a dominant 


propagation path. 


The result at 70 Hz is a special case. While an investigation of the complex 
time series output from the preprocessor shows a signal of reasonable strength to be 
present, neither the MLM nor the least squares beamformers produce acceptable 
outputs. The inability to beamform this data effectively is attributed to the sensor 


displacement problem discussed in Section 3.7. 


7.2.4 Verification of the Modal Structure 


The high resolution of the single beam MLM beamformer allows experimental 
observations of the first mode horizontal phase speed to be made for signal 
frequencies below 40 Hz. These measurements can then be used to verify the modal 


structure derived in Chapter 4. Figure 


7-7 summarizes this data. Here, both the 
computed and the observed first mode phase speeds are plotted as a function of 
frequency. The error bars that are plotted with the experimentally determined 
values indicate the approximate width of the first mode peak seen im tlre MIM 
processing output. With the exception of the 27 Hz result, the agreement between 
prediction and observation is more than adequate. The 27 Hz data point is 
obviously out of line with both the computed results and the other experimental 
data, but the cause of this disagreement is not understood. Since the phase speeds 
are much more sensitive to environmental variations than the mode shapes, one can 
conclude from Figure 7-7 that the results presented m Chapter 4 adequately 


describe the actual modal structure at the receiving arrays. 
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Figure 7-8: Conceptual Relationship of the Different Types 
of Vertical Array Processing Outputs 


7.3 Presentation of the Vertical Array Processing Results 


In keeping with the discussion of Chapter 5, the data taken from the vertical 
array have been processed using a multiple beam least squares beamforming 
algorithm. Because this beamformer is a linear processor, no bias issues occur, and 
the output levels are quantitatively accurate. Of the two most critical aspects of 
this processing, the first, which is the number of modes that can be included in the 


beamformer at any particular frequency, has been fully resolved in Section 2.10 (see 
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Table 5-1). The other is the estimation of an effective array tilt angle for each of 


the different signals recorded on the vertical array, the practical details of which are 


discussed at some length here. 


Three different forms of output have been generated for each signal in the 
vertical array data set. One is a plot of the various modal amplitude estimates 
versus time for any assumed tilt angle. The second is a plot of the modal amplitude 
estimates at one particular time versus a range of assumed tilt angles. These two 
outputs can be thought of as being orthogonal cross-sections of a three-dimensioual 
amplitude function, the axes representing time, mode number, and array tilt. This 
is illustrated in Figure 7-8. The third orthogonal cross-section would typically be a 
plot of the amplitude of one particular mode versus both time and array tilt. 
Instead, a somewhat more useful output is presented here to convey information in 
the third dimension. That output ts the beamformer's residual error versus time 
and assumed tilt angle. This parameter is more appropriate because it may be 
considered to represent a projection or an Integration of the amplitude data from all 
the modes into a single value which can then be plotted on the time-tilt surface. It 
is this mtegration effect across all the modes included in the beamformer that allows 
the residual error plot to provide direct information concerning the tilt of the array 


orsus tini. 


7.3.1 The Sensitivity of Vertical Array Data to Array Tilt 


The sensitivity of the modal decomposition. process to vertical array. tilt. has 
been demonstrated theoretically and bv simulation in Sections 6.2 and 6.3. Similar 
sensitivity 1s experienced when dealing with actual data. as seen in Figures 7-9 and 
7-10.) Figure 7-9 shows a plot of the modal amplitude outputs from the multiple 


beam least squares beamformer plotted versus assumed array hs between -10 and 








Figure 7-9: Output of Multiple Beam Least Squares Modal Beamformer 
versus Array Tilt for 17.75 Hz Tone (Vertical Array) 
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Figure 7-10: Output of Multiple Beam Least Squares Modal Beamformer 
versus Array Tilt for 47.00 Hz Tone (Vertical Array) 
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10°. These results are computed from one particular 120 sec window of data 
centered on the 18,710 sec mark In the signal time series (the amplitude and phase 
time series of this signal are displayed for the different vertical array hydrophones 
in Figures 3-11 and 3-12). The estimated modal amplitude levels for all three 
resolvable modes are included on the plot. Even over just the range from -5 to 5 
(the most likely region for the actual effective tilt angle). significant. variations in 


the modal amplitude estimates can be found. These variations are as much as I dB 


for mode 1 and 10 dB for mode 3. 


Figure 7-10 shows similar results for the 47 Hz tone from the vertical array 
data set; here, a total of seven modes can be resolved. The 2 minute window of 
data used for this figure is centered on the 1760 sec point of the signal time sertes 
(the time series for this signal can be found in Figures 3-9 and 3-10). The problem 
of tilt angle sensitivity obviously increases at higher frequencies; at any particular 
frequency, the estimates of higher order mode amplitudes appear to vary more with 
array tilt than do those of lower modes. Note the deep nulls displaved by some of 
the modes at selected tilt angles; variations in amplitude of more than 15 dB can be 
seen over tilt variations of Just 1. in the 17 Hz result, for example. It Is apparent 
that if one wishes to obtain accurate measurements of the various modal 
amplitudes, then estimates of the array tilt accurate to about +0.25 must be 
avallable. At lower frequencies and lower mode numbers, this accuracy. requirement 
may be relaxed somewhat. Ht is also important to observe that the sensitivity of 
amplitude measurements to assumed array tilt errors is quite dependent on both the 
mode being considered and the actua: tilt value; a 0.1. change in array tilt is much 


more significant for some mode number/tilt angle combinations than for others. 


Figures 7-9 and 7-10 do provide qualitative information about the various 


modal amplitudes, even without determining the actual arras tdt. This information 
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can be extracted through the gross averaging of results over all possible array tilts. 
The 17.75 Hz results indicate a strong first mode arrival relative to the second and 
third mode amplitudes; this is entirely consistent with previous conclusions for data 
at this frequency made in Sections 3.12 and 7.2. Similarly, the 47 Hz outputs of 
Figure 7-10 are in agreement with previous observations of a very low level 
amplitude for the first mode. This consistency generates confidence in both the 


operation of the beamformer and the validity of the observations. 


7.3.2 Determination of Vertical Array Tilt Angles 


Having demonstrated the sensitivity of the experimental results to the tilt. of 
the vertical array, the next major step in processing the data is to make sufficiently 
accurate estimates of the array tilt at various times. If rn situ measurements of 
both the magnitude of the tilt angle and its direction in azimuth were available. 
then they could be employed to resolve the issue directly. Since no such 
measurements are available for the FRAM [IV Experiment, the alternative procedure 
outlined in Section 6.6 is employed here to make an estimate of the array tilt for 
each signal in the vertical array data set. This procedure involves minimization of 
the beamformer's residual error over the range of physically reasonable array tilts. 
Recall that the multiple beam least squares modal beamformer computes the full set 
of complex modal amplitudes which best fit the observed sound pressure field, given 
one assumed effective tilt angle. The normalized residual error 1s the performance 
measure for this fit (1.e.. a measure of how good the fit really is). Even though thre 
residual error is a scalar value, it should be remembered that its calculation is 
affected by the number of modes included im the multiple beam algorithm. By 
selecting. the. assumed. tilt. angle to. minimize the residual error, the outputs 
generated by the beamformer then form the set of amplitudes that best fit the 


observed field for any possible array tilt. 
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Figure 7-11: Residual Error of Multiple Beam Least Squares Modal 
Beamformer versus Array Tilt for 17.75 Hz Tone 
(Vertical Array with 3 Modes Included) 








с о о © © о O © 3 КӘ 
о о O о O O O O O O 

E e 

- Ci го) + Те о ~ 00 c — 

c | | | | | | | | | 

+ + ж ж ж ж ж ж ж ж 

END 10.00 DG 
00 00 4 00 2 00 


ТІК 


CONT 


-230- 

: (^ FR on (3 ` T о о O 
= С С C Sd сз <=) e =) 5 
= N г“) ч 42 © ГЕ О о СЭ 
— -- =— e— e e- — co — ON 
| | | | | | | | | | 
их ж. ж ж 
START: -10 00 06 

0 00 2 00 00 5 00 00 


10 00 DG 
1 00 DG SPAN 20 00 
1 00 DB МАХ 0 00 


MODAL BEAMFORMER RESIDUAL FRROR AT 


17 





START -10 00 DG 
DG CHAN | 
ов MIN -20 00 DB 
75 HZ {3 MODFS) 


Figure 7-11: Residual Error of Multiple Beam Least Squares Modal 
Per а е Pilt for 17.7» Hz Pone 
(Vertical Array with 3 Modes Included) 


оа саан Е Ете г 


= 
о 


20000 00 


TM 


19000 00 


TM 


18000 00 


TM 


17000 00 


TM 


Ем 20560 00 


16659 99 






ші Майн | 
іти тИ 
| ПИТИ AMINY 


) 20 


+ 


-230- 
С ^ (X 47 y y O 
су Ce O O (= о о сз mM =) 722) C C) C (y O O O 
O O O O сз о O O O 
=) == CN 122) E 42 (O r a 
= C1 T» ч Шо о ~ 02 Cn = = = теш = uc = m s 
| | | | | | | | | | | | | | | | | | 
к + x X ж ж ж ЖК X X X X X X X X ж X 
END 10. 00 DG START -10 00 
00 8 00 00 2 00 00 2 00 00 00 
4 
| 
00 
END 10 00 DG START -10 00 
TIK 1 00 DG SPAN 20 00 DG CHAN ! 
CONT ı 00 pa MAX 0 00 06 MIN -20 00 


MODAL BEAMFOPMER RESIDUAL ERROR AT 17 75 HZ (3 MODES) 


И oe O 





O 

= 

O 

ON 

| 
DG 
10 

о 
= 
о 
o 
о 
о 
еч 
= 
> 
о 
о 
о 
© 
о 
о 
= 
= 
о 
o 
o 
e 
о 
ао 
= 
h 
с 
о 
о 
= 
о 
m 





ER 


DG 


D8 


Figure 7-11: Residual Error of Multiple Beam Least Squares Modal 


Beamformer versus Array Tilt for 17.75 Hz Tone 
(Vertical Array with 3 Modes Included) 


Ем 20560 00 


i6659 99 





-231- 





| 
| 


Figure 7-12: Residual Error of Multiple Beam Least Squares Modal 
Beamformer versus Array Tilt for 47.00 Hz Tone 
(Vertical Array with 7 Modes Included) 
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Figure 7-12: Residual Error of Multiple Beam Least Squares Modal 
Beamformer versus Array Tilt for 47.00 Hz Tone 
(Vertical Array with 7 Modes Included) 
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Indications exist that the array may not necessarily be completely stationary 
over periods of longer than about twenty minutes; the discussions in Sections 3.12 
and 7.3.4 are germane. For the purposes of determining array tilt, however, any 
possible dynamic shifts in array shape over the signal’s duration have been ignored, 
so that one average value of the array tilt can be computed for the full length of 
each signal. This is done primarily for ease in implementing the software necessary 
to perform the processing. А secondary consideration is that the tilt angle 
estimation scheme appears to lack the resolution necessary to accurately track much 


more than the largest variations. 


Figure 7-11 shows a plot of the normalized residual error of the multiple beam 
least squares beamformer on the time-tilt surface for the 17.75 Hz data. Here. as in 
Figure 7-9, the three modes resolvable by the array at 17.75 Hz have been included 
in the multiple beam algorithm. The error values have been computed at every 
0.1° of array tilt in the range from -10 to 10°. The array tilt generating the 
minimum residual error is -0.15°. Data taken when the signal was off, as well as 
the data near the 17,100 see point (When the airgun was active), have been excluded 
to avoid biasing the time average. The tilt that provides the minimum error is 
reasonably consistent over the length of the signal, and errors as small as -15 JB 
(fitting 97° ofthe data) are achieved. The broad nature of the minimum presents 
some difficulty. While a determination of the minimum point can be made to any 
desired accuracy, some question exists as to the physical relevance of that accuracy. 
From the nature of Figure 7-11, it appears that the minimum tilt angle is probably 
known to better than +1 , but that the 0.25 — accuracy desired has not been 


achieved. 


Hieure 7-12 illustrates a similar set of results for the 17 Hz tone. Tle residual 


error computation now Includes the seven modes resolvable at this higher Frequency. 
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No value of array tilt in the expected range does a very good job of fitting the 
observed data; even the best fits account for less than 75°%o of the signal energy 
(equivalent to a -6 dB error level). This result is believed to be due to the presence 
of significant amounts of signal energy in modes bevond the seven included in the 
beamformer. In the horizontal array data of Figure 7-6, such a distribution of 
energy in the higher modes has already been demonstrated; the preliminary analysis 
of Section 3.12 supports the same view. Further, in Section 6.6, it has been shown 
that such a distribution would effectively disrupt the ability of the scheme to 
properly estimate array tilts. Thus, the results of Figure 7-12 should be interpreted 
as proving that the first seven modes form an inadequate model for explaining the 


structure of the received field. 


The two error plots just presented are typical of those generated by the 
remaining signals in the data set. In general, the received fields for signals below 
frequencies of about 40 Hz are adequately modeled by the modes that cau be 
resolved with the FRAM IV vertical array; therefore, reasonable estimates of array 
tilt can be generated from tlie residual error analysis for these tones. Minunuin 
residual errors of less than -10 dB (better than a 90f( fit) are typically achieved for 
these signals. The minimum can normally be determined to an accuracy of better 
than +1 for these data sets. but the desired 0.25 accuracy appears unachievable 
with this method. Based on Table 5-H, the technique probably provides accurate 
first mode results and partially accurate second mode results. The reliability. of 
second aud third mode amplitude estimates cam be expected to improve somewhat 
at lower frequencies for two reasons: the sensitivity of the estimates to tilt angle 
decreases; and the residual error plots for these frequencies tend to provide better 


mise tilt «aes 


Dy. contrast, all of the signals above 10 Hz in frequenes provide results similar 
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to those shown in Figure 7-12, making even a rough estimate of array tilt 
impossible. Additionally, analysis of the 35.50 Hz signal, while producing a 
respectable minimum error, generates a double minimum several degrees apart. 


Therefore, no reliable tilt angle estimate can be made for this tone either. 


The array tilts measured by the residual error minimization technique and the 
corresponding minimum errors achieved are listed in Table 7-I. The best available 
fit generally becomes poorer with increasing frequency. The excellent fits achieved 
at the lowest end of the frequency spectrum must be at least partially attributed to 
the range dependence of the channel, which tends to suppress propagation of many 
of the modes that cannot be resolved by the array at the lower frequencies. Note 
that the tilt angles are reported to the nearest 0.05° only because the beamforming 


algorithm requires an input of infinite precision. 


Arras lt 





Minimum Corresponding Assumed in 
Frequency Error Array Tilt Processing 
15.00 Hz -17 dB ШЕШ 3.05 
liane: Tz -14 dB -0.15 -0.15 
20.00 Hz -12 dB 2.95 2.05 
23.50 Hz -10 dB 2.20, 2.30 
30.00 Hz -11 dB 12105 2702 
32.25 Н? -12 dB -1.65 - 
35.50 Hz -9 dB - - 
47.00 Hz -4 dB - - 
ООЗ -4 dB - - 
55.00 Hz -5 dB - - 
71.00 Hz -7 dB - - 


Table 7-I: Results of the Array Tilt Estimation Procedure 


As mentioned above, the residual error analysis does not provide the til angle 
precisely enough to accurately estimate all of the modes that are resolvable by the 


A mates has һе refined further by comparing and 
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averaging the angles estimated from tones broadcast consecutively. Such an effort 
also contributes by providing an independent means for checking the accuracy of 
the tilt angle estimation process. A review of Table 3-IV, for example, indicates 
that the 15, 20, and 30 Hz signals were all broadcast within one hour of each other. 
It is reasonable to expect that their array tilts should be approximately equal. This 
is in fact the case, with the largest difference between any pair of the three being 
only about 1.5°, implying an accuracy of roughly 0.8 ° in these three tilt angle 
measurements. Assuming that the array did not move significantly over the course 
of the hour during which the three tones were broadcast, the accuracy of the tilt 
estimate can be improved by averaging the three individual findings. This vields a 
tilt angle estimate of about 2°, which is the value that has been used in the 


subsequent processing for all three signals. 


The 17.75 Hz data was not taken concurrently with any other signal in the 
data set, so there is no way of verifying its estimated tilt angle. The 35.25 Hz data 
was taken just after the 23.50 Hz signal: a comparison between the array. tilt 
estimates for these two cases is not very favorable, though, since a difference of over 
4° is found. Both residual error plots show minimum error tilt angles that are 
reasonably stable over the length of the signal. Tt is hardly likely that the array 
demonstrated a constant tilt for the 55 minute duration of the earlier signal. and 
then moved 4 in tilt during the 5 minute pause between signals before once more 
becoming stationary for the full duration of the second signal. [t is more reasonable 
to suspect that at least one, and possiblv both, of the estimates are being biased by 
signal energy in the unresolvable modes. Such an effect has been demonstrated 
with svnthetic data in Figure 6-8. The horizontal array results of the previous 
section imply that the 35.25 dz result is probably the less reliable of the two 


estimates. Therefore, no tilt angle has been assigned for subsequent processing of 
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this tone. In consideration of the 35.25 Hz result, however, a small correction has 


also been made in the tilt angle assumed when processing the 23.5 Hz data. 


7.3.3 Time Series of the Modal Amplitudes 


Figure 7-13 displays the modal amplitudes as a function of time for the three 
modes resolvable by the array at 17.75 Hz. These levels were computed assuming a 
tilt angle of -0.15°, as indicated in Table 7-I. The results shown here are in 
agreement with the amplitudes shown in Figure 7-9 for this tilt angle. The 
amplitude of the first mode ts about 113 dB re 1 Pa; the second mode amplitude is 
7 dB below this. The relative amplitudes of the first and second modes thus agree 
with predicted values, assuming that both modes are attenuated equally. Тіс 
amplitude of the third mode when compared with the first two is somewhat more 
problematic. The measured amplitude is 5 dB lower than that of the second mode. 
while approximately equal levels are expected. This error can best be explained by 
referring back to Figure 7-9. At an assumed tilt angle of -0.15 . the third niode 
estimate lies In a region of high sensitivity to tilt angle, so that even small changes 
in the assumed array tilt can greatly alter the beamformer output. These results 
support the earlier conclusion that at these frequencies, the array tilt is known well 


enough to accurately determine at least the first two modes, 


Only about 1 dB oof fluctuation is seen in Figure 7-13 for both the first and 
second mode estimates over the 55 minute duration of the signal (not meluding the 
period of the airgun blast). These levels are in agreement with those seen In the 
horizontal array data. The seemingly large Е у ама поп over tine are due to 
the fact that onlv three data points are represented along the vertical asis. 
"luetnations in the third mode amplitude are onh a little larger. perhaps 2 AB. 


Thus observation provides a strong Indication that the array was essential 
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Figure 7-13: Output of Multiple Beam Least Squares Beamformer 
versus Time for 17.75 Hz Tone (Vertical Array) 
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stationary for the duration of the signal, since even small amounts of array 
movement would generate large fluctuations in the highly sensitive third mode 


estimate. 


For the signals with frequencies of 30 Hz and below, reasonably reliable array 
tilt estimates have been made; hence, reliable time series of the modal amplitudes 
ean be obtained. All the time series show a dominant first mode structure, which is 
in general agreement with the horizontal array results. The fluctuation levels 
observed in these time series are also often comparable to the levels seen in the 
horizontal array data. In certain instances, more significant fluctuation levels are 
observed. The time scales associated with these larger fhictuations are almost 
always larger tlran 10 or 20 minutes. Since it is hypothesized that the cause of these 


variations Is vertical array movement, they are taken up in the following section. 


Since the array tilt cannot be determined for the vertical array signals with 
frequencies greater than 30 Hz, valid time series results cannot be generated. The 
best that can be done is to assign a number of different tilt angles on a provisional 
basis and then to examine the results, realizing that the absolute modal amplitudes 
are inaccurate, but possibly allowing. some general conclusions to be drawn. 
particularly about signal fluctuations over time. The modal time series generated 
by this procedure show a pattern similar to that seen at lower frequencies. Some 
signals show only very little fluctuation, and when larger variation is observed over 
time, it i5 invariably associated with the time scales of longer than 10 minutes. 


These also are discussed in the following section. 


7.3.4 Possible Array Movements Seen in the Vertical Array Data 


пр пе аге юм lave tidications of possible array 


movement over the duration of the signal. One such example is the 23.5 Hz tone. 
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Figure 7-14: Residual Error of Multiple Beam Least Squares Modal 
Beamformer versus Array Tilt for 23.50 Hz Tone 
(Vertical Array with 4 Modes Included) 
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Figure 7-15: Output of Multiple Beam Least Squares Beamformer 
versus Time for 23.50 Hz Tone (Vertical Array) 
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The residual error plot for this signal is presented in Figure 7-14, while Figure 7-15 
provides the modal amplitude time series data for the signal for an assumed array 
tilt of 2.30°. The amplitude time series show a small (but consistent} trend over 
time; the first and second mode amplitudes tend to increase slightly, while the third 
mode amplitude decreases. This drift ts also reflected by an equivalent trend in the 
residual error data, which shows about a 1^ change in the minimum error tilt angle 
over the hour duration of the signal. These variations may be due to changes nir tlie 
array tilt over time, or they may be actual fluctuations caused by chauges in the 
acoustic transmission path. Unless actual array tilt measurements are available. the 
two effects are fundamentally inseparable. Certainly, the data presented here does 
not preclude either explanation. However, the consistency of the residual error 
result over time suggests that at least some of the variations observed here are due 
to array movement during the signal. It is important to observe that the variatious 
under discussion appear to have time scales in the 10 to 20 minute range, and that 
there is no evidence of significant fluctuations possessing significantly shorter time 


scales. 


Although it is not possible to obtain accurate modal amplitude results in an 
absolute sense for the signals that are above 30 Hz in frequency, the stability of 
these signals over time can be assessed in a qualitative fashion by computing the 
modal decomposition for an arbitrary tilt angle. For such a situation, the outputs of 
the beamformer can be expected to be constant over time if no significant temporal 
fluctuations are present in the sound field. When this is done, many of the signals 
above 30 Hz demonstrate a considerable amount of long term variation in the 
resultant modal amplitudes. И is even more difficult to assess the meaning of these 
long term fluctuations than it is for those present it the lower frequency data. since 


meaningful residual error analvsis ls not available. 





-242- 

In summarizing the issues of fluctuations and array movement im the modal 
amplitude time series generated by the vertical array, it can be said that the short 
term amplitude fluctuations observed (those with time scales on the order of 
minutes) are small, and roughly equivalent in size to those seen in the horizontal 
array data. Some longer term fluctuations, of time scales that are greater than 10 
or 20 minutes, have been observed in the vertical array data. .\ definitive 
separation of these varlations Into array Movement effects and propagation path 
fluctuation effects ıs not possible, due to the lack of field measurements of vertical 
array tilt angle over time. Two observations provide circumstantial evidence that 
at least the fluctuations in the lower frequency data (that at 30 Hz and below) may 
be primarily related to array movement. The first Is the nature of the residual 
error analyses available for the low frequency data: these analyses are consistent 
with the array movement explanation. The other is the lack of any similar 
fluctuations in the processed outputs of the horizontal array data, which indicate 
that the phenomena are probably array specific. Long term fluctuations In vertical 
array modal amplitudes for frequencies greater than 30 Hz are similar to those 
below 30 Hz; reliable residual error analyses and comparative horizontal array data 


are lacking for these signals, though. 


7.4 Modal Amplitude Analysis 


The absolute modal amplitude estimates available from the beamformers can 
be exploited in two ways. The values generated from the vertical array data by 
modal beamforming can be compared to those made with more standard plane wave 
techniques employed with the horizontal array data in order to quantitatively verify 
the consistency of the two methods. This must be accomplished with first mode 


data, since the horizontal array cannot provide accurate amplitude estimates for 
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higher order modes. Also, the comparison of observed levels to those predicted from 
propagation theory provides insight into the unique nature of long range 
propagation of low frequency sound in the central Arctic, and allows experimental 


measurements of modal attenuation coefficients to be made. 


7.4.1 Comparison of First Mode Amplitude Estimates from the 


Horizontal and Vertical Arrays 


Figure 7-16 summarizes the estimates of the first mode amplitudes versus 
frequency for both the horizontal and vertical array data sets. When the two 
groups of results are compared, it can be seen that the horizontal аггау values 
uniformly exceed the vertical array levels by about t dB over the 15 to 30 Hz range. 
the region in which both sets of estimates are most accurate. The lack of accurate 
array tilt estimates make the vertical array results outside this range less reliable: 
the values shown for the horizontal array at 47 and 53.25 Hz are also questionable. 
since the first mode no longer dominates in these signals. For this reason, the data 
points at these frequencies have been connected with a dotted line. The £ dB 
difference between the two sets of results appears to have two causes. A slight 
mismatch between horizontal and vertical array hydrophone sensitivities accounts 
for about 1 or 2 dB of the difference. The remaining 2 dB error is due to the 
shehtly different propagation paths over which the two sets of signals were 
transmitted. Although the signals in each data set Were grouped fairly closely in 
time, eight days elapsed between the last horizontal array data and the first vertical 
array data studied here. Note, however, that the issue is not merely one of range, 
as the PRISTEN/FRAM range for the horizontal array data (265 ki) was actually 


acia top tiie lower level vertical array data (25! kin). 








First Mode Amplitude (dB re 1 uPa rms) 
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Figure 7-16: Summary of First Mode Amplitude Measurements 
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Figure 7-18: Summary of Third Mode Amplitude Measurements 








7.4.2 Modal Attenuation Coefficient Measurements 


In addition to the data described above, Figure 7-16 also shows two other 
types of information. The amplitudes of the first mode estimated from vertical 
array data at higher frequencies are presented in a different fashion than those at 
lower frequencies. Since no accurate estimate of array tilt is available at. these 
frequencies, the range of amplitudes generated by all possible array tilts between 
-4° and 4° is displayed instead. The actual markers that are connected by the 
dotted line indicate the first mode levels for one subjective estimate of the actual 
array tilt. While this data is much less accurate than the lower frequency results, 
the ranges do provide valuable information about the general trend of the amplitude 


with frequency. 


Also presented is the amplitude of the first mode as a function of frequency at 
a range of 254 km for an ideal horizontally stratified channel (ie. geometric 
spreading with no attenuation mechanisms present). Figure 3-3 has been used as 
the source spectrum in computing these levels. The difference between this curve 
and the experimental results indicates the attenuation loss for the first mode as a 
function of frequency. It is obvious that an attenuation mechanism Is at work ìn 
the sound channel. and that it gets quite severe at higher frequencies. The 
attenuation is the cause of the unexpectedly low first mode amplitudes seen above 


40 Hz in the data. 


Figures 7-17 and 7-1& summarize similar results for the second and third 
modes, respectively. Again, the effect of the attenuation mechanism can be seen, 
particularly at the higher frequencies. Note that the measured modal amplitudes 
for the higher modes display greater scatter than do those for the first mode: this is 
indicative of the increased sensitivity to array tilt errors that the modal 


beamforming process exhibits at higher order modes. 
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Figure 7-20: Second Mode Attenuation - Experimental and Predicted 
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Figure 7-21: Third Mode Attenuation - Experimental and Predicted 





То characterize the modal attenuation further, Figures 7-19, 7-20, and 7-21 
have been included. In these plots, the actual attenuation experienced by the first. 
second, and third modes is displayed as a function of frequency. Vertical array 
results at higher frequencies are again plotted as ranges connected by a dotted line. 
The data is plotted on a scale of total attenuation at 254 km; equivalent attenuation 


coefficients can be recovered simply by dividing the dB levels by this range. 


Experience [25] suggests that the cause of the attenuation is the rough ice 
canopy; a comparison of the three plots supports this supposition by demonstrating 
directly that the attenuation mechanism is located in the surface duct of the central 
Arctic profile. Such a conclusion can be reached by noting the relationship between 
the onset of severe attenuation and the mode number. In Figure 7-19, the 
experimentally determined attenuation for the first mode is seen to be strongly 
increasing, even at the lowest frequencies in the data set. By comparison, the point 
where attenuation becomes significant is delayed somewhat in frequency for the 
second mode (Figure 7-20), and even more for the third mode (Figure 7-21). 
Typical frequencies for the onset of severe attenuation are 30 Hz for the second 
mode and 50 Hz for the third mode. These values correlate well with the 
frequencies at which the various modes transition into the surface duct, as discussed 
in Section 4.5.2. From this, one can conclude that the attenuation mechanism must 


be located in the surface duct, making the rough ice surface a prime candidate. 


Figures 7-19 through 7-21 also include plots of predicted modal attenuation 
versus frequency for two different rough surface scattering theories. Both cases 
assume a rough pressure release boundary at the ice-water interface. The statistics 
of this boundary are assumed to be those given by DiNapoli and Mellen [25], as 
discussed in Section 4,3. Тһе simpler case is the well known Kirchhoff 


th 


approximation, for which the attenuation coefficient: of the 77 mode can ђе 


2 





computed from 


(7.1) 





Here ø is the surface roughness (taken to be 2 m) and c¿(0) is the sound speed at the 
surface, while 2(0) is the slope of the mode’s normalized shape function evaluated at 
the surface and c, is its horizontal phase speed. The other case is the full method of 
small perturbations (MSP) solution, as developed by Kuperman and Ingenito [17]. 


In this development, the modal attenuation 1s shown to be 
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In equation (7.3), PX) is the two-dimensional spectrum of the rough surface (here 
given by equation (4.1) with an assumed correlation length of 44.8m), and wis a unit 
vector in the direction of propagation. The integration Is carried out over a full 
range of vectors 5 in the horizontal plane. Details of the derivations of both these 
results can be found in [47]. It is important to note that both contain factors that 
include the square of the slope of the mode shape evaluated at the surface, a term 
which is not present tn the more familar plane Wave scattering results, and, hence. 
often overlooked. These factors reflect the ability of the ice to scatter modes with 
their energy concentrated near the surface more strongly than those with their 


energy concentrated at deeper depths. 


The attenuation. predictions of both models for the first three modes have 


been plotted in Figures 7-19 through 7-21, respectively. The agreement of the 
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experimental data with the Kirchhoff approximation is uniformly poor. On the 
other hand, the attenuation results predicted by the MSP match the first mode data 
quite well, and the second mode data somewhat less exactly. It is possible to 
speculate that the poorer second mode fit is caused by the scattering of first mode 
energy into the second mode by the ice. In general, anv scattered energy 
redistributes itself across the modal spectrum. If significant amounts of motion were 
present at the surface, the scattered energy could be expected to be temporally 
incoherent when compared with the transmitted signal. For the Arctic. though, it 
appears reasonable to model the rough surface as being essentially stationary. over 
time, so that an scattered energy would still be coherent with the specular 
transmission. Much of the scattered energy is lost into higher order modes that do 
not propagate well in the channel; however, a considerable percentage can be 
expected to end up in the lower order modes. In fact, the second mode will 
normally receive the largest share, since modal coupling is typically strongest 
between neighboring modes. In the 30 to 60 Hz frequency band, the first mode 
excitation exceeds that of other modes by 15 to 30dB (see Figure 4-9). Thus. energy 
scattered from the first mode into the second mode may be significant when 
compared to second mode excitation levels, even though the coupling coefficient 
Itself is negligibly small. Under such circumstances, the attenuation experienced by 


the second mode would appear to be reduced, as ts the case in Figure 7-20. 


The MSP also appears to fit the third mode data (Figure 7-21) reasonably 
well, although there are hints of a phenomenon similar to that of the second mode 
in {һе 71 12 data point as well. It must be remembered that the third mode data mr 
probably only partially resolved at best, due to the inability to estimate the array 
tilt angle to sufficrent accuracy. As might be expected, the third mode data display 


less consistency than those of the other two modes. 
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Figure 7-22: Sonogram of Explosive Shot Transmitted Over the 
TRISTEN/FRAM Channel as Received on 90 m Hydrophone 
of the Vertical Array (300 ft Source Depth) 
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7.4.3 Comparison with Shot Data 


It is of interest to contrast the horizontal and vertical array data results Just 
presented wıth data taken from an impulsive source broadcast over the same 
channel. This exercise confirms the horizontal array modal received level 
measurements directly and the vertical array results indirectly. It also generates 


insight into the relationship between the results of two data sets. 


Figure 7-22 displays the sonogram of a shot set off at 300 feet and then 
transmitted over the TRISTEN/FRAM sound channel. The sonogram is a plot of 
received level versus frequency and time, and is generated by computing a sequence 
of very short length spectral estimates from the time series output from a single 
hydrophone, in this case the 90 m hydrophone from the vertical array. The results 
shown here were generated using a Burg spectral estimation technique to provide 
greater frequency resolution in the 250 msec window length used. A new spectral 
estimate was completed every 100 msec. The original time series for this shot has 


already been presented in Figure 1-2. 


The peak marked with a dotted line is the first mode arrival, corresponding to 
the long tail of the signal displayed in Figure 1-2. This can be shown from a 
comparison of the measured dispersion characteristic with that expected of the first 
mode. Note that the received level of the first mode arrival is quite strong and ds 
reasonably flat versus frequency up to about 40 Hz, after which it almost completely 
disappears. This 15 exactly the effect that has been demonstrated in the horizontal 
array tonal data which, in turn, have been shown to be in quantitative agreement 
with the vertical array tonal results in Figure 7-16. In retrospect, the strong drop in 
the received level of the first mode in both the horizontal array data and the shot 
data can be seen to be a combination of two effects. The first mode exhibits some 


attenuation even at the lowest frequencies in the data set, and this attenuation 
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increases rapidly at higher frequencies. Below 40 Hz, this attenuation is offset bv 
the increase in size of the mode shape at 91 m, which causes both the source and 
the receiving hydrophones to couple to the first mode with increasing efficiency. 
The two effects almost cancel themselves, making the first mode received. level 
appear relatively flat over this frequency range. Above 40 Hz the size of tlie first 
mode shape at 91 m stops growing with frequency (as is seen in Figures 4-8 and 
4-9). The balance between the two effects ıs lost, and the now severe first mode 
attenuation dominates the remaining results, causing the very quick drop in first 


mode received level that ıs seen in the data. 


7.5 Mode Coherence Measurements 


The results of the previous section indicate that the modal beamformer ts 
performing reasonably well for frequencies of 30 Hz and below, where reliable array 
tilt estimates are available. This means that meaningful estimates of the coherence 
between the various modes can also be made. The results of these coherence 
measurements for the vertical array signals below 30 Hz are given in Tabie 7-1. 
The window of data extracted from each signal to make the coherence estimate ts 
the longest uninterrupted period available; the lone exception to this is the 23.5 Hz 
data, where the processing software limits the number of data points used. Air gun 
blasts, such as the one imbedded in the 17.75 Hz data, and other anomalies have 


been avoided so that the results would not be biased by invalid data. 


The resultant coherence estimates are striking because of their uniformly large 
values. For all intents and purposes, these results show that the central Arctic 
channel can be considered to be a completely deterministic medium mp time, and 
that all the different source-receiver propagation paths are completely plrase locked. 


The slight variation that does exist between results appears to be strictly am issue of 











Frequency (Hz) 15.00 17.75 20.00 23.50 30.00 
Length of Average (sec) 400 2725 1120 2299 550 

Mode 1-2 Pair 0.99 0.98 0.99 0.99 0.99 
Mode 1-3 Pair 0.99 0.93 0.99 0.99 0.99 
Mode 1-4 Pair - - - 0.98 0.99 
Mode 1-5 Pair - - - - 0.99 
Mode 2-3 Pair 0.99 0.93 0.99 0.99 0.99 
Mode 2-4 Pair - - - 0.98 0.99 
Mode 2-5 Pair - - - - 0.99 
Mode 3-4 Pair > - - 0.97 0.99 
Mode 3-5 Pair - - - - 0.99 
Mode 4-5 Pair - - - - 0.99 


Table 7-1: Mode Coherence Estimates for the Vertical Array Data 


SNR. The low levels of coherence in all pairings that include the third mode at 
775 Hz, for example, can be explained by realizing that the third mode amplitude 
is far enough below that of the first mode to allow the background noise level to 
exert a small amount of influence in the coherence estimate. Since confidence 
intervals for coherence estimates very close to unity are also «uite tight [13]. the 


estimates are guaranteed to be quite reliable from a statistical viewpoint. 


While true coherence measurements cannot be made for vertical array signals 
above 30 Hz, gross averages of observations over many different assumed array tilts 
are agam possible. Like the lower frequency results, the modes at the higher 
frequencies also appear to be almost completely coherent. Modes that appear 
Incoherent for a given tilt angle can normally be explamed by an SNR argument. 
On the average, some small amount of coherence degradation can be seen at the 
higher frequencies, indicating that upper mits in both frequency and averaging 
length exist, beyond which the present results are not valid. Typical coherence 


СО ООО Гог Ще № Hz data and roughly 0.92 for the 71 НА data, 
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compared to values of 0.99 seen at lower frequencies. None the less, these coherence 
levels are close enough to unity to conclude that the present experiment only covers 


a fraction of the frequency/duration space over which these conclusions are 
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Figure 7-23: Time Series of the Relative Phases of Various 
Mode Pairs for 17.75 Hz Data 


7.6 Relative Phase Measurements for the Vertical Array Data 


The high levels of coherence just demonstrated between the various modes 
indicate that their relative phases should be very stable over time. This is the case, 
as is seen in Figure 7-23. The three time series displaced represent the variations of 
the relative phases with time for the three possible pairings of the modes resolvable 
at 17.75 Hz. As expected, the relative phases are seen to be extremely stable over 


the full 55 minute duration of the signal. The mode 1-2 pair is the most stable. 
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demonstrating not more than a few degrees of fluctuation over the course of the 
signal. The somewhat larger fluctuations in the pairs involving the third mode 
reflect the same SNR issue discussed in the last section. Since the value of the 


mode 1-3 pair is very near to 180°, a large amount of phase wrap is visible. 


Table 7-III gives the experimentally measured relative phases for all possible 
pairs of resolvable modes at frequencies below 30 Hz. The signal durations used to 
compute these values are identical to those used in the coherence estimates. 
Assuming a source-receiver range of 254 km and a range independent channel. à 
comparison of the phases measured at 17.75.Hz with the values expected can be 
made. The results are presented in Table 7-IV. The significant differences between 
the theoretical and experimental data indicate two things. First, the calculation of 
the expected relative phase is sensitive to inaccuracies in the assumed horizontal 
phase speeds of the modal pairs. For example, a 1 m/sec difference in the first 
mode phase speed results in about a 40° in the predicted data of Table 7-IV. 
Second, the measurement is also quite sensitive to range variations in the channel, 
since the total phase represents an integration of all the КЕНЕС encountered 
in propagating through the channel. These appear to be non-negligible for this 
situation. Thus, while the relative phases themselves are quite stable and well 
behaved, there is some question as to whether it ıs possible to predict them 


accurately enough to accomplish passive source ranging in this fashion. 


7.7 Summary 


onal and vertical array (onal data sets from the FRAM IN Arctic 
Acoustie Experiment that are described in Sections 3.9 and 3.10 have been 
processed using the algorithms discussed in Chapter 5. Fhe results for both have 


been presented, compared, and analyzed in a variety of ways. 
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Frequency (Hz) 15.00 HN. 20.00 23.50 30.00 
Mode 1-2 Pair 171^ 102° 120° 48° BS 
Mode 1-3 Pair -153° -166 ° -292 ° -65 ° 174 
Mode 1-4 Pair - - - -154 ^ 85^ 
Mode 1-5 Pair - - - - 38° 
Mode 2-3 Pair 36 92° -142* lilo. 39° 
Mode 2-4 Pair - - - 158 -50 ° 
Mode 2-5 Pair - - - - -97 ° 
Mode 3-4 Pair - - - - 59 -89 
Mode 3-5 Pair - - - - -136 ° 
Mode 4-5 Pair - - - - -47 


Table 7-HI: Relative Phase Estimates for the Vertical Array Data 


Measured Expected 
Mode 1-2 Paır 102° о 
Mode 1-3 Pair -166 ° -16 
Моде 2-3 Рап 92° -170° 


Table 7-IV: Comparison of Measured and Expected Relative Mode 
Phases for 17.75 Hz Data 


The horizontal array data have been beamformed in azimuth using the single 
beam MLM algorithm in order to find the correct azimuth angles for the TRISTEN 
source. This processing results in azimuth angle estimates that are accurate to 
about +0.1°. The azimuth angles obtained agree completely with less precise 
estimates available from the navigation data. The MLM algorithm performs quite 
well in this application, producing a much narrower main lobe and significantly 


lower side lobe levels than an equivalent conventional beamformer. 


The horizontal array has also been steered in the vertical direction at the 
azimuth angles previously determined to study the modal structure of the received 
field. “Flis has been accomplished by varying the assumed horizontal phase speed 


with which the steering vector Is computed. The processing has been accomplished 
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twice: once with the single beam MLM algorithm to provide high resolution in pliase 
speed; and once with the conventional beamformer to recover the correct received 
levels. The horizontal array data below 30 Hz exhibit a strong peak at horizontal 
phase speeds closely matching those predicted for the first mode. Above 40 Hz, this 
first mode arrival vanishes; instead, a much more evenly distributed set of modal 
received levels is observed, possessing a very small peak at phase speeds 
corresponding to deep RSR path propagation. This occurs even though analysis of 
the modal structure suggests that near 40 Hz, the received level of the first mode 
should dominate those of other modes by as much as 60 dB. The region between 30 
and 40 Hz might be considered to be a transition region between the two 
propagation regimes. Although a strong first mode arrival is seen in this band. it is 


perhaps not as great as might be expected from studies of the modal structure. 


In the 27 Hz tone taken from the horizontal array, a strong fade similar to 
that seen by Mikhalevsky [57] has been identified. A study of the time series of this 
fade seems to indicate that it is a phase discontinuity which, in turn, causes the 
instantaneous frequency of the signal to walk out of the pass band of the low pass 
fitter. Analysis in both the horizontal and the vertical shows no observable 
differences between the structure of the received field before and after the fade, 


suggesting that the cause is probably source related. 


The vertical array data have been processed using the multiple beam least 
squares beamforming algorithm to produce modal amplitude and phase estimates. 
Three different types of output have been generated for the vertical array data set. 
One is a plot of the amplitude estimates for all the resolvable modes versus time for 
one given array tilt angle. The second is a plot of these same amplitude estimates, 
now displayed at one given time and across all reasonable values of array tilt. The 


third is a plot of the residual error versus time and assumed tilt angle. This last 
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may be considered to be a condensation or integration into a single number of the 
all the various modal amplitude estimates made at one time and one assumed tilt 
angle. The three outputs represent different cross-sections along mutually 


orthogonal planes of an amplitude function that depends on mode number, time, 


and assumed tilt angle. 


The modal decomposition of real data has been shown to be sensitive to the 
assumed array tilt in a fashion similar to the demonstrations made with synthetic 
data in Section 6.3. This sensitivity is seen to increase at higher frequencies and 
higher mode numbers. A tilt angle accuracy on the order of 40.25 appears to be 
needed if one wishes to properly estimate all of the modes resolvable by the vertical 
array at any of the frequencies of interest. This requirement can be relaxed 


somewhat for lower frequencies and lower mode numbers. 


As shown in Section 6.6, plots of residual error on the time-tilt plane allow 
estimation of the actual array tilt from the acoustic data. For the purposes of array 
tilt estimation, the residual error plots generated for the vertical array data have 
been used to compute one average tilt angle for the full duration of each signal, 
despite indications that the array may not be completely stationary over periods of 
longer than about 20 minutes. Where possible, the tilt angle estimates made bv this 
procedure have been verified by comparing estimates from consecutively broadcast 
signals. For frequencies of 30 Itz and below, it appears that the estimation 
technique is capable of generating estimates to an accuracy of about +0.7 , with 
somewhat better results at the very lowest frequencies in the data set. This 
precision allows valid first mode estimates and partially valid second mode estimates 
to be made; at the lowest frequencies, third mode estimates may also be reasonably 
trustworthy. All tilt angle estimates made by this technique fall inside the expected 


limit of 5 
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The residual error plots for signal frequencies above 30 Hz generally do not 
demonstrate the single deep minimum necessary to make accurate array tilt 
estimates. Even when estimates are possible, they do not appear to be consistent 
over consecutive signals. This problem is attributed to the increase at these 
frequencies of the proportion of energy carried by modes that are fundamentally 
unresolvable by the array. This effect has been demonstrated in the results of the 
horizontal array data. The outputs of vertical array modal decompositions versus 
both time and tilt angle also verify this effect. Above 30 Hz, the first mode 
amplitude also appears to be abnormally low when compared to those of higher 


modes in the vertical array data. 


The temporal stability of all the data examined in this chapter appears to be 
similar to that seen in the preprocessed time series of Section 3.12. In general, short 
term fluctuations for both horizontal and vertical array data are in the I dB range. 
Some larger short term variations have occasionally been noted on specific sensors 
in the horizontal array, but these do not appear to be coherent over the full array. 
Longer term variations (on the order of about 20 minutes) are seen in selected 
vertical array signals. The variations appear to increase as the frequency of the 
signal increases. No such long term fluctuations are observed in the horizontal 
array data. The lack of actual array tilt measurements makes the separation of 
channel fluctuation effects from array movement effects impossible, but nothing 1n 
the data precludes the latter from explaining the great majority of the long term 


Varlations seen. 


The first mode amplitudes, as computed from the vertical array data. have 
been shown to be consistent with those derived from the horizontal array data over 
the range of frequencies (below 30 Hz) where accurate estimates of each can be 


made. Both the horizontal and vertical array data indicate significant. attenuation 
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of first mode energy. Similar results are seen ın the second and third mode 
amplitude results taken from the vertical array data. The nature of the attenuation 
curves versus frequency for the three modes demonstrates that the source of the 
attenuation is located in the Arctic surface duct. A good fit to the measured first 
mode attenuation can be generated by applying the method of small perturbations 
to a rough pressure release surface located at the ice-water interface, where the 
statistics used for the rough surface have been determined experimentally from ice 
profiles. A simpler approach using the Kirchhoff approximation yields only a very 
poor match to the measured attenuation values. Second and third mode 
attenuation estimates made with the method of small perturbations match 
experimental data somewhat less perfectlv, predicting larger attenuation than is 
actually observed at higher frequencies. This result could possibly indicate that 
modal coupling between the first mode and neighboring modes cannot be ignored. 
If this is indeed the case, then the observed stability of the-received signals E 
the scattered field to be completely coherent with respect to the specular field, 


which, in turn, presumes à rough surface that is essentially stationary over time. 


Coherence measurements for the various mode pairs have been made for 
signals with frequencies 30 Hz and below, where reliable tilt angle estimates are 
available. Average modal coherence over periods of as long as 45 minutes is 
typically 0.99; all the measurements made exceed a value of 0.93. The few lower 
results observed appear to be caused by noise interference. By averaging results 
over a number of arbitrarily selected tilt angles, it is possible to get an Indication of 
the coherence levels hkely to be found at higher frequencies, Where actual array tilt 
estimates cannot be made. Typical values are 0.96 at 47 Hz and 0.92 at 71 lz. 
Therefore, the various propagating modes for any signal in this data set can be 


considered to be completely coherent over the full duration of the transmission. 
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The large coherence values observed also guarantee that they are highly reliable 
estimates. The slightly lower values seen at higher frequencies lead to the 
supposition that there is some frequency/duration limit beyond which these 


conclusions would no longer be valid. 


The high levels of coherence observed guarantee that estimates of the relative 
phase of various mode pairs are meaningful. These relative EM estimates appear 
to be extremelv stable, exhibiting only very small variations over time. However, 
the phases measured do not closely match those that would be expected if the 
channel were totally range independent. Thus, while the relative modal phases are 
quite stable, there appears to be some question about the ability to accurately 


predict them. 
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Chapter 8 


Summary and Conclusions 


This chapter is divided into three sections. The first summarizes the gencral 
topics covered throughout the thesis, while the second discusses the major 
conclusions, locating them in the larger framework of signal processing theory and 
underwater acoustics propagation theory. A short list of some ideas for further 


work along these lines completes the thesis. 


8.1 Summary 


The twin issues of modal beamforming and mode coherence have been 
Investigated at some length. The prime motivation for studying these topics comes 
from a desire to exploit the modal nature of sound propagation in a waveguide more 
fully, so that the source's range and depth mav be estimated directly from sound 
field measurements. To this end, two necessary contributions have been made. 
First, beamforming techniques that address the Waveguide nature of low frequency 
sound propagation in the oceans have been studied extensively and implemented on 
real data, including data taken from a vertical array; this has not been 
accomplished previously. Second, the relative amplitudes and phases of the various 
modes have been measured and their stability studied through the use of modal 


coherence estimates. This also is an original contribution of this thesis. 


Some of the practical differences between horizontal and vertical arrays have 
been analyzed. The most important difference from a signal processing view point ts 
the increased importance of vertical array sensor displacement, This is due to the 


increased size of the sensor displacements, Which is caused by the greater scope of 
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vertical arrays; and to the fact that the displacements occur in an organized pattern 


across the face of the array. 


A preprocessing scheme involving two stages of quadrature demodulation, PUR 
low pass filtering, and decimation has been implemented to reduce the data set to a 
usable size prior to further processing. This scheme results a complex output time 
series for each channel that is demodulated to within +1 mHz of the nominal signal 
frequency, filtered in a passband of 12.5 mHz, and downsampled to a sample 
interval of 1.6 sec. The resulting output time series contains about 100 independent 


degrees of freedom in a typical 55 minute signal. 


The modal structure of the sound channel encountered during the FRAM IV 
Experiment has been studied extensively, and. predictions of the mode shapes and 
associated horizontal phase speeds have been made for all frequencies of interest 
here. The results have been analyzed both with respect to relative mode shape sizes 
at 91 m (the deployment depth for the TRISTEN source and the FRAM horizontal 
array) and with respect to horizontal phase speed variations as a function. of 
frequency. The arrangement of the experiment makes both the source and the 
horizontal receiving array highly tuned to first mode propagation: first mode 
excitation typically exceeds that of other modes by more than 15 dB. The 
variations of phase speed with frequency allow cone to observe the transition of 
modes from the deeper sound speed profile into the strong Arctic surface diet. The 
effect of variations in the sound speed profile is shown to be small compared to the 


effect of the overhead ice canopy for central Arctic surface duct propagation, 


Four different algorithms for computing modal decomposition of the observed 
sound field have been developed theoretically. These four algorithms include single 
beam and multiple beam variants of both the least squares and. MEM . algorithms 


The least squares algorithm is shown to be a true. least. miean squares linear 





estimator for the complex modal amplitudes; its single beam variant is exactly the 
equivalent of the conventional beamformer. The MLM algorithm is shown to be a 
weighted least mean squares estimator, with the weighting essentially решо the 
inverse of the sensor noise. The performance of all four methods has been 


compared. 


In general, the multiple beam least squares technique can always be expected 
to provide better performance than its single beam counterpart, as long as the 
beams included together in the multiple beam approach are all different enough to 
avoid matrix singularity issues. Problems of this type are encountered when one 
begins to deal with modes not fundamentally resolvable to the array. The modal 
resolution of a vertical array may be hmited in two ways, by its length and by 
inadequate inter-sensor spacing. An overly short array makes it difficult to resolve 
neighboring modes, while inadequate iInter-sensor spacing causes Nyquist-like 
aliasing for modes widely separated in mode number. It has been demonstrated 
that the FRAM IV vertical array is array length limited for the low frequencies of 


interest in this thesis. 


The single beam MLM algorithm can be expected to provide resolution 
superior to that of any of the other approaches in situations where tlie assumptions 
underlying the method are met. The most important of these assumptions requires 
that all the energy in the various modes be completely incoherent. The presence on 
coherent energy in several different beams causes severe distortions of the single 
beam MLM algorithm output, though. This coherent interference problem can be 
theoretically eliminated by the use of the multiple beam variant of the MLM 
algorithm, but the practical realities ef this approach make it unattractive. Thus. 
the multiple beam least squares processor has been chosen for use with the vertical 


array data analyzed here. 1t. has been. demonstrated. that the single beam MEM 





- 269- 
algorithm still provides higher resolution than either of the least squares approaches 
when properly used with horizontal arrays, even though the coherent interference 
problem corrupts the output levels. Therefore, the horizontal array data have been 
processed twice, once with the single beam MLM algorithm for resolution and once 
with the single beam least squares algorithm for the purpose of quantitative 


estimation of the received levels. 


The effect of vertical array tilt on the modal decomposition process has been 
investigated. An examination of the fundamental physics of the process suggests 
that the modal beamformer results are very sensitive to variations m array tilt over 
the range of tilt angles most likely to be encountered with the FRAM IV vertical 
array, which ıs from -5 105 away from vertical. This examination also leads to 
a simple method for estimating such tilt angle sensitivity. Upon application to the 
FRAM IV vertical array, it is found that the tilt angle must be known to about 
+0.25° in order to properly resolve anything more than the first mode. This 
sensitivity to tilt angle increases at higher mode numbers and higher frequencies. 
Since no tilt angle measurements are available for the FRAM IV Experiment. a 
scheme for estimating the array tilt from the acoustic data has been developed and 
implemented in this thesis. The method involves the minimization of the гем па! 
error associated with the multiple beam least squares beamforming algorithm over a 
reasonable range of assumed array tilt angles. This scheme has been shown capable 
of tracking the array in tilt angle by simulation: its performance Is dependent upon 
the total amount of energy present im the modes that are included ш the 
beamformer relative to the amount present im modes not resolvable by the array. 
The approach is subject to bias when a large percentage of the signal energy that ts 


present resides in fundamentally unresolvable modes. 


Two extensive tonal data sets taken during the ТАМУ Хебе Acoustic 
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Experiment have been processed and analyzed in this thesis. One set was recorded 
from the hydrophones of the horizontal array available at the FRAM ice camp, 
while the other set was recorded from eighteen sensors of the vertical array 
deployed at the same camp. The horizontal array data have been analyzed to 
provide experimental amplitude and horizontal phase speed measurements for the 
first mode; these, in turn, have been used to successfully test the consistency of the 
modal beamformer output and the validity of the mode structure determined in 
Chapter 4. The vertical array data has been used to measure modal amplitudes and 
relative phases for various low order modes, and to make estimates of the temporal 
stability and. coherence of the complex modal amplitudes over time. The modal 
amplitudes measured have been analvzed with respect to the attenuation observed, 
and measured modal attenuation levels have been compared with. predicted. levels 
using both the Ixirehhoff approximation and the full method of small perturbations. 
Good agreement with experimental attenuation results has been obtained from the 
latter method using a simple rough free surface model and experimentally derived 
estimates of the ice canopy statistics. The stability and coherence measurements 
have been used to investigate the nature of the phase relationship between various 
modes. It has been concluded that the central Arctic sound channel can ђе 
considered to be completely deterministic in terms of signal propagation over 
periods of time well in excess of one hour, All propagating modes appear to be 
phase-locked, indicating that no independent, phase-random propagation paths exist 


ш the Arctic for these time scales and operating frequencies. 





8.2 Conclusions 


The important conclusions of this thesis fall into the same three general 
categories into which the contributions discussed in Section 1.4: results pertaining 
primarily to the discipline of signal processing; results dealing with the practical use 
of vertical arrays; and results affecting the field of underwater acoustics. Each 


category 15 discussed in turn. 


t 


8.2.1 Signal Processing Conclusions 


The most important signal processing conclusion Is that a method for 
decomposing the sound field generated by a distant tonal source and received at a 
vertical array into Its component normal modes has been proven to be effective by 
theorv, bv simulation, and by practical application to real data. This last as 
particularly significant, as no previous work dealing with narrowband modal 
beamforming has involved application to actual data taken in the field. Thus. tls 
thesis demonstrates the feasibility of estimating narrowband modal amplitudes and 


phases directly from data taken with a vertical array. 


A number of conclusions have been reached concerning the relative 
performance of the four beamforming algorithms investigated, primarily along the 
lines of comparisons between the least squares and MEM approaches and between 


single beam and multiple beam algorithms. 


The critical issue involved in any comparison of MLM and least squares 
beamforming algorithms for narrowband modal applications Is the coherence of the 
various modal amplitudes. What ts defined as the signal for the MILNE algorithin 
must be incoherent with respect to everything else in the sound field. which ts 
automatically defined as noise. Violation of this requirement has been shown to 


have severe consequences for MEM performance, particularly when the steering 
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vectors that represent the various beams of interest approach orthogonality; such is 
the case for the low order modes that are of interest here. If it is desired to use the 
MLM algorithm when highly coherent signals are simultaneously present on 
different beams, then a multiple beam approach that includes all possible beams 
upon which the coherent energy can arrive must be used. Alternatively. it may be 
possible to create an effective steering vector from a properly weighted sum of the 
beams. The former approach is generally infeasible. because the number of beams 
that must be included almost always exceeds the number of sensors available. so 
that the MLM inversion becomes singular. In the rare cases when an adequate 
number of sensors is available, one can still expect to encounter problems with the 
processing stability of the algorithm. On the other hand, the latter approach 
generally requires one either to index the beamforming over an unacceptably large 
number of parameters (all possible relative modal altes and phases), or to 
somehow solve the modal propagation problem a priori, such as through the use of 
normal mode or parabolic equation predictions. Although it is clearly not of ise for 


the present effort, this final method does offer some promise. 


It is interesting to note that in situations where the MLM algorithm has been 
traditionally emploved, namely, for horizontal arrays when high resolution is 
needed, the effects of the coherent interference problem are minimized. Higher 
resolution than that of conventional approaches is in. fact achieved, at the expense 
of acenracv m the absolute levels computed. This effect has probably been 
incorrectly attributed to the issue of MLM bias ın the past. maskmg the coherent 


interference problem. 


By comparison. the performance of the teast squares algorithms ts not sensitive 
to the issue of coherence between signal and noise. making it much more robust. 


Additionally, because the technique ds hnear, И does not suffer from the bras issues 
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known to plague MLM algorıthms. Still, the results of Chapter 5 indicate that least 
squares approaches provide sıgnificantly less resolution than MLM techniques when 


the latter can be applied. 


The major difference between multiple beam algorithms and their single beam 
equivalents stepped over the same set of steering vectors is that the multiple beam 
processors force the decoupling of the various beams while the single beam variants 
do not. Stated another way, each beam included in a multiple beam processor has 
its nulls arranged in such a fashion to prohibit the leakage of energy detected on the 
other beams included in the processor. Whether or not this null placement 
improves performance over that of the equivalent single beam processor depends on 
two things. First, the way in which the degree of freedom that each specially placed 
null represents is employed in the single beam approach must be considered. For 
the least squares algorithm, the single beam processor does not use null placement 
to improve the results; thus, the multiple beam least squares processor can be 
expected to provide at least marginally better performance by eliminating the cross- 
talk between beams. In the MLM algorithm, however, the single beam) processor 
deploys the nulls in such a fashion as to minimize the total noise power allowed into 
the beam. Thus, the multiple beam MLN > processor decouples the included beams 
only at the expense of less satisfactory noise rejection performance. If the energy in 
the different beams is incoherent, this is always a losing proposition in terms of total 
performance, as is shown in Appendix B. On the other hand, if the energy is 
coherent across different beams, then the single beam MLM processor experiences 
serious difficulties due to the coherent interference, and. the multiple beam MLM 


approach must be used. 


The other thing that must be considered. when discussing s mele beam versus 


multiple beam approaches is the processing stability issue. Tf the different beams 


| 
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that one is trying to decouple with the multiple beam algorithm look too much 
alike, then the decoupling process becomes singular; the beamformer must work so 
hard to separate the closely spaced beams that 1t becomes sensitive to even small 
amounts of background noise. The conclusion made here is that the issue can be 
quantified for the least squares algorithm through the use of the singularity 
coefficient defined in Section 5.10. Significant processing stability problems begin to 
occur at singularity coefficient values of roughly 0.4. For modal beamforming with 
vertical arrays, the problem arises from the fundamental inability of the array to 
resolve higher order modes. Two causes are possible: insufficient array length. whieh 
causes adjacent modes to look too much alike; and insufficient inter-sensor spacing, 
which causes Nyquist-like aliasing between modes widely spaced in mode number. 
The FRAM IV vertical array has been shown to be array length limited for the low 


frequencies of interest here. 


The final conclusion of this section deals with the usefulness of the restdual 
error as a performance measure for the least squares algorithm. Since the residual 
error essentially provides a goodness of fit measurement. rt can. If used carefully. be 
successfully employed to estimate critical unknown parameters from the acoustic 
data. One example is tts use in finding the array. tilt angle. Its reltability та 
accomplishing such a task is fundamentally related to the percentage of signal 
energy In the modes resolvable to the vertical array. When a large proportion of 
the stgnal energy resides in fundamentally unresolvable modes, then the method can 
be expected to work poorly, since the implicit signal model is inadequate. “This 
effect has been demonstrated with real data, Where reliable estimates of array. tilt 


for signals over 40 Hz cannot be made. 





8.2.2 Conclusions Concerning the Practical Use of Vertical Arrays 


Three major conclusions deal with the practical employment of vertical arrays 
in a real ocean environment. The first, which has already been mentioned, Is the 
modal resolution issue. The techniques developed in Chapter 5 are more than 
adequate for investigating this question for arbitrary vertical arravs, and can even 
be inverted for use in a design role, such as selecting sensor depths to maximize the 


singularity coefficient. 


From the comparisons made in Sections 3.5, 3.6, and 3.7, ıt can be concluded 
that the proper design and implementation of vertical arrays is a much harder 
process than the equivalent development of horizontal arrays. — Variations in 
hydrophone sensitivity across the array must be very carefully controlled, since 
typical applications of vertical arrays depend more heavily on accurate relative 
signal levels between different sensors than do those of horizontal arrays. The 
greater scopes of vertical arrays significantly increase the effect of ocean currents on 
array shape and, therefore, on sensor location, to the point where the net sensor 
displacements can no longer be neglected im the signal processing. This increased 
importance is due both to larger physical sensor Offsets and to the fact that these 


offsets now occur in an organized pattern across the face of the array. 


The results of Chapter 6 and Chapter 7 leave little doubt about the conclusion 
tliat modal beamforming techniques are extremely. sensitive to. variations in array 
tilt. This sensitivity increases for higher order modes and with higher operating 
frequency. The combination of this conclusion with the previous one shows tlrat 
improved methods of either controlling or measuring vertical array shape. and 
particularly the effective tilt angle of the array, are needed for proper exploitation 
of the imique capabilities of vertical arravs. From Section 6.4. it can be concluded 


that the effective array tilt must be known to better than about +0.25 in the 





Arctic Ocean if one wishes to accurately measure the amplitudes and phases of just 
the lowest order modes in the 0 to 80 Hz frequency range. A simple method of 
predicting the array tilt sensitivity has been presented there. It is important to note 
that the results of Chapter 3 and Chapter 7 indicate the shortest time scales for 
significant tilt angle variations in the FRAM IV vertical array are between 10 and 


20 minutes. 


8.2.3 Underwater Acoustics Conclusions 


The results of Chapter 7 paint a fresh and unique picture of low frequency 
sound propagation im the central Arctic, from which a number of important 


conclusions сап be drawn. 


The validity of all the other results obtained tn this area is dependent upon 
the fundamental conclusion that the modal structure developed in Chapter 4 Is a 
reasonably -accurate representation of that present during the FRAM IV 
Experiment. This inference can be readily drawn from Figure 7-7. since. in general, 
modal phase speeds tend to be much more sensitive to channel variations than mode 


shapes. 


The next significant conclusion to be drawn concerns the importance of the 
Arctic surface duct in sound propagation, even at the lowest frequeneres. This duct 
influences the FRAM IV sound propagation problem in three ways: by. controlling 
the modal distribution of source energy; by controlling the sensitivity of the 
horizontal array to the various modes; and by concentrating large amounts of sound 
energy very near the rough тсе canopy. It is obvious that the overwhelming 
majority of sound energy generated by апу source located in the surface duct is 
channeled mto modes that are effectively trapped in that duct. A stmiiar conclusion 


ix true for horizontal receiving arrays deployed within the surface duet. Тог the 
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frequencies of interest here, this means that the first mode has a tremendous 
relative advantage over all other modes; first mode source excitation is often more 
than 15 dB above that for other modes. It is also interesting and important to note 
the frequencies at which various modes transition into the surface duct. The first 
mode crosses over from the deeper Arctic profile to the surface duct at about 20 Hz. 
the second mode at about 40 Hz, and the third at roughly 65 Hz, although thev are 
not completely trapped in the upper layer until frequencies of about 30 Hz, 60 Hz. 


and 90 Hz, respectively. 


Given the great relative advantage that the first mode possesses compared to 
the higher order modes, the observation of only very small amounts of first mode 
energy at frequencies above 40 Hz is extremely important. It highlights the 
importance of modal attenuation effects in central Arctic Ocean propagation. In 
terms of what 15 observed from horizontal array data, two very different 
propagation regimes exist. One of these regimes occurs below 40 Hz, where strong 
first mode arrivals are observed on the horizontal array. The other regime. above 
40 Hz, shows a much more equal partition of energy in the higher order modes, with 
the deepest RSR paths being marginally dominant. Comparison with vertical array 
results explains this sharp change. At frequencies below 40 Hz, the increase in 
attenuation of the first mode with frequency is approximately compensated by 
increasing source excitation and horizontal array sensitivity. Only above 40 117. 
when these effects no longer mask it, does the strong frequency dependence of the 
first mode attenuation become apparent in the horizontal array data. The 10 117 
breakpomt between the two regimes can be considered to be an effective upper 
bound to the range of frequencies that. propagate efficiently in. the Aretic surface 
duet, and may thus be taken as a high frequeney limit for long range acoustic 


svstems designed to be employed m the Arctic environment. 
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The modal attenuation measurements made for the first three modes from the 
vertical array data verify the dominant role that the rough surface of the ice canopy 
plays in Arctic Ocean sound attenuation. The frequency dependence of the 
attenuation for the different modes correlates well with the frequencies at which 
they transition into the surface duct, suggesting that it is the ice cover that provides 
the attenuation mechanism. Moreover, the agreement of the measured modal 
attenuation with predictions using the full method of small perturbations and 
concurrent experimental estimates of the ice statistics Is excellent for a model 
consisting solelv of a rough free surface at the ice-water interface. This result is at 
variance with that of DiNapoh and Mellen [25], in which they cannot reach 
adequate agreement between experiment and theory using measured values of the 
ice statistics. The reason for the difference mn results is not presently understood. 
The Kirchhoff approximation is shown to be poor by comparison in this situation. 
The result demonstrates the need for the careful inclusion of surface correlation 
effects in rough surface scattering models. Additionally, it is possible to interpret 
the second mode attenuation results as indicating that multiple scattering effeets, 
or, equivalently, modal coupling effects, are significant in the rough surface 
scattering problem posed by the ice canopy. Certainly, a great deal more work i5 
needed to verify this interpretation; if it is correct, however, it demonstrates a 
unique case where the modal coupling is significant even though the coupling 
coefficients themselves are negligible. What would make the coupling important iu 
this case is the great difference in modal excitation levels. Even very low level 
scattering from the first mode into the second or third modes may be пороге 
when the original excitation levels of these modes are 20 or 30 dB smaller than that 


of the first mode. 


The signal stability observed. throughout. this investigation has been 





remarkable. Many of the temporal variations observed can be attributed either to 
background noise effects or to vertical array movement, although assignment of 
observed fluctuations to the latter cause cannot be done conclusively. This is due to 
the lack of concurrent vertical array tilt measurements for the experiment. The 
only exceptions to these findings are some signal fades observed only on selected 
hydrophones of the arrays used; the cause of these variations is not presently well 
understood. The observed stabihty of the central Arctic sound channel and the 
high values of the modal coherence measured strongly support the conclusion that, 
for the low frequencies investigated here, the central Arctic channel сап be 
considered to be completely deterministic m a temporal sense for periods in excess 
of one hour. These results agree closely with those of Mikhalevsky [57], ђиг ако 
break new ground, since the various propagation paths have now been at least 
partially separated for independent inspection. Not only are the relative phases of 
the various modes very stable; it appears that the environment causes almost no 
phase fluctuations of any kind. Thus, the different propagation paths appear to be 
completely coherent, not so much because they all exhibit the same phase 
fluctuations, but more because none of them exhibit any significant phase 
fluctuations at all. This observation has a number of important implications. First. 
it means that any energy scattered due to rough surface effects (particularly from 
the ice) remains coherent with respect to the transmitted field. This is why sound 
transmissions above 40 Hz. while obviously suffering major amounts of scattering bv 
the ice, still display such stable phase traces at the different sensors. Second, it 
means that totally unrelated tonal sources can still exhibit very high cross-coherence 
in the central Arctie: the primary hmiting factor ıs the ability to mateh the different 
sources In frequency. This, dm torn. makes the use of high resolution beamforming 


algorithms that are subject to coherent interference (such. as the. MEM algorithm) 
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highly suspect for stationary (non-moving) tonal sources in the central Arctic, unless 
either only a single source is known to be present, or the averaging length 15 
extended to periods of at least several hours, if not several days. Most importantly, 
the observations show that the channel does not destroy the relative phase 
information from which target range information might be extracted. Therefore. 
target range estimation Is feasible with respect to the temporal fluctuations in the 
central Arctic. What is not so clear from the results of this investigation Is whether 
relative phase predictions can be made accurately enough to support such a passive 
ranging effort, especially if multiple scattering or mode coupling effects are. in fact. 


significant. 


In retrospect, It Is interesting to apply the insight gained from the Arctic data 
presented here back to the open oceans of the world. Clearly, these warmer waters 
do not exhibit anything near the same stability and coherence as observed in the 
central Arctic. Since the most glaring difference between the two environments Is 
the type of free surface involved, one must at least suspect that the temporal 
instabilities more regularly observed in oceanic sound transmission are intimately 
related to the time variations in this free surface from wind and wave action. 
Certainly, these effects generate the strongest temporal variations im channel. The 
Arctic results prove that temporal channel variations are needed to generate 


temporal fluctuations in the propagating sound field. 


8.3 Some Thoughts on Further Work 


Every scientific Investigation answers certain questions while raising other, 
previously unseen, Issues. This thesis is certatnly no exception. Before concluding, 
then, Ц 1% useful to discuss at least a few areas that may prove fruitful for further 


research. 
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From a signal processing Viewpoint, two valuable investigations come to mind 
immediately. The first is the issue of coherent interference in the MLM algorithm. 
For a situation such as that found in the central Arctic Ocean, no high resolution 
beamforming alternative that presently exists could be expected to work well. 
because of the coherent nature of the channel. An investigation of how one might 
successfully employ high resolution concepts in this type of environment would be of 
great interest, not only for its practical value, but also because it would probably 
provide a great deal of insight into how present high resolution algorithms nught be 
made more reliable. It ts this lack of robustness which often inhibits potential users 
from emploving high resolution techniqnes. — Additionally, such. an. investigation 
would probably remove much of the mystery surrounding the often misunderstood 


MLM bias issue. 


A second useful investigation would deal with better methods of inferring 
vertical array tilt from a received acoustic field. If a highly aceurate method of 
accomplishing this task can be found (one much more reliable than the technique 
used here), it would provide a rather cheap and simply implemented solution to the 
most important problem facing anyone Wishing to use vertical arrays for modal 
decomposition. Ву comparison, any other type of solution, such as an in situ 
measurement scheme or a modified deployment method. is likely to be both 
cumbersome and expensive. A related effort to this would be an mvestigation of the 
effects of higher order array shape variations on modal beamformmeg. Although 
assumed negligible here, there is likely to be some frequency /mode number linut 
bevond which even the small amounts hkely to be encountered m= practice are 
significant. Such an Investigation rm. fairly straightforward, and. would. add 
confidence to. the beltef. that. the shape of a vertical array сап be adequately 


modeled by a straight line at some angle to the vertical. 
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The most interesting effort that could be undertaken in the underwater 
acoustics domain to follow the present work would be a repeat of this study in an 
open ocean situation. Lacking that, two other related investigations come to mind. 
First, the indications that mode coupling might be significant for the rough surface 
scattering observed need to be studied more closely. Intimately coupled to that 
question is the more difficult issue of whether or not relative modal phases can be 
predicted from the general characteristics of the medium with enough accuracy to 
allow the source range to be estimated. Such an investigation would obviously have 


to include an assessment of the effects of mode coupling. 


Finally, it would be both useful and enjoyable to repeat the exact same 
experiment with an adequate array tilt measurement system m place. The lack of 
tilt angle measurements has caused more than a little frustration m this thesis, and 
In too many ways the results are still tenative as a result. In particular, it would be 
interesting to know just how effective the array tilt estimation scheme used here Is 
In practice, since this is probably the greatest remaining unknown affecting. the 


results presented here. 
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Appendix A 


On the Inclusion of Signal 
Components in Generalized Least 
Squares Weighting Matrices 


Consider a generalized least squares estimation problem with an Ух ~ 
Hermitian positive definite weighting matrix УУ, апі ап JN X M complex steering 
matrix E. The M x 1 complex vector À, which minimizes the scalar error Q 


=) 


defined in equation (5.14) is given in equation (5.15), which is repeated for 


convenience 
A Sse ЕЕ ЕЕ lot Ее 
A Е ЕМЕ (4.1) 


Here, P is the observed complex pressure vector. 


Typically, the desired weighting matrix is the inverse of the sensor cross- 
coherence matrix for the noise. However, it is usually more convenient to use the 
total sensor cross-coherence matrix, Including both the signal and noise components. 
since estimation of this matrix does not require separating the signal from tlie nore. 
Thus, one must be concerned with the effect of including signal components in the 
inverse of the weighting matrix. Let W be such a modified weighting matrix. Its 


form can be defined by 


W!=W-!+EBB+E*, E 


where the .M Xx 1 complex vector B consists of arbitrary constants. The second 
term of equation (4.2) represents the inclusion of the stgnal components. Using the 


Identity 





(Z+UVt)} =z! = ——_—_ КЕ 


NE NN 

— — W,EBBtEtW, 

Y= I wen ЕШ 

1+В Е W,EB 
The optimal estimator for W 15 

А —(Е' МЕЈ ТЕ МР. (1,5) 
Using equation (A.4), it is possible to write 

___ _____ (E*W,E)BB*(E* W. E) 

E'WE —(E' W, E) DE ш эы (1.6) 

1+В ЕУ ЕВ 

which can again be inverted using equation (.3.3) to yield 

(E'WE) —(E'W,E) -BB'. (17) 
Substituting eqnations (A.4) and (A.7) into equation (1.5) vields 

5 кезу БОБА, Р— 4 : 

I EWR ENGE IN (1.8) 


It can be seen that the least squares estimator produces the same result for 
either of the two weighting matrices. If W, Is taken to be the inverse of the noise 
portion of the sensor coherence matrix, then, from equation (4.2), W Is the mverse 
of the full sensor coherence matrix (including both signal and поке) {ога созрее 
coherent signal that Is uncorrelated with the noise. Note that this is the nature of 


the requirement for independent signal and noie im the MEVE algorithm. 


The above result can easily be generalized to the case of à randon signal that 
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is less than fully coherent. In this situation, the relation between the two weighting 


matrices is of the form 


ЕЕЕ ЕЕ Б Б 

ХМ => +ЕЗ Е’, (9) 
where Ss is the M X M modal (or signal) cross-coherence matrix | 

= _ A | 

ВВ. (-1.10) 


This is a complete generalization, since for complete signal coherence the vector B 
becomes non-random, and the result then reverts back to that given in equation 


(A.2). 


By singular value decomposition, the matrix S,, can be written as 


M | 
— u a | 
Su = 5 A; е; S Е (A. | | ) 


=> 


where the A, and e, are its eigenvalues and eigenvectors, respectively. The more 
general result is then proven by applying the approach presented above M times In 
succession, once for each eigenvector in the sum. It is important to observe that 
while this generalization relaxes the deterministic signal assumption, It in no wav 


alters the requirement that the signal and noise be completely uncorrelated, 
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Appendix B 
The Performance Relationship Between Single 
Beam and Multiple Beam MLM Algorithms 


In the first 


Consider two related multiple beam MLM estimation problems. 
From equation 


problem, M modes (or beams) are included in the beamformer. 


(5.20), the estimate of the M x Af modal (or beam) cross-coherence matrix for this 


situation Is 
(13.1) 


Ел 
Sy = (Ey, Sx Е\ . 


Where E q ìs ће № х Af steering matrix and Sv is a sensor cross-coherence matrix 


of size У Х N. 
As the second problem, consider a multiple beam beamformer having a total 
and let the first M of these modes be identical to those 


of M + 1 modes included 

from above. The steering matrix for this problem is then № х A/41, and may be 

written as 
(1.3) 


Еда = | Еп | Е]. 


where E is the .V X I steering vector of the additional mode. The estimate of tli 






MEL XEM modal ¢ross-coherence matrix is 
гс-ір i Rte! -1 

E Sy E| EuS E 
Со) 





р = 
ETS E, | ET Sy E 


EN. 


ESTEE! 
zx pL (Ey, SEES 


Np 


An estimate of the modal eross-coherenee matrix for the original M modes 


may be obtained from the second problem by extracting the upper left. hand 


ЕХ 





атр 


(BA) 


—M+] 





= | 
S,, may then be identified by using the formula for the inverse of a partitioned 


Hermetian matrix. If 














T Y|Z 
x = | (1.5) 
Ztl v” 
and 
" A|B 
x! = F (B.6) 
Bt с“ 
then 
2 ¡EE ib 
= a = E 
и 
В Y'z (BR) 
Пи А | 
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: (8,9) 
= om. Jf 
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From equation (B.7), 
SES Eya + (B.10) 
PRESS EE SE MESE," 
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Equations (B.1) and (B.10) can be interpreted to show that the beamformer 
containing Af modes always provides better performance than the one containing 
M -- 1 modes. Consider a situation where the sensor cross-coherence matrix 
consists solely of noise with no signal. For these conditions, the two different 
estimates of the M X M modal cross-coherence matrix indicate the amount of noise 
that is not rejected in each beamformer. In particular, the diagonal terms indicate 
the noise power that leaks into the amplitude estimate of each mode. A comparison 
of equation (B.1) with (B.10) shows that the diagonal terms in the latter are never 
less than the diagonal terms in the former. The additional contributions indicated 
in (B.10) are always positive. since the numerator of the fraction is a matrix formed 
bv taking an outer product of a vector with itself, while the denominator ts a scalar 
that is also guaranteed to be positive. This last conclusion comes from a 


о Е at? : : pem M Ж 
consideration of &3,, , in equation (B.4). If one assumes S.. ta be positive definite, 
4 А 


+1 
39 


then it follows that S must also be positive definite and, hence, c must be 


—\1+1 44-1 
positive. But, from equation (B.9), 
= (BAL) 
1 
+ с—1 + <-In т-с-ір \ı-1l/7+ o-Im\' 
(E S\ E) = (E 5. I (Ex SN E, y) (Ex, SN E) 


so that the denominator (which is identical to that in. equation. (B.10]) must be 


positive. 


More insight can be gained by considering a situation where the sensor cross- 
coherence matrix is a combination of two components, one consisting of noise and 
the other being a signal of amplitude a in the A/F PS mode. From the viewpoint of 
the first M modes, the signal in the M+ mode represents additional noise against 


which the processors must discriminate. For this situation, the sensor cross- 
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coherence matrix Is 


$1 = 5, + ||" ЕЕ*. (В.12) 


After much algebra, it can be shown that 


2 (B.13) 


(Е 5, En)” (ENS, ' E) (ES, ' Es) (Es, S,' Ey)” 
1/|е| + (Е* 6 РЕ) —(Е*5 Е Ј(Е 5" Е, ) (Е 5" Е) 
while 
S — —_.л— = 
Su = (Es S, | E i "E (Ш.Е 
(EX $,' Ej)! (Ex S; ' EJ (E* $,' EA) (Es S, EXE 
(E'S.!E) -(E* S! E.) (EAS 5 Е) | (Bt s-!E) 


Note that the result in equation (B.14) is completely independent of the amplitude 
of the M+1% mode, and is in fact identical in form to equation (B.10). This is 
indicative of the fact that the beamformer containing M + 1 modes automatically 
generates a null in the direction of the A/+1” mode for the beam patterns of all the 
other modes, whether or not any signal is actually present there. The result ın 
equation (13.13), on the other hand, balances the amount of noise against the 
amount of A/+1”" mode signal, optimally deploying the extra null to best 
discriminate against the combination of the two. When the amplitude goes to zero, 
the second term in equation (B.13) disappears entirely, and the beamformer 
completely ignores the M-- D mode. Conversely, as the amplitude gets larger, the 
beamformer increasingly diseriminates against the mode: at very large amplitudes. it 
Is forced to devote a full null to the mode, and result approaches that of equation 


(13.14). 
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If a multiple beam MLM beamformer containing Af modes always produces 
better performance than one containing M -- 1 modes, then à .M— 1 mode 
beamformer always performs better than one of M modes. "Thus, it is obvious that 
a single beam MLM beamformer (which contains only one mode) always provides 
better performance than ary multiple beam MLM beamformer. However, this 
statement must be qualified by the assumptions inherent in the derivation of the 
MLM algorithm, particularly the assumption that the signal and noise be 
independent of each other. For a single beam MLM beamformer operating in a 
sound field consisting of many modes, this assumption is met only if the different 
modes are incoherent with respect to each other. If this is not the case, then 
coherent interference adversely affects the performance of the single beam 


algorithm, a degradation not accounted for in the analysis presented here. 
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